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Introduction 

Magnetic  resonance  is  normally  very  much  broadened  if  a  ferrimagnetic 
sample  is  immersed  in  a  spatially  non-uniform  field.  Experimentalists  measur¬ 
ing  fundamental  resonance  parameters  take  great  pains  therefore  to  employ 
ell  ipsoidal  sample  shapes  (usually  small  spheres)  that  are  positioned  in  fields 
of  very  high  uniformity.  Because  surface  roughness  is  known  to  cause  scatter¬ 
ing  from  the  uniform  mode  to  degenerate  spinwaves  of  short  wavelengths,  there¬ 
by  increasing  the  resonance  linewidth,  additional  effort  is  expended  in  polish¬ 
ing  the  surfaces  to  optical  tolerances. 

Commercial  manufacturers  of  tunable  microwave,  yttrium  iron  garnets  (YIG) 
filters  avail  themselves  of  this  knowledge  and  employ  uniformly  magnetized, 
highly  polished  spherical  single  crystals  in  their  designs. 

From  this  perspective  it  is  therefore  remarkable  that  we  at  MIT  observed 
extremely  sharp  resonances  of  a  very  localized  character  in  single  crystal  YIG 
slabs  and  films  that  encounter  highly  uniform  bias  fields.  On  the  other  hand, 
it  has  been  known  for  some  time  that  magnetoelastic  waves  can  be  highly  focussed 
by,  and  propagate  with  low  loss  in,  steep  magnetic  field  gradients.  One  view  of 
the  high  Q  resonance  is  that  magnetostatic  mode  patterns  are  formed  for  which 
the  resonant  energies  are  highly  confined  to  certain  regions  or  "tracks"  wi thin 
the  crystal  that  allow  wave  propagation  around  them.  If  the  mode  amplitudes 
are  very  small  at  the  edges  and  corners  of  the  sample,  the  surface  scattering 
(which  one  would  expect  to  be  enormous)  is  largely  prevented;  consequently  the  Q 
of  the  resonance  is  governed  primarily  by  the  intrinsic  linewidth  of  the  bulk 
crystal  together  with  normal  circuit  loading  considerations.  In  effect, 
appropriately  designed  magnetic  field  profiles  create  surfaces  of  discontinuity 
where  there  are  no  actual  surfaces;  surface  wave  propagation  at  such  "surfaces" 
should  be  free  of  many  of  the  drawbacks  and  loss  mechanisms  encountered  at  true 
surfaces.  It  has  been  our  intent  to  learn  how  to  characterize,  control  and 
efficiently  couple  to  such  modes  so  that  one  can  create  a  new  class  of  microwave 
magnetically-tunable  resonance  filters. 

Our  research  goals  concerned  Magnetoelastic  Delay  Line  and  Magnetostatic 
Mode/Wave  Synthesis;  we  separately  enumerate  our  results  for  each  of  these  two 
major  topics. 
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MAGNETIC  FIELD  SYNTHESIS  PROCEDURES 
FOR  MAGNETOSTATIC  AND  MANETOELASTIC  DEVICES* 


F.R.  Morgenthaler  and  A.  Placzker 


Department  of  Electrical  Engineering  and  Computer  Science, 
Massachusetts  Institute  of  Technology,  Cambridge,  Mass.  02139  and 
Chu  Associates,  Littleton,  Mass.  01460 


ABSTRACT 

Certain  types  of  magnetostatic  and/or  magneto- 
elasclc  devices  require  nonuniform  dc  bias  fields 
of  sufficient  strength  to  locally  saturate  the 
active  ferrite  element. 

We  here  review  a  synthesis  procedure  for  cy- 
lindrically  symmetric  geometries  that  allows  pre- 
specification  of  the  field  either  on  the  symmetry 
axis  or  on  a  plane  perpendicular  to  it. 

The  method  is  then  applied  to  the  cases  of 
both  a  thin  film  disk  magnetized  normally  to  its 
plane  and  microwave  magnetoelastic  delay  line 
designed  for  linear  frequency  dispersion  over 
wide  bandwidths. 

Finally,  we  report  construction  details  and 
test  data  on  an  actual  packaged  device  having 
greatly  Improved  characteristics.  Measured  para¬ 
meters  for  two-port  operation  at  S-band  over  a 
1  GHz.  bandwidth  include  a  linear  dispersion  fac¬ 
tor  D-.3  nsec/MHz.  and  an  insertion  loss  (untuned) 
of  27-30  Db. 


INTRODUCTION 

Although  many  magnetic  resonance  devices  ideal¬ 
ly  operate  with  a  dc  magnetic  field  bias  that  is 
spatially  uniform,  there  are  a  number  of  magneto¬ 
static  and  magnetoelastlc  wave  devices  that  either 
require  or  benefit  from  nonuniform  dc  fields.  Al¬ 
though  the  equations  governing  the  dc  magnetic 
field  within  a  locally  saturated  ferrimagnet  are 
well  known  and  can  be  solved  numerically  for  any 
set  of  boundary  conditions  that  force  uniqueness, 
the  solution  when  found  may  not  be  the  field  re¬ 
quired  for  proper  device  operation.  Fortunately, 
a  synthesis  procedure  has  been  formulated1-  that  al¬ 
lows  one  to  specify  the  desired  field  along  axes  or 
planes  of  symmetry  and  then  work  outward  to  find, 
first,  the  entire  field  within  the  magnetic  mater¬ 
ial,  second,  appropriate  high  permeability  pole 
pieces  or  current  windings  capable  of  sustaining 
chat  field. 

We  first  describe  the  synthesis  procedure  and 
show  how  to  produce  prescribed  nonuniform  radially 
symmetric  dc  fields  in  a  normally  magnetized  thin 
film.  We  then  apply  the  method  to  high  perfor¬ 
mance  magnetoelastic  linearly  dispersive  delay 
lines  that  employ  cylindrical  crystals  of  yttrium 
iron  garnet  (Y1G) . 

SASIC  EQUATIONS 

Within  a  currenc  free  magnetic  material,  and 


neglecting  crystalline  anisotropy,  the  dc  H-field 
is  both  curl  free  and  parallel  to  the  magnetiza¬ 
tion  vector,  M.  If  the  latter  is  assumed  locally 
saturated  to  a  uniform  value  M,  it  follows  in  terms 
of  the  scalar  magnetic  potential  ■.</ 

Chat  H  -  Vip  (1) 

atld  M  -  MV  ip/ |  Vi|i|  (2) 

The  final  constraint  is  (H+M)  »  0,  there¬ 
fore  ° 

v-[(i-m/|v<i»  1)7^]  -  o  (3) 

Although  this  equation  can  be  solved  numeric¬ 
ally  subject  to  any  sec  of  boundary  conditions 
that  force  a  unique  solution,  the  analysis  of  such 
a  boundary  value  problem  is  often  tedious.  More¬ 
over,  and  of  greater  importance,  the  H-field,  when 
found,  may  not  be  that  which  is  desired.  Then, 

Che  boundaries  must  be  modified,  the  field  recom¬ 
puted  and  so  on  in  an  iterative  "cut  and  try"  man¬ 
ner.  It  would  appear  to  be  much  more  difficult  to 
carry  out  field  synthesis  rather  than  analysis. 
However,  the  synthesis  of  a  desired  field  is  actu¬ 
ally  a  much  simpler  problem. 

In  essence,  the  procedure  starts  by  assuming 
the  specified  field  to  exist  along  an  axis  or 
plane  of  symmetry  for  which  the  direction  of  M 
can  be  deduced.  The  potential  p is  then  expanded 
in  an  appropriate  series  within  Che  magnetic  re¬ 
gion.  This  portion  of  the  synthesis  is  termed  the 
"lnner"-f leld  determination  and  is  carried  out 
without  regard  to  the  boundaries  of  the  magnetic 
material.  The  next  stage  of  the  procedure  Involves 
expanding,  in  a  convenient  series,  the  Laplaclan 
potential  in  the  nonmagnetic  region  outside  of 
some  assumed  boundary  surface.  The  coefficients 
of  the  "outer"-f ield  potential  are  next  matched  to 
those  of  the  "inner"-field  so  as  to  satisfy,  in  a 
least-squares  sense,  the  proper  boundary-conditions 

Finally,  we  design  high-permeability  pole- 
pieces  or  current  sheets  that  when  suitably  ener¬ 
gized  will  create  the  outer  and  hence  also  the  in¬ 
ner  field. 

AN  OVERVIEW  OF  THE  SYNTHESIS  PROCEDURE 

Although  more  general  situations  also  can  be 
treated  by  similar  methods,  ve  here  restrict  our 
discussion  to  cylindrically  symmetric,  fields  when 
the  material  shape  is  a  cylinder  or  disk.  There 
are  then  two  cases  of  particular  interest.  They 
result  from  specification  of  the  H-field  along 
either  the  axis  of  symmetry  H  (o,z)  or  an  appro¬ 
priate  plane  perpendicular  to  that  axis  (H^tr.z^). 
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The  synthesis  procedure  can  be  divided  Into  the 
following  steps: 

1.  Find  the  H-field  that  Beets  Che  on-axis  or  on- 
plane  field  requirement  Inside  of  the  cylinder 
of  lengch  L  and  radius  R  for  specified  magne¬ 
tization  M. 

2.  Choose  a  convenient  series  expansion  for  the 
Laplacian  scalar  potential  outside  of  the  cy¬ 
linder  and  match  boundary  conditions  over  the 
entire  surface.  Solve  for  the  coefficients  of 
Che  outer  potential.  It  may  prove  useful  to 
subdivide  the  outer  region  and  carry  out  a  se¬ 
parate  a  match  for  each  subregion. 

3.  Plot  the  equlpotentlals  of  the  outer  field  and 
choose  two  (or  more)  that  are  appropriate  to 
serve  as  surface  contours  of  high  permeability 
magnetic  pole  pieces.  If  subdivision  of  the 
outer  region  has  been  utilized,  pole  pieces 
must  in  general  couch  Che  macerlal  ac  Che  divi¬ 
sion  poincs  so  as  to  isolate  the  various  outer 
subregions. 

4.  Alternatively,  if  a  solenoid  is  to  be  used  to 
generate  the  field  without  benefit  of  magnetic 
pole  pieces,  Che  field  outside  of  Che  winding 
surface  must  be  chosen  with  continuous  normal 
flux  and  vanish  as  r-*®  .  The  discontinuity  in 
che  tangential  field  then  determines  the  sur¬ 
face  electric  current  density  and  hence  wind¬ 
ing  design. 

5.  If  a  satisfactory  design  does  not  result  from 
the  synthesis  procedure,  one  can  alter  the  as¬ 
sumed  material  boundary  surface  and/or  the 
"outer"-field  expansion  and  try  again. 

It  is  important  to  realize  chat  although 
an  infinite  number  of  combinations  of  material 
size,  shape  and  pole  piece  designs  exist,  all 
of  which  would  create  the  desired  field,  once 
one  has  been  chosen  the  field  it  creates  is 
unique. 1 

INKER  FIELD 

Because  symmetry  dictates  that  both  H^o.z) 
and  M^Co.z)  must  vanish,  and  the  axis  (r-0)  is  as¬ 
sumed  nonsingular,  an  appropriate  expansion  of  che 
potential  within  the  magnetic  macerlal  is 
r  ,  .  2n 

*  "  n-0a^(z)  r  (4) 

One  could  substitute  Eq.  (4)  into  Eq.  (3)  and  ex¬ 
pand  che  result  so  as  to  find  a-  in  terms  of  a 

2  o 

and  its  derivatives  and  so  on.  However,  it  is 
more  convenient  to  separately  expand  M  as 

5"MfIz  nI.Ob!^r\^Ob2n>l(z>r2n+15  <-S> 

and  require  |M|  -  M,  MxH  •  0,  and  V-B  ■  0.  The 
result  is  three  sets  of  constraints,  respectively 

Il3'2n  b [2(s~n)+i  ] ~2(s+i-n) *2 (s+1  -n) b2n  ^  *  0  (ba) 

I^b2n+lb2(s-r.)+l  *b2ob2(s+l-n) 1+bob2 (s+1)  -  0  (6b) 

4(s+l)2a2(s+1)  +2(s+l)MbJs+1  +  *  0  (6c) 


where  the  primes  denote  differentiation  with  res¬ 
pect  to  z  and  s«0,  1,2,3,...  In  addition,  b~-  1 

0 

and  o  a’  >C.  In  order  to  ensure  local  saturation. 

is  either  positive  or  negative  definite  over 

the  incerval:  we  take  a'>0  and  b  «+l  without  109a 
o  o 

cf  generality.  Eqs.  (6)  yield,  for  each  value  of 

3:b  , ,  a  ,  and  b  ,,  it  terms  of  the  lower  order 

s+1  s+2  s+2 

a  and  b  coefficients  and  their  derivatives.  There¬ 
fore  one  starts  with  s-0  and  proceeds  upwards  to 
generate  as  many  terms  as  are  needed.  For  s»0 


-  a"/2(a  ’+M) 

0  0 

(7s) 

a'b-f/2 

0  1 

(7b) 

-bjr/2 

(7c) 

Provided  a^(z)  Is  either  specified  or  can  be  de¬ 
duced,  che  various  functions  a2n  constitute  the 

desired  solution  of  Eq.  (3).  The  number  of  terms 
required  to  satisfactorily  approximate  che  field 
depends  upon  the  excent  of  r,  the  value  of  M  and 

the  particular  function  a', 
o 

If  the  material  has  uniaxial  magnecic  anis- 
cropy  oriented  along  the  z-axis,  the  procedure  can 
be  generalized  by  replacing  in  Eq.  (6a) 

2K 

a2n  a2n~  u"m  b2n 
o 

where  Kq  Is  the  uniaxial  anisotropy  conscanc. 

(F.  <0  easy  axis;  R>0  ea3y  plane).  If  the  anis- 
o  o 

tropy  is  not  uniaxial  with  respect  to  the  z-axis, 
(as  for  example  a  cubic  material  with  [100]  or 
[ill]  orientation)  the  formulation  may  still  be 

used  by  replacing  Kq  with  an  appropriate  effective 

value  provided  the  radial  component  of  M  is  not 
too  large.  If  M^/M^  is  large,  the  field  is,  of 

course,  not  strictly  cylindrically  symmetric. 

The  expansion  of  Eq.  (4)  is  very  helpful  when 
che  on-axis  field  is  what  is  specified  because 

then  a’(z)  is  known.  On  the  other  hand,  if 
o 

Hz(r,zQ)  is  specified^n  alternace  approach  Is  pre¬ 
ferable.  In  this  case  we  expand  ^  as 
p»a  (r)+a  (r)(z-z  )+a-(r)(z-z  )2+.  .  .  (8) 

0  1  0  4  o 

with  a,(r)*H  (r,z  )  and  assurance  chat  a'(0)-0 
a.  z  o  n 

for  all  n. 

THE  OUTER  FREE  SPACE  POTENTIAL 

If  the  z-axis  passes  through  an  outer^ field 
subregion,  the  Laplacian  outer-potential  for 
that  subregion  may  be  taken  nonsingular  over  all 
z  and  expanded  in  the  form 

-  „  a"(z)  ,  a"”(z)  , 

*  ‘  ao(z)-  mr2-  ¥  +  ¥~  •  •  •  (9) 

»  9  in 

Naturally,  if  is  taken  to  be  (cos)  (kz)  or 
( co  aft)  (kz)  'l>  factors  into  the  product  of  aQ 
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and  either  I  (kr>  or  J„(kr) .  However,  th*  usual 

cylinder  function*  are  not  especially  convenient 
beceuse  our  boundary  specification  do  a  a  sot  1m  d 
readily  to  iedentificacion  of  eigenvalue*  of  k. 

For  *0“*ni  Eq.  (9)  generate*  eh«  eat  of  poly¬ 
nomial*  PJ(r,e)  that  satisfy  the  recurrence  formula 


?•-  dlSli  i  f*  -  Y  C.V)**  . 

n  n  n-1  n  n-^ 


With  PJ-l,  it  follow*  that  P*«(r1+*2)n^2?5*eo*0) 


wh*r*  tin?1  r/x  «nd  P®  i*  th*  ***oclat*d  L*g»ndr* 

function  of  d*gr««  n  and  ord*r  zero. 

For  subregions  of  th*  outer  field  that  do  not, 
contain  th*  x-axls,  eolutlon*  with  a  logarithmic 
elngularlty  at  r«*0  ar*  often  helpful.  In  euch 
ea*«e,  w*  «oploy  th*  **t  of  polynomial*  QJ(r,*) 

that  eatiefy  th*  recurrence  formula 


Q*« 


HT  *  «S-i  - 


n-1 

a 


th*  radial  match  at  r»R  and  th*  *nd  match  at  *“0 
or  z“l  eeparataly. 

Radial  hatch 


We  wish  to  expand  H  (z,R)  and 

(Hr+Mr)(z,R)  in 

power  series  form 

»,<*.«  *  S 

(11a) 

*nd  N-1  - 

(Ht+Mr)(*,R)  -  n:0  3b(z/L) 

(lib) 

where  N  1*  th*  nvszber  of  term*  that  will  produc* 
a  tolerable  error.  Although  an  and  ^  could  be 

found  from  matching  term*  of  a  Taylor  aerie*,  w* 
often  find  it  preferable  to  calculate  then  from 
a  laaet  squar**  fit,  uelng  the  invar**  Hilbert 
Matrix  of  order  S. 

In  term*  of  an  *«d  8a»  th*  coefficient*  and 

are  then  found  to  »ati»fy 

C  -a  ./  Cpl1”1)  -  (inR+l)C '  (11*) 

p  p-i  P 


with  Q*  •  (1  +  inr) 

It  1*  *l*o  permieeable  and  oftan  advantageou* 
to  utlllae  in  th*  expansion  axial  multipole*  of 
the  form 

ii-  ((*-*  )/[(*-*„  )  W/J)  n-0,1,2 . 

.  B  0  0 


aa  long  a*  their  location*  r-0,  z»*0  are  anywhere 


within  the  materiel  boundary. 

We  plan  to  ue*  iron  pole  piece*  to  energize  all 
of  th*  eyntheeized  field*  described  in  thi*  paper, 
therfor*  th*  behaviour  of  th*  polynomial*  at  large 
distance*,  from  th*  origin  i*  no  deterrent  to  ex¬ 
panding  t|/ln  any  outer  region  aa 


*■  V 


N  II"  1* 

nil  VS  +  n«0  VS 


(10a) 


or  what  1*  exactly  equivelent, if  C^-0 


ilm 


S  n*/2 

n«0  kA  IVVlor+V>An.k* 


n-2kf2k 


where 


and 

*k’ 


iiiliaL. 

<n-2ttl(kl>*4fc 

1 

+  1/2  +  1/3  +  ,.,.1/k) 


k-0 

K2l 


and  n*  i«  the  even  integer  n  or  n-1. 


(10b) 


BOtHtDART  COMOITIOKS 

Secauae  th*  **mpl*  *h*p*  1*  her*  r**trlct*d  to 
be  a  right  circular  cylinder  (or  diek)  of  radio*  S 
end  length  L,  it  ie  beneficial  to  *ubdivid*  the 
outer  region  at  th*  cornqr*  and  carry  out  three 
•operate  expanelone  for  tl>  .  Tho*e  for  *4  and 
*>L  cannot  contain  Q*  function*.  • 

Ptoa  th*  "inner"  field  *ynth**l*  w*  know  H 
end  M  everywhere  within  th*  cylinder .  W*  consider 


r21c  V*VVv*anR+Vc;V 


k-1 

C’-Rfl  /LP 
P  P 


lh  r2"  Vu!2i!c^ 


(12b) 


Notice  that  for  a  given  value  of  n,  and  de¬ 
pend  only  on  the  value*  of  Cp  and  with  ?>n. 
Therefore,  *inc*  by  design  C^-0,  Eq.  (12a)  immedia¬ 
tely  yield*  Cjj .  Thereafter,  alternating  between 

Eq*.  (12b)  and  (12a)  for  successively  decreasing 
value*  of  p,  produce*  an  unravelling  that  caieu- 
latae  in  order  C^[  C^_2,  C^_2,  C,,_2i  •  •  •  !  (,■ 

The  value  of  C  la  immaterial  and  can  be  eet  to 
0 

any  convenient  velu*  including  zero. 

End  Match 

In  a  like  manner,  th*  outer  potential*  for 
z<0  or  *>L  can  be  expanded  a* 

„  n*/2 

^end"  S  k«0 


z<0 


(13a) 


*‘ni  n-2k  2k 

*„d-  S  kio  cA.k<-“  1 


*EL 


(13b) 


For  the  z-0  end  face,  matching  4 and  H^+M^  gi\*e 


a  (05 

n 

Ln,a 


n  iv«tv 


(14) 


*n-l^+Kb„  ,(0)  n  odd 

_  U" 1  .  . 

*n,n-l 
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vhere  a  and  b  satisfy  Zaa.  (6).  Similar  aqua- 
n  Q 

tiona  result  for  cha  i*L  and  fac«. 

A a  an  example,  aaauma  chat  cha  on-axia  fiald 
in  a  cylindar  of  radiua  R* 1 . Jma  and  langch  L“5mo 
la  required  co  ba 

Hi(o,*)  -  300  +  300*  Oa. 

whan  cha  aaturatlon  magnetization  expressed  in 
gauss  ia  4ttm  -  1780C. 

A  plot  of  both  aquipotantials  chat  paaa  chrotgh 
Cha  cornar  (r»R,  *»0)  of  cha  cylindar  dacaralnaa 
a  raglon  chat  can  ba  made  cha  aama  pocancial  by 
forming  lc  of  high  paraaabilley  iron.  The  aan- 
nar  in  which  cha  ouCar-flald  la  aubdlvldad  by  cha 
pola-placa  that  Couchaa  Cha  r-R,  z»0  adga  la 
ahown  in  Fig.  1.  Of  couraa,  Chla  cachniqua  ia 
noc  faaalbla  should  Cha  and  macch  aquipocandal  lia 
balow  chac  of  cha  radlal-oacch  pocancial. 


llnaar  fiald  profila  on  Cha  axle  of  a 
YIG  cylindar. 


SYNTHESIS  07  Hi(r)  IS  A  VERY  THIS  DISK 

In  ordar  Co  normally  magnaciza  a  vary  Chin  disk 
of  magnadzadon  M  ao  aa  Co  produca,  within  it,  a 
praacribad  fiald  Hr(r),  ic  is  nacaaaary  co  firae 
axpand  it  in  tha  form 


Hx(r<ll) 


n«0  °2n  r 


2n 


(IS) 


If  Cha  film  la  locatad  at  tha  plana  z«0  and  Cha 
oucar-pocancial  la  axpandad  aa 


*  "  n»0  C2n+1  nr*!1*’1' 


(16) 


cha  boundary  condltiona  can  ba  maechad  (naglaceing 
fringing  at  tha  rim)  providad  Cj-  ♦  M  and 


C2n+l" 


(2n+l) I 


32n 


nil 


(17) 


SYNTHESIZES  LI.VEARLY  DISPERSIVE  WAgKETOElASTIC 
DELAY  LINES 

Ua  hava  uaad  cha  aynehaala  procadura  ouclinad 
abova  to  daalgn  and  build  a  apaclal  claaa  of  sag¬ 
as  coals  sdc  dalay  linaa,  naaaly,  llnaarly  dispar- 
alva  dalay  linaa  with  wida  inatancanaou*  bandwidth. 
Savaral  aecampta  2-4  hava  baan  aada  in  Cha  paae 


co  raallza  such  davicas  but  chair  auccaaa  waa  vary 
lialead.  Tha  difficulty  llaa  in  eha  face  that  cha 
raquirad  fiald  profila  ia  noc  achiavabla  la  aiapla 
YIG  gaosacrlaa  imoaraad  in  a  uniform  bias  magne¬ 
tic  fiald.  In  eontraet,  wa  hava  pravloualy  rapor- 
cadS  auccaaa  in  aynchaalzlng  laboratory  dalay  lSiee 
and  wish  hara  to  raporc  our  furchar  pregrass. 

First,  wa  ravlaw  eha  undarlying  physical  prln- 
ciplas  of  aagnatoalaacic  dalay  linaa  and  amphasiza 
cha  araas  of  importance  for  improved  davica  opera¬ 
tion  which  require  further  understanding.  Next, 
wa  give  design  parameters.  Third,  and  laat,  va 
supply  details  of  davica  construction  and  report 
measured  characteristics. 

DEVICE  PHYSICS 

In  cases  where  cha  internal  nagantic  field 
varies  slowly  enough,  spatially,  an  expansion  of 
cha  rf_  fields  and  magnetization  in  tarma  of  plana 
waves  is  still  possible.  Howavar,  tha  propagation 
constant  k  associatad  with  tha  wave  pgekat  Chan 
varies  with  position.  At  each  internal  point  z 
of  tha  single  crystal  YIG  rod,  the  frequency  of 
tha  z-dlrectad  magnetic  spin  waves  la  given  by 

w/YU  ■  H(z)  +  A/k2  +  Dk2  (18) 

0 

where  H(z)  is  tha  internal  magnetic  field,  y  is 
tha  gyrooagnatic  ratio,  and  A,S  are  constancs  ap¬ 
propriate  co  eha  material.  The  abova  relation  la 
a  simplification  only.  In  reality  chare  is  a 
aultitida  of  modes  Inversely  proportional  to  k,  as 
shown  in  Fig.  2a.  Alao  ahown  in  cha  figure  is  cha 
u/k  ralation  of  Cha  shear  elastic  wave  which  ia 
"coupled  to  tha  spin  wave*.  This  coupling  glvaa 
rise  to  cha  splitting  at  tha  crossovar  point;  l.s. 
tha  point  whars  tha  unperturbed  frequency  values 
ara  equal.  Tha  longitudinal  elastic  wave  la  not 
important  to  tha  dalay  line  operation  and  is  there¬ 
fore  omitted  from  tha  figure.  In, the  region  of 
■mall  k,  tha  dominant  term  is  A/k4  and  Cha  group 
velocity  Vg«ou/3k  is  nagatlva.  Tha  wave  in  this 

region  is  called  a  backward  magnacostatlc  wave 
(BMW) .  In  tha  high  k  region  tha  group  velocity 
ia  positive  and  Cha  wave  thara  ia  tarmad  exchange 
spin  vave  (SW) , 

To  understand  tha  operation  of  tha  single  end¬ 
ed  dalay  llna,  we  follow  the  path  of  a  particu¬ 
lar  vava  packet  of  frequency  a_  inside  tha  magne¬ 
tic  medium  as  shown  in  Fig.  2b.  Tha  pech  of  eha 
wava  stares  vary  near  cha  face  of  cha  rod  (point 
A  in  Fig.  2b)  in  cha  form  of  a  vary  low  k,  back¬ 
ward  magnetic  wava.  This  wave  la  ganaratad  by  an 
electromagnetic  signal  applied  co  an  sntsnna  plac¬ 
ed  in  close  proximity  to  tha  face  of  tha  rod.  As 
tha  wava  packet  rapidly  penetrates  into  tha  rod 
it  encounters  a  monotonlcally  increasing  magnetic 
field.  Tha  initial  group  (or  anargy)  velocity  is 
vary  high  but  decreases  rapidly  until  k  raachas 
tha  value  ky  •  (A/D)1'4  where  tha  velocity  la 

zero.  This  point  termed  tha  turning  point  is 
tha  point  of  furthest  penetration  of  tha  wava  into 
cha  rod  (point  B  in  Fig.  2b).  A  schematic  repre¬ 
sentation  of  tha  group  velocity  vs.  time  Is  shown 
in  Fig.  2c.  The  time  spent  in  tha  BMW  raglon  is 
vary  short,  on  the  ordar  of  a  few  nanoseconds  and 
will  henceforth  ba  naglactad  in  tha  calculation  of 
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che  total  delay  time. 

This  omission  should  not  be  construed  as  an  in¬ 
dication  that  the  BMW  region  is  of  little  impor¬ 
tance  in  the  operation  of  the  delay  line.  On  the 
contrary,  this  region  where  the  wave  packet  initi¬ 
ally  forms  is  of  crucial  importance  to  efficient 
delay  line  operation.  Unfortunately,  the  com¬ 
plicated  processes  associated  with  the  formation 
and  propagation  of  the  wave  in  this  region  are 
understood  only  qualitatively  and  in  scant  detalL 
If  che  magnetic  field  gradient  is  comparatively 
small6  at  the  turning  point,  k  will  suffer  no 
discontinuity  and  its  magnitude  will  continue  to 
increase.  Except  for  a  small  reflection,  the 
bulk  of  the  energy  continues  along  a  path  of  in¬ 
creasing  (k)  inco  the  exchange  spin  wave  region. 

The  wave  trajectory  is  now  coward  the  front 
face  of  the  rod  in  che  direction  of  decreasing 
magnetic  field.  In  this  SW  region,  lc  is  neces¬ 
sary  to  consider  both  the  radial  and  axial  compo¬ 
nents  of  the  magnetic  field.  Certain  radial  dis¬ 
tributions  tend  to  bend  the  wave  away  from  Che 
axis,  that  is  defocus  che  coherenc  energy  beam 
while  other  profiles  tend  to  focus  it  toward  che 
axis.  It  is  very  important  to  ensure  the  exls- 
tance  of  focusing  conditions  since  very  high  propa¬ 
gation  losses  occur  otherwise.  Too  much  focussing, 
on  the  other  hand,  is  counterproductive  since  it 
saturates  the  propagation  channel  by  creating  loc¬ 
al  regions  of  extremely  small  cross-section  and  this 
of  high  energy  density.  This  tendency  to  premature 
saturation  limits  the  power  handling  capabilities 
of  the  delay  line  and  hence  its  useful  dynamic 
range. 

The  determination  of  whether  a  focussing  or  a 
defocusslng  condition  exists  may  be  made  by  evalua¬ 
ting  a  dimentionless  quantity  Q  that  is  a  function 
of  the  H-field  and  its  first  and  second  spatial  der¬ 
ivatives,  evaluated  at  the  axial  point  z  under  con¬ 
sideration6.  A  concave  field  profile  (H"<0)  auto¬ 
matically  ensures  focussing  but  a  slightly  convex 
profile  is  also  allowed. 

Provided  a  focussing  condition  exists,  the  beam 
propagates  to  the  left  until  k  reaches  the  value 
k^-u/v,  where  v  is  the  velocity  of  the  elastic 

shear  wave.  This  is  defined  as  the  cross-over 
point  where  the  unperturbed  frequencies  of  the  elas¬ 
tic  and  spin  waves  are  equal  (point  C  in  Fig .2b)  . 
The  distance  between  the  turning  point  and  the 
cross  over  point  traversed  in  this  region  is  very 
small,  but  the  wave  on  the  other  hand  is  a  slow  one 
with  its  maximum  velocity  reached  at  the  cross  over 
point.  As  a  result,  che  time  spent  in  the  SW  re¬ 
gion  accounts  for  a  major  portion  of  the  total  de¬ 
lay  time  (more  than  half  in  the  case  of  a  linearly 
varying  field  profile) .  Representative  values  of 
the  distance  and  velocity  far  TIG  at  4  GHz  and  a 
magnetic  field  gradient  of  106-  Oe/cm  are  20mm  for 
che  distance  and  1.15xl04  cm/ sec  for  the  maximum 
velocity.  This  velocity  is  only  3Z  of  the  shear 
elastic  wave  velocity. 

The  coupling  between  the  spin  wave  and  the 
elastic  wave  causes  a  substantial  conversion  of  the 
spin  wave  packet  inco  an  elastic  shear  wave  at  the 
cross-over  point.  Under  the  same  stipulations  as 
before,  namely  that  the  field  gradient  is  not  too 
steep,  very  high  conversion  efficiencies  are  pos¬ 
sible.  The  elastic  wave  continues  to  propagate 


unaffected  by  the  magnetic  field  towards  the  face 
of  the  rod  where  it  undergoes  a  reflection  In  the 
highly  polished  surface  (point  D  in  Fig.  2b).  From 
this  point  the  wave  packet  retraces  Its  former 
steps  through  the  SW  and  3MU  regions  and  is  picked 
up  at  the  surface  in  the  form  of  an  electromagne¬ 
tic  wave.  This  may  be  done  by  the  same  antenna 
used  to  launch  the  input  signal,  or  by  another 
antenna  placed  in  close  proximity  to  it  when  two 
port  operation  is  desired.  The  delayed  signal 
whose  rather  complex  trajectory  we  have  just  fin¬ 
ished  describing,  is  the  main  output  of  the  delay 
line,  the  so  called  1st  echo.  Mot  all  the  wave 
energy,  however,  is  extracted  by  the  receiving  an¬ 
tenna  and  the  unextracted  portion  is  launched  back 
into  the  crystal  to  arrive  back  at  the  rod  face 
as  a  2nd  echo.  This  spurious  3lgnal  is  delayed  by 
twice  the  delay  of  the  1st  echo  and  is  of  reduced 
amplitude.  Higher-order  echoes  exist  as  well  at 
further  reduced  amplitudes. 

When  two  port  operation  of  the  single  ended 
delay  line  is  required,  the  direct  electromagne¬ 
tic  feed  the  rough  between  the  input  and  output 
antennae  which  are  at  close  proximity  causes  an¬ 
other  spurious  output  signal.  This  sometimes 
bothersome  spurious  whose  high  amplitude  had  been 
reported,^  could  be  reduced  substantially  by  im¬ 
proving  the  wave  coupling  antennae . 

The  total  transit  time  of  the  1st  echo  is  the 
Sum  of  the  times  spent  as  an  exchange  spin  wave 
and  as  an  elastic  shear  wave.  This  is  given  by 


T 


(19) 


where  z T>  z^  are  the  positions  of  the  turning  poirt 

and  the  cross-over  points  respectively,  v  is  the 
velocity  of  the  elastic  shear  wave,  and  v^«3oj/3k 

is  the  group  velocity.  In  the  above  expression 
for  evaluating  T  we  neglect  the  transit  times 
through  the  BMW  and  the  cross  over  regions. 

We  therefore  assume  u)/y«H(z)+Dcxk<-,  and  by 

expanding  H(z)  in  a  Taylor  series  around  z^,  ob¬ 
tain  the  approximate  expression  for  T 


T 


2w 

U  YH'vb 
o 


tan  LS 


(20) 


where  3-  ;jJ2D  H"  /vH'  is  a  dimensionless  para¬ 
meter,  H',  H"  are  the  first  and  second  spatial  dsri- 
vatives  of  the  axial  magnetic  field  evaluated  at 
the  cross  over  point  z^,  and  is  the  exchange 

constant.  For  TIG,  when  converted  to  CCS  units, 

D  ■  5xl0-’  Oe  cm  ,  v  *  3.84x10^  cm/oc  and 
ex 

Y“2.8MHz/0e:  In  evaluating  T  in  the  above  equa¬ 
tion,  we  assumed  kj,=0  and  k  "Wv. 

Equation  (20)  for  T,  '  can  be  used  to  com¬ 
pute  the  delay  time  for  a  given  field  profile  or 
alternatively  to  synthesize  the  necessary  profile 
for  achieving  a  desired  delay  characteristic.  In 
the  latter  case,  T(uO  is  a  given  function  of  fre¬ 
quency.  When  a  field  synthesis  is  attempted,  a 
word  of  caution  is  appropriate.  For  a  specified 


T(ui),  Che  synthesized  field  profile  nay  turn  ouc  to 
be  physically  unrealizable  or  If  physically  possible, 
mighc  turn  out  to  be  a  defocusing  profile,  render¬ 
ing  It  Impractical.  For  example,  the  profile  pro¬ 
posed  for  a  constant  delay8  is  a  highly  defocus¬ 
lng  one. 

In  the  limits  of  Dg^-8-0  we  obtain  from  Eq.  (20) 
the  position  z  as 


U  YH 
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T(ia)du  +  C 


(21) 


where  C  is  an  Integration  conscant  whose  value  may 
be  chosen  at  will. 

The  field  profile  as  a  function  of  position  z 
along  the  axis  of  the  rod  can  be  found  by  Invert¬ 
ing  the  above  expression  for  z^.  When  higher  ac¬ 
curacy  of  the  prescribed  delay  characteristics  owr 
wide  instantaneous  bandwidth  is  required,  a  better 
determination  of  the  field  profile  may  be  warranted. 

An  Important  delay  characteristics  Is  a  linear 
dispersion  where  che  dealy  time  Increases  linearly 
with  frequency.  In  this  case  T«To+Dui  .  By  solving 

Eq.  (21),  and  when  00  is  chosen,  the  required  mdal 
field  profile  Is  a  linear  one,  given  by 
2T 

H  -  -r^  + 


YD 


4z 

YvD 


(22) 


The  linear  dispersion  coefficient  D,  is  therefore 
given  by  D-4/V  YvH'  whre  H'  is  the  field  gradient. 
When  Ci*0  Is  chosen,  the  resulting  profile  is  non¬ 
linear. 

DESIGN  PARAMETERS 

The  most  Important  operational  parameters  in 
this  class  of  devices  together  with  our  design  goals 
are  listed  below: 

Insertion  loss  and  flatness 

For  the  device  to  be  of  practical  use,  it  should 
have  no  more  than  a  moderate  amount  of  Insertion 
loss  and  reasonable  variation  of  such  losses  across 
the  frequency  band.  Obviously  any  matching  arrange¬ 
ment  to  minimize  the  loss  has  to  be  broad  band 
since  no  external  tuning  at  each  frequency  is  al¬ 
lowed  . 

Linear  dispersion  value  and  accuracy 

Potential  uses  of  and  interest  in  devices  poss- 
sslng  a  wide  range  of  dispersion  factors  exist. 
Based  on  our  laboratory  experience,  we  feel  chat 
values  from  D-. INS/MHz  to  2  NS/MHz  are  currently 
achievable,  not  necessarily  over  the  same  band¬ 
width.  For  certain  applications  a  useful  criter¬ 
ion  is  the  time  bandwidth  product  defined  as  DJlu) 
and  a  practically  achievable  number  of  >500.  The 
deviation  from  linearity  is  also  dependent  on  che 
bandwidth  and  values  in  the  1-3Z  range  over  band¬ 
width  of  up  to  1500  MHz  are  achievable.  There  is 
a  wide  room  for  tradeoffs  among  the  three  para¬ 
meters  dispersion  value,  allowable  percentage  de¬ 
viation  from  linearity  and  frequency  bandwidth. 
Frequency  range  and  lnstantaneoua  bandwidth 

Aside  from  bandwidth  limitations  which  can  be 
traded  off  as  discussed  above,  there  is  an  upper 
limit  imposed  by  che  appearance  of  che  spurious  2nd 
echo.  As  was  earlier  explained,  the  1st  echo  gen¬ 
erates  its  own  echo  and  che  2nd  echo  free  bandwidth 


is  therefore  defined  as  the  Lw  over  which  the  only 
delayed  output  present  is  the  1st  echo.  This  band¬ 
width  increases  upon  increasing  the  frequency  (an 
octave  being  the  theoretical  maximum)  and  a  1GHz 
or  more  at  s  band  was  sec  as  a  goal.  Device  opera¬ 
tion  in  the  frequency  range  of  .5-5CH z  has  been 
verified  with  }GBz  being  a  laboratory  instrumenta¬ 
tion  limitation.  We  belive  that  useful  operation 
at  C  and  possibly  x  bands  in  feasible  with  instan¬ 
taneous  bandwidth  exceeding  2CHz. 

Input/Output  isolation 

Since  che  nondisperslve  electromagnetic  leak¬ 
age  appears  instantly  at  the  output,  a  high  ampli¬ 
tude  may  be  tolerable  for  certain  applications. 

For  seme  important  applications,  however,  this  is 
not  so.  Although  it  may  seem  that  due  to  the 
close  proximity  of  the  input  and  output  antennae, 
high  isolation  could  not  be  achieved,  our  design 
goal  was  to  limit  the  leakage  to  a  level  at  which 
its  amplitude  is  no  higher  than  che  amplitude  of 
the  1st  echo  signal. 

Dynamic  range 

Our  delay  line  is  a  magnetic  device  and  like 
all  such  devices,  its  behavior  changes  drastically 
when  high  power  rf  signals  are  present.  These 
high  power  effects  limit  the  dynamic  range  for  use¬ 
ful  operation.  Since  high  energy  densities  are  the 
culprits  an  obvious  way  to  increase  the  dynamic 
range  is  to  increase  the  effective  cross  sectional 
area  of  the  signal.  This  may  be  accomplished  by 
optimizing  the  coupling  to  the  BMW  and/or  by 
decreasing  the  amount  of  focussing  of  the  signal 
while  it  traverses  the  exchange  spin  wave  region. 

As  large  a  dynamic  range  as  possible  is  obviously 
desirable  for  many  applications.  For  our  ealier 
work  we  achieved  ranges  of  Z5-30db,  and  any  im¬ 
provement  was  deemed  welcome. 

DEVICE  CONSTRUCTION  AND  MEASURED  CHARACTERISTICS 

An  exploded  view  of  a  packaged  device  is 
shown  in  Fig.  3D  Notice  the  position  of  the  YIG 
rod  between  the  synthesized  soft  iron  pole  pieces. 
The  pole  pieces  were  designed  to  implement  the 
radial  match  only.  Other  surfaces  were  omitted 
to  allow  access  for  the  wave  coupling  structure. 

The  resultant  internal  field  is  therefore  only  an 
approximation  to  the  exact  requirement.  Its  axial 
position  is  adjustable,  and  it  is  held  securely  by 
the  tension  supplied  by  the  adjustment  screw  in 
the  front,  and  the  locking  cylinder,  which  presses 
on  the  back  pole  pike,  at  the  back.  The  input 
and  output  antennae  which  are  loops  of  5  mil 
enameled  wire  soldered  to  semirigid  cables,  are 
pressed  to  the  front  face  of  the  crystal  by  a  snail 
cylindrical  teflon  piece.  The  diameter  of  the 
loops  is  30  mils  and  they  are  placed  side  by  side 
off  the  center  of  the  face  of  the  YIG  rod.  The 
parallel  feed  lines  joining  each  loop  to  the  semi- 
rized  cable  intersect  at  approximately  90*  on  the 
face  of  the  red.  The  magnetic  flux  is  supplied 
by  Alnico  8  permanent  magnet  rings.  The  dimens  lavs 
of  the  package  are  2"  dia.  by  2  1/2"  long  and  it 
accomodaces  a  120  ails  dia.  by  200  mil  long 
(IDO)  YIG  rod.  Smaller  size  packaged  devices  have 
been  built  and  successfully  tested. 

With  the  improved  coupling  structure  described 
above  we  observed  simultaneous  improvement  in  sev¬ 
eral  device  characteristics.  Specifically,  overall 
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lnsartion  loss  vaa  reduced  by  up  to  lOdb;  loaa  rip- 
pla  vaa  limited  co  1-3  db ;  dynamic  ranga  maaaurad 
between  the  noiaa  floor  and  Cha  onaac  of  nonlinaar 
high  povar  affacca  incraaaad  by  at  laaat  lOdbj 
alactromagnaele  laakaga  faadthrough  could  ba  ra¬ 
ducad  to  a  laval  at  which  it  vaa  equal  to  or  lovar 
than  tha  laval  of  tha  frequency  diaparaiva  time- 
dalayad  primary  aignal.  Input/output  laolation 
incraaaad  aa  a  raault  by  up  to  20db.  Tha  apurioua 
2nd  ache  aignal  could  alao  ba  auppraaaad  to  lavala 
lovar  than  tha  aignala'  by  up  to  20db. 

In  Figure  4  va  aaa  oacllloacopa  tracaa  which 
highlight  tha  charactarlatlca  of  a  rapraaantatlva 
packaged  davlca  operating  at  S  band.  Tha  davlca 
contained  aynthaalaad  pola  placaa  daalgnad  for  a 
field  on  axis  of  3004300c  (c  in  am).  In  all  four 
parta  of  tha  figure  ahort  input  pulses,  40  NSec 
wide)  warn  applied. 

Part  (a)  ahowa  tha  output  obtained  at  dlacrate 
input  fraquanclaa  aa  marked.  Tha  let  and  2nd  achoaa 
era  aean  to  ba  delayed  in  proportion  to  tha  frequen¬ 
cy.  the  latter  at  twice  tha  rata  of  the  formate. 

Tha  laakaga  faadthrough  la  not  tlma-dalayad  and 
lta  poaltlon  doaa  not  vary  with  frequency.  We  aaa 
that  tha  amplltudea  of  tha  apurioua  laakaga  and  2nd 
echo  aignala  are  lower  than  that  of  tha  primary  lat 
echo  aignala.  tha  latter  by  aa  much  aa  20  db.  By 
awaaplng  tha  ahort  input  pulaaa  acroaa  tha  frequency 
band  and  applying  tha  awaaplng  voltaga  ramp  to  tha 
vertical  axia  of  tha  oacllloacopa,  the  frequency  va. 
dalay  time  dependence  la  obtained.  Thla  ia  ahovn  in 
part  (b).  Tha  dlaperalon  la  .3  NSec/MHa  and  tha 
frequency  band  lOOOMHa.  Tha  linearity  of  tha  dla- 
paralon  la  quite  good  with  a  2X  deviation  from  a 
atralght  line  obeervad.  By  removing  tha  awaaplng 
ramp  from  tha  vertical  axia  of  cha  oacllloacopa 
while  continuing  to  awaap  tha  input  aignal,  wa  ob¬ 
tain  part  (c)  which  ahowa  tha  amplitude  of  the  out¬ 
put  aignal  aa  a  function  of  lta  frequency.  The  in- 
aartlon  loaa  of  tha  primary  lat  echo  la  27-30  db 
acroaa  tha  1000  MHx  frequency  band  and  a  negligible 
amount  of  aignal  dlatortlon  la  obaarvad.  Tha  apur¬ 
ioua  aignala  era  at  raducad  amplltudea  with  input/ 
output  laolation  of  33-40  db,  and  2nd  acho  lavala 
of  13-29  db  below  tha  lavala  of  tha  primary  lat 
achoaa. 

Bart  (d)  ahowa  operation  at  tha  fixed  fraquaney 
of  2.3  CHt  where  cha  input  povar  lncraaaaa  from 

tha  upper  trace  downward a . _ Tha  relative  input 

power  lavala  are  marked  in  tha  figure.  Tha  lovaae 
trace  claarly  ahowa  the  nonlinear  high  povar  effst 
on  Cha  aignal;  tha  output  pulaa  ahowa  signs  of  a 
breakdown  and  la  no  longer  a  delayed  replica  of 
tha  input.  4  complete  breakdown  occura  upon  a  fur¬ 
ther  lncreaaa  In  tha  input  power  (not  ahovn  in  the 
figure) .  The  threahold  level  at  which  nonlinaar 
raaponaa  aaca  in  la  app.  -lObm  and  tha  dynamic 
range  for  linear  operation  la  44db.  Thla  waa 
maaaurad  in  a  aatup  which  included  a  47  db  S  band 
amplifier  with  a  3  db  noiaa  figure.  It  can  ba 
aean  In  part  (d)  that  tha  2nd  acho  aaturataa  at 
lower  power  lavala  than  cha  lat  echo.  Thla  auggaata 
that  the  energy  denalty  of  the  aignal  tanda  to  ln¬ 
creaaa  in  proportion  to  ita  dalay  time  and  hence 
that  amallar  dynamic  ranges  era  to  ba  axpectod  in 
dalay  linea  exhibiting  long  dalaya.  Thla  tendency 
may  ba  offaac  in  low  focuaalng  internal  magnetic 
field  profllea  where  tha  aignal  energy  denaltlaa 
are  lower. 


Tha  magnetic  field  on  tha  axia  of  tha  TIG  rod, 
calculated  by  fitting  tha  data  repraaented  In  Fig. 

4  to  Bq.  (20)  for  tha  dalay  time  T,  ia  B;(r>0)  » 
49041240a  -  434a2  0a  (a  in  nan) .  Thla  field  la 
qulta  different  from  tha  aaaumed  field  at  tha  out- 
aat  of  the  ayntheela.  Two  pointa  ehould  ba  borna 
in  mind  in  connection  vlth  thia  diacrapeney.  Tha 
flrat  one  la  that,  aa  mentioned  above,  only  tha  ra¬ 
dial  match  of  tha  aynthaala  waa  implemented.  Thia 
will  diatort  tha  field  in  tha  active  region  of  tha 
dalay  line  which  axtenda  to  a  depth  of  only  l-2mm 
into  tha  TIG  rod.  Tha  field  on  cha  axia  of  a 
naked  TIG  rod;  whan  placad  in  a  uniform  DC  blaa 
field,  haa  a  vary  email  linear  and  a  larga  quadra¬ 
tic  apatlal  coafflcianta.  Tha  required  axial  fidd 
on  cha  other  hand,  ia  draatlcally  different,  hav¬ 
ing  no  quadratic  and  a  larga  linear  eoefflclenta; 
thla  difference  haa  to  ba  overcome  by  cha  eynthe- 
aiaad  pola  piacaa. 

Tha  aacond  point  ia,  that  a  act  of  pola 
plaeaa  daalgnad  for  a  apeclflc  field  profile  la 
capable  of  aypportlng  a  whole  family  of  profllea 
by  adjuating  the  flux  paaalng  through  them.  Thla 
ia  analogoua  to  tha  eaaa  of  a  capacitor  where  the 
electric  field  may  ba  varied  by  adjuating  tha 
charge  accumulated  on  lta  plataa.  Tha  moat  advan- 
tagaoua  field  profile  for  efficient  device  opera¬ 
tion  aa  determined  by  rf  conalderatlona,  may  not 
ba  tha  original  profile  from  which  the  pola 
plaeaa  ware  originally  aynthaslzed. 
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Fig.  2  Diagrams  explaining  magnetoelastic  delay 
line  operation;  (a)  w-k  relation, 

(b)  axial  cross-section  of  the  YIG  rod 

(c)  group  velocity  vs.  time. 
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Figure  4  Measured  delay  line  performance  (refer  to  text  for  discussion). 
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The  following  sets  of  figures  provide  further  information  concerning  the 
characteristics  of  the  Time  Prism  Filters  that  were  developed  by  employing 
the  field  synthesis  techniques. 


FREQUENCY  (  GHz 


Figure  1 .  Characteristics  of  packaged  device  after  its  removal  from  the 
magnetizing  field,  (a)  Frequency  as  a  function  of  delay  time  with  a 
straight  line  superimposed  for  comparison .The  dispersion  is  . 35NSec/MHz 
(b)  Signal  amplitude  as  a  function  of  frequency . Insertion  loss  39-^5  db. 
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Figure  2.  Time  Prism  characteristics  after  coupling  optimization. 

(a)  Device  output  at  four  distinct  frequencies  as  marked . Leakare  amplitude 
equal  to  or  lower  than  first  echo,  (b)  Frequency  as  a  function  of  delay 
time . Dispersion  is  .3N3ec/MHz.  (c)  Si/paal  amplitude  as  a  function  of 
frequency . Insertion  loss  32~35db, leakage  amplitude  lower  than  or  equal  to 
first  echo,  second  echo  amplitude  lower  by  at  least.  !5db. 
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Figure  3.  Dynamic  range .Relative  input  power  is  as  marked,  (a)  Fixed 
frequency  2.8  GHz. Lowest  trace  shows  nonlinear  high  power  effects, linear 
dynamic  range  44db.  (b)  Full  band  sweep  of  1000MHZ  in  each  trace, linear 
dynamic  range  44db 
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Figure  4.  Characteristics  of  a  filter  in  a  laboratory  setup,  (a)  Frequency 
as  a  function  of  dela^  time.  Dispersion  is  .38  NJec/KHz.  (b)  Simal  ampli¬ 
tude  as  a  function  of  frequency. Very  flat  insertion  loss  of  35+. 5db 


Figure  5*  Filter  charge  teris ',i  in  -■  laboratory  setup,  (a)  Frequency  as  a 
function  ol  delay  time,  low  c  i  sp-'r  ;io,  ■  .  1*7  A'Sec/lv.hz .  (fc)  Signal  amplitude  as 
a  function  of  frequency  under  some  conditions  as  (a), low  insertion  loss 
25-28  db  across  1200  GHz.  (c)  different  settings  from  those  in  ( a  )&(  b ) ,  same 
pole  pieces  and  YIG  rod .  Disne  'tier;  .-,3'  .  IdNSec,/  Hz ,  deviation  from  linearity 
of  less  than  1.51  across  a  12  .  -Mir  oandwidtn, very  fiat  insertion  loss  33*.5db 
T(  2 . 75GHz  )  -  680  NGec. 
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FIGURE  2 


OUTPUT  OF  TIME  PRISM  MICROWAVE  PULSE  ANALYZER,  (a)  FREQUENCY/DELAY 
TIME  CHARACTERISTICS,  SWEPT  INPUT,  (b)  RESPONSE  TO  THREE  DISTINCT  IN¬ 
PUT  FREQUENCIES,  VERTICAL  POSITION  PROPORTIONAL  TO  FREQUENCY,  (c) 
AMPLITUDE  AS  A  FUNCTION  OF  FREQUENCY,  SWEPT  INPUT,  (d)  SUPERPOSITIPN 
OF  SYSTEM  RESPONSE  TO  THREE  DISTINCT  INPUT  FREQUENCIES.  ALL  INPUTS 
ARE  CW  AND  APPLIED  TO  INPUT  1,  GATING  IS  INTERNAL. 


RESPONSE  OF  TIME  T WAV  P..I  SF.  ANALYZER  TO  TWO 
SIMULTANEOUSLY  AT.  ;  : i  ;  :  ..  '  Art  V'VAi.S.  fa)  f1  =  2.9  GHz, 

Af  =  f2  -  f1  --  Si”.  ‘-V;  V  V  -  M“'z  (c)  Af  =  2.5  MHz 
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Thin  Film  Input/Output  Coupling  Structures 

The  antennae  are  the  major  components  in  determining  both  isolation  from 
input  to  output  and  the  coupling  efficiency.  Previously  optimization  of  these 
two  specifications  was  difficult  at  best  due  to  a  number  of  problems  inherent 
in  using  hand  formed  loops,  including:  (1)  inability  to  make  reproducible 
loops,  (2)  questionable  reliability  of  wire  loops  due  to  embrittlement  from 
work  hardening,  (3)  limited  isolation  due  to  overlapping  of  input  and  output 
loops.  The  proposed  design  attempts  to  overcome  these  difficulties.  The  thin 
film  antennae  will:  (1)  be  reproducible,  given  a  design  mask,  (2)  not  be 
subject  to  embrittlement  through  work  hardening,  because  the  antennae  will  he 
sputtered  and  etched  on  a  rigid  substrate  of  alumina,  (3)  improve  isolation 
with  nonoverlapping  loops  and  perhaps  ground  separators.  In  addition  to 
accommodating  for  proposed  antennae,  the  proposed  delay  line  fixture  should 
allow  for  interchanging  of  the  iron  pole  pieces,  which  shape  the  magnetic  field 
applied  to  the  YIG  rod.  Ultimately  this  ability  to  interchange  pole  pieces  and 
antennae  will  open  the  way  for  a  future  series  of  experiments  to  maximize  power 
handling  capabilities,  minimize  attenuation  and  maximize  antennae  isolation  and 
coupl ing. 

The  SM  thesis  of  Leslie  Itano,  now  nearing  completion,  has  as  its  main  goal 
the  design  of  a  magnetoelastic  delay  line  utilizing  thin  film  input  and  output 
antennae.  Previous  delay  lines  utilized  resonant  cavities  or  hand  formed  wire 
loops  to  provide  the  electromagnetic  input  pulse  and  recieve  the  output  pulse. 
Although  the  delay  line  performance  with  the  hand  formed  wire  loops  was  excellent, 
the  results  were  extremely  difficult  to  reproduce,  requiring  many  hours  of  pain¬ 
staking  initial  assembly  work  and  many  more  hours  of  tuning  at  the  test  bench  for 
each  antennae  pair.  Then,  with  all  this  loop  adjustment,  the  wires  sometimes 
became  work  hardened  and  embrittled,  resulting  in  highly  tuned  out  unreliable 
antennae.  In  contrast,  the  thin  film  antennae  are  made  photol i thographi cal  1 y 
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using  a  mask  with  a  specific  design.  This  method  produces  consistent,  easily 
reproducible  and  reliable  results.  The  fabrication  time  for  thin  film  circuits 
is  also  significantly  less  due  to  several  factors.  First  of  all,  a  pair  of 
antennae  is  produced  on  a  single  substrate  so  that  several  substrates  can  be 
processed  simultaneously.  Secondly,  once  a  substrate  has  been  processed  and  the 
antennae  pair  is  etched  onto  the  substrate,  the  circuit  is  fixed.  The  thin 
film  antennae  are  not  malleable  like  the  wire  loop  antennae,  eliminating  the 
tuning  step  and  further  reducing  fabrication  time.  Thus,  the  thin  film  antennae 
is  not  only  more  reproducible  and  reliable,  but  it  requires  less  fabrication  time 
than  the  hand  formed  loop  antennae. 

A  secondary  achievement  of  this  thesis  is  the  design  of  a  flexible  del  ay 
line  fixture.  By  "flexible",  we  mean  that  the  fixture  is  designed  to  allow  for 
a  range  of  antennae  of  substrate  thickness,  pole  pieces  which  shape  the  internal 
rod  profile,  and  thin  film  antennae  designs.  This  ability  to  interchange 
components,  and  to  interchange  them  quickly  makes  this  design  well  suited  to 
research . 

This  new  fixture  design  also  retains  the  advantages  of  previous  designs.  The 
new  fixture  is  compact;  it  can  be  self-contained  since  allowance  is  made  for  the 
use  of  ring  magnets.  The  pole  pieces  have  been  synthesized  to  provide  linear 
dispersion . 


Description  of  Fabrication  Sequence  (Summary) 


This  section,  which  will  form  a  portion  of  itano's  thesis, 

describes  the  fabrication  of  the  YIG  delay  line.  First,  an  over¬ 
view  of  the  fabrication,  and  associated  problems  will  be  pre¬ 
sented.  The  three  major  categories  which  will  be  included 
are:  (1)  producing  the  thin  film  antennae,  (2)  manufacturing 

the  fixture,  and  (3)  assemblying  the  delay  line.  Second, 
a  summary  of  the  fabrication  sequence  will  be  provided  as  an 
aid  for  subsequent  work. 

4. A.  Producing  Thin  Film  Antennae 

The  first  step  in  producing  a  thin  film  antennae  is  mask 
fabrication.  This  step  involves  designing  the  antennae, 
reducing  the  pattern  and  transferring  it  to  a  glass  photo¬ 
graphic  plate.  This  glass  plate,  complete  with  antennae 
pattern,  is  the  mask.  In  the  meanwhile,  an  alumina  substrate 

has  been  sputtered  by  an  outside  recorder  with  thin  film 

o 

chrome-gold.  The  thin  layer  (10C-200A  )  of  chromium  acts  as 
an  adhesion  layer  between  the  substrate  and  the  thicker  gold- 
layer  (200  micro  inches) .  The  mask  is  now  used  as  a  photo¬ 
lithographic  negative  to  selectively  etch  both  chromium  and 
gold  layers;  the  pattern  remaining  on  the  substrate  is  the 
thin  film  antennae. 

Developing  Etching  Process 

The  complete  "etching  process"  for  thin  film  chromium 
gold  circuits  has  been  well  established.  In  this  overview. 


we  will  first  review  the  general  process,  explaining  the 
necessity  of  each  step.  Next,  we  describe  the  special 
"challenges"  presented  by  the  antennae  designed  for  this  thesis, 
and  the  attempts  to  overcome  these  challenges.  Lastly,  we 
summarize  the  final  approach  used  in  this  thesis. 

The  "etching  process"  begins  by  spinning  on  a  thin 
layer  of  photosensitive,  etch  resistant  liquid  to  the  sub¬ 
strate.  The  spin  rate  and  viscosity  of  this  "photoresist" 
liquid  determine  the  thickness  of  this  layer.  The  substrate 
is  then  "prebaked"  to  harden  the  photoresist,  making  it  more 
resistant  to  breakdown,  i.e.,  local  voids  or  weak  points 
in  the  photoresist  caused  by  handling  in  the  exposure  step 
to  follow.  "Breakdown"  in  the  photoresist  will  allow  the 
etchant  to  seep  in,  etching  the  gold  and/or  chromium. 

The  glass  mask  already  prepared  is  a  negative;  that  is, 
the  design  area  is  blackened,  blocking  out  light,  while 
the  remaining  area  is  totally  transparent.  This  mask  is 
aligned  to  the  substrate  edges,  contact  is  made  with  the 
substrate,  and  the  photoresist  is  exposed  to  ultra  violet 
light  in  those  areas  where  the  mask  is  transparent.  The 
photoresist  is  chemically  developed,  toughening  all  of  the 
unexpected  areas.  The  exposed  areas  are  washed  away,  leaving 
behind  the  design's  protective  photoresist  pattern.  The 
substrate  is  postbaked'  to  toughen  the  remaining  photoresist. 

Since  the  design  area  is  completely  protected  by  photo¬ 
resist,  the  remaining  unprotected  area  can  be  etched  away. 


Figure  4.1a  -  Thin  film  antennae  designed,  set  in  delay  line 
housing 


Figure  4.1b  -  Kn  la  r«  foment  of  couplin')  region 
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The  top  gold  layer  is  etched  first  with  a  potassium  iodide 
solution,  while  the  lower  chromium  layer  is  etched  with  a 
KMnO^  and  sodium  hydroxide  solution.  Finally,  the  protective 
photoresist  layer  is  stripped  away  with  an  industrial  strip¬ 
per,  leaving  behind  the  intact  design  pattern. 

The  antennae  designed  for  this  thesis  (see  Figure  4.1) 
requires  several  sputtered  through  via  holes  to  ground  in 
close  proximity  to  narrow  (3-4  mils  wide)  lines.  The  chal¬ 
lenge  was  to  protect  the  via  holes  without  bridging  over 
to  the  narrow  lines.  Using  the  established  processes  and 
taking  no  special  precautions  to  protect  the  via  holes,  gold 
etch  seeped  into  the  holes,  partially  etching  the  insides 
of  the  holes.  Also  large  voids  were  present  in  those  regions 
where  the  photoresist  layer  was  particularly  thin,  such  as  the 
edge  of  the  via  holes.  The  first  attempt  to  solve  this  problem 
was  a  fourfold  increases  in  the  photoresist  thickness  (to  4 
micron) .  As  hoped,  this  change  resulted  in  improvement  co¬ 
verage  for  the  via  holes;  however,  voids  still  appeared 
at  the  edges.  In  addition,  the  etch  factor  has  increased 
significantly,  causing  a  dramatic  change  in  line  width  or 
effective  line  impedence.  The  etch  factor,  as  shown  in  Figure 
4.2,  is  the  difference  between  the  top  surface  or  apparent 
edge  of  the  photoresist  and  the  edge  of  the  etched  line.  In 
general,  the  etch  factor  can  be  due  to  breakdown  or  softening 
of  the  photoresist  edge  which  allows  the  etchant  to  seep 
underneath.  In  the  case  of  the  thicker  photoresist,  the 
increase  in  etch  factor  is  likely  due  to  diffraction  of  the 
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exposure  light.  Figure  4.3  illustrates  how  this  diffraction 
causes  extra  exposure  of  the  thicker  photoresist  near  the  con¬ 
ductor  surface,  exactly  where  its  protection  is  needed. 

Figure  4.4  shows  how  a  vertical  visual  inspection  completed 
before  etching  would  not  detect  the  diffraction  caused  dif¬ 
ference  between  conductor  surface  or  effective  photoresist 
edge  and  top  surface  or  apparent  photoresist  edge.  So  this 
modification  brought  some  improvement,  but  proved  to  be  more 
problematic  than  helpful. 

Returning  to  the  original  1  micron  thick  photoresist  layer 
the  etch  factor  was  further  improved  by  increasing  the  post¬ 
bake  temperature  from  90°C  to  150°C.  This  seemed  to  further 
harden  the  photoresist  edges,  resulting  in  less  etchant  seep¬ 
ing  underneath  and  etching  the  metallic  layer. 

The  remaining  problem  of  protecting  the  sputtered  through 
holes  was  solved  by  manually  applying  a  thick,  acetone  based 
resist  to  the  holes  with  a  few  strands  from  a  camel  hair 
brush.  Since  the  holes  were  only  3  mils  in  diameter  with  ad¬ 
jacent  conductor  lines  only  2-3  mils  away,  this  task  reauried 
steady  hands,  but  resulted  in  totally  protected  sputtered 
through  holes . 

4.B.  Manufacturing  the  Fixture 

Manufacturing  the  delay  line  fixture  involves  designing 
the  components  of  the  housing  and  then  machining  them.  Since 
rusting  of  the  pole  pieces  had  altered  the  field  profile  of 
previous  efforts,  these  parts  were  flashed  with  Ni  for  pro¬ 


tection. 
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Figure  4.4  -  Shows  how  vertical  inspection  does  not  detect 
diffraction  caused  etch  factor 
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4 . C .  Assemblvinq  the  Delay  Line 

The  final  delay  line  assembly  starts  by  putting  together 
the  connectors  and  main  housing  (see  Figure  4.1).  The 
next  step  is  to  fix  all  of  the  flexible  parameters,  including 
raising  the  front  pole  piece  to  compensate  for  the  substrate 
thickness,  installing  the  specially  designed  pole  pieces,  and 
epox.ying  the  thin  film  antennae  to  the  front  pole  piece. 

The  electrical  connections  between  the  substrate  and  coaxial 
connectors  are  made  with  a  1  mil  x  25  mil  x  100  mil  gold 
ribbon  which  has  been  gap  welded  to  the  substrate,  then 
epoxied  to  the  connector.  Gap  welding  is  a  general  technique 
for  thermal  compression  bonding  gold  ribbon  to  gold  metalli¬ 
zation  (including  thin  film  gold  metallization) .  After  fin¬ 
ally  installing  the  remaining  components,  including  YIG  rod, 
front  pole  piece  and  retainers,  the  delay  line  is  complete 
and  ready  to  test. 

The  following  detailed  summary  is  provided  as  aid  for 
subsequent  work: 

Description  of  Antennae  Fabrication 

A.  Mask  Fabrication 

1.  Design  antennae 

2.  Transfer  design  (at  20x)  to  rubvlith 

3.  Photoreduce  patterns 

4.  Transfer  (xl)  patterns  to  glass  photographic 
plates 

B .  Sputtered  Substrate  Procurement 
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1.  Order  alumnina  to  size — 15  mil  thickness  ordered 
to  minimize  gap  between  rod  and  pole  piece. 

Width  and  length  determined  by  fixture. 

2.  Order  holes  drilled  for  sputtered  through  holes 

3.  Sputter  substrates  with  chrome  (for  adhesion)  and 
200  y  inches  of  Au. 

C .  Substrate  Etching  Summary 

1.  Spin  AZ  1350J  photoresist  on  to  metallized 
substrates . 

2.  Bake  substrates  at  90°  for  25  min.  to  dry  photo¬ 
resist. 

3.  Align,  glass  photographic  mask  to  both  substrate 
edges  and  sputtered  through  hole  sites. 

4.  Expose  photoresist  with  ultraviolet  light, 
using  glass  photographic  mask  as  a  negative. 

5 .  Apply  protection  paint  to  sputtered  through 
holes . 

6.  Develop  photoresist  with  Shipley  AZ  Developer, 
removing  exposed  portions  of  photoresist. 

7.  Postbake  substrates  at  150°C  for  25  min.  to 
toughen  photoresist. 

8.  Make  batches  of  both  chromium  and  gold  etch 
solutions  as  follows: 


Gold  Etch  Solution 


0.324  grams  I  (solid) 

+  1.134  grams  KI 

50  milliliters  water 

Chromium  Etch  Solution 
3.25  grams  KMnO^ 

3.25  grams  NaOH 
50  milliliters  water 

9.  Etch  away  exposed  gold  and  chrome  metallizations. 

10.  Strip  remaining  photoresist  with  Allied  Chemical 
A-20,  a  phenol  based  industrial  stripper. 

D.  Housing  Construction 

1.  Design  housing 

2.  Machine  housing 

3 .  Apply  Ni  flash 

E .  Final  Delay  Line  Assembly 

1.  Assemble  coaxial  connectors 

2.  Gapweld  Au  ribbon  to  antennae  input  and  output 

3.  Screw  front  pole  piece  into  main  housing,  ad¬ 
justing  height  so  that  substrate  surface  is 
5-10  mils  below  connector  height. 

4.  Attach  substrate,  antennae  side  up,  to  front  pole 
piece  with  Ag  point.  Antennae  input  and  output 
must  line  up  with  connectors. 

5.  Attach  Au  ribbon  to  connectors  with  Ac  epoxy 
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6.  Assemble  balance  of  delay  line,  including  back- 
pole  piece,  YIG  rod,  teflon  rod  retainer,  fixture 


retainer. 
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Appendix  A  -  Pole  Piece  Design 


I.  Overview 


This  appendix  outlines  the  synthesis  procedure  de¬ 
veloped  by  Moraenthaler  and  Platzker  [19].  Although  more 
general  cases  can  be  treated  using  similar  methods,  this 
discussion  will  be  restricted  to  the  synthesis  of  cylin- 
drically  symmetric  fields  where  the  material  shape  is  a 
cylinder  or  disk.  The  logic  behind  the  pole  piece  design 
is  as  follows.  The  magnetic  potential  inside  the  cylin¬ 
drical  material  is  determined  by  an  appropriate  series 

expansion  once  the  axial  magnetic  field,  H  (0,z),  is  spec- 

z 

ified.  The  space  outside  the  cylinder  is  then  divided 
into  separate  regions.  V/ithin  each  outer  region,  an  ap¬ 
propriate  selection  of  functions  is  made  with  which  to 
expand  the  Laplacian  magnetic  potential  in  that  region. 

The  coefficients  for  each  outer  region  expansion  are  then 
determined  by  matching  boundary  conditions  with  the  inner 
potential  so  that  the  error  is  minimized  in  a  least  sauares 
sense . 

Once  the  outer  potential  has  been  determined,  the 
edges  of  pole  pieces  which  surround  the  cylinder  can  be 
designed  to  coincide  with  lines  of  equipotential .  With 
the  application  of  an  appropriate  dc  magnetic  field,  the 
pole  pieces  become  energized  so  as  to  produce  the  correct 
outer  field  and  consequently  the  correct  axial  magnetic 
field  in  the  cylinder.  It  is  important  to  realize  that 
although  an  infinite  number  of  combinations  of  material 
size,  shape,  and  pole  piece  design  exist,  all  of  which 
would  create  the  desired  field,  once  one  has  been  chosen 
the  field  it  creates  is  unique  [18] . 


II.  Inner  Field  Determination 


The  equations  governing  the  qc  H  field  in  a  ferro- 
magnet  that  is  everywhere  locally  saturated  are 


M*  <H  +  H3)  =  0 

(A . 1 .  a) 

M  •  H  >  0 

(A.l.b) 

|  M  J  =  M 

(A. 2) 

V  x  H  =  0 

(A. 3) 

V  •  (H  +  M)  =  0 

(A. 4) 

where  Ha  is  an  effective  magnetic  anisotropy 

field  which  may 

be  ignored  initially.  In  terms  of  the  inner 
potential  ^ ,  the  relations 

scalar  magnetic 

H  «  -  Vtfr 

(A. 5) 

and  M  =  -  ^  •  M 

(A. 6) 

satisfy  equations  (A.l),  (A. 2),  and  (A. 3).  These  equations 

may  be  substituted  into  equation  (A. 4)  to  arrive  at 

»  •  Id  +,,£,>*•;]  -  o 

(A. 7) 

In  the  limit  as  Jv^|  the  non-linear  portion  of 

equation  (A. 7)  goes  to  zero  and  the  equation  becomes  Laplace’s 
equation.  Thus  if  the  magnitude  of  H  grows  as  distance 
from  the  cylindrical  axis  increases,  then  at  the  radial 
boundary  the  inner  field  potential  should  be  approximately 
Laplacian.  This  aids  the  matching  of  boundary  conditions 
with  the  outer  field,  which  is  Laplacian. 


■■V 
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The  axis  of  the  cylindrically  symmetric  ferromagnet 
is  assumed  nonsingular,  so  an  appropriate  expansion  of  .  is 

00  2 

-  \l>  =  E  a2n(z)  r  (A. 8) 

n=0 
00 

so  that  H(r,z)  =  i  la'  (z)  r2n 

2  n=0"  -  . 

+  *r  =  2(n+1)  a2(n+l)(z)  r  <A*9> 

where  in  is  the  unit  vector  in  the  n  direction  and  the 
primes  denote  differentiation  with  respect  to  z.  Expanding 
M  separately  as  » 

M(r,z)  =  M [i  E  b„  ( z )  r2n 

n=0  oo  2n  , 

+  x  E  b2n+l(z)  r"n  ]  (A. 10) 

n=0 


allows  the  use  of  equations  (A.l),  (A. 2),  and  (A. 4)  to 

determine  constraints  on  the  a^  and  coefficients. 
Use  of  equation  (A.l. a)  results  in 


M(  a2n(z)  r2n) ( 


OO 

E 


n=0 


.2n, 


n=0  2n+l 


(z)  r2n+1) 


-  M(  j  b2n(z)  r‘“)  (  E  2  (n+1)  a2(n+1)(z)  r2n+1  =  0 
n=0  n=0 

where  M  may  be  factored  out  as  assumed  non-zero.  Multiply¬ 
ing  out  and  collecting  terms  yields  a  polynomial  in  r, 


E 

s=0 


OO 


{[  r 

n=0 
-  2 


(a2n  b2 (s-n) +1 
(s+1~n)  a2  (s+l-n) 


b2n” 


s+1 , 
r  } 


0 


(A. 11) 


Since  this  equation  must  be  satisfied  for  all  r,  we  must 
have  that 


(a2n  b2 ( s-n) +1 

n=0 


2 (s+l-n) 


a2  (s  +  l-n) 


b2n}  =  0  (A'12) 


for  s  =  0,  1,  2,  3 .  This  gives  the  first  constraint 

on  the  coefficients.  However,  equation  (A.l.b)  also  re- 
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stricts  the  si<~n  of  bQ  to  be  bQ  >  0. 

Use  of  equation  (A. 2)  results  in 

(  i  b2n(z)  r2n) 2  +  (  £  b2n+1(z)  r2n+1)2  =  1 

n=0  n=0 

Multiplying  out  and  collecting  terns  yields  a  polynomial 
in  r. 


b0  +  *  [nfn  (b2n+l  b2  (s-n)  +1  +  b2n  b2(s+l-n)} 


Tk  „  l  r  2  ( s+1 )  _  . 

b0  ^(s+l)1  “  1 


(A. 13) 


Setting  r  =  0  reveals  that  bg  =  1.  Furthermore,  from  re¬ 
peatedly  differentiating  the  equation  with  respect  to  r 
and  then  evaluating  at  r  =  0  we  find  that 


Z.  (b2n+l  b2  (s-n) +1  +  b2n  b2(s+l-n)) 
n=0 

+  b0  b2s+l  *  ° 


(A.l-'l 


for  s  =  0,  1,  2,  3,  ...  .  This  gives  the  second  constraint 

on  the  coef f icients . 

Use  of  equation  (A. 4)  results  in 
00 

’  •  J0Ja2„  +  Mb2„>  r2" 

+  5r  E  (M  b2n+l  +  2<n+1>  a2 (n+1 ) 1 
n=0 

00 

or  [4  (s+1)  2  a2(s+1)  +  2  (s  +  1)  M  b2s+1 

+  a2s  +  M  b2SJ  r2S  =  0  <A*15) 
This  can  hold  true  for  all  r  if  and  only  if 


4(S+1)  a2  (s+1)  +  2(s  +  1)  M  b2s+l 


+  a2s  +  M  b2s 


0 


(A. 16) 


X 
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for  s  =  0,  1,  2,  3,  ...  .  This  gives  the  final  constraint 
on  the  coefficients. 

In  addition  to  equations  (A. 12),  (A. 14),  and  (a. 16)  we 
have  the  constraints  that  bQ  =  1  and  bQ  a^  >  0.  Since  a^ 
must  be  either  positive  or  negative  definite  over  the  inter¬ 
val  to  insure  saturation,  without  loss  of  generality  we  may 
take  a^  >  0  and  b^  =  +1.  The  solution  of  equations  (A. 12), 
(A. 14),  and  (A. 16)  is  then: 


s  =  0 , 


(A. 17. a) 


(A. 18 . a) 
(A. 19. a) 


s  >  1, 

b2s+l  =  2  (s+lj  (a2s  +  M  b2s5 


(A. 17.b) 


s 

I 

k=l 


[2  (s-k+1) a 


2 (s-k+1) “2k~a2k"2 (s-k) +1 


1  )/(a’+M) 


a2  (s+1)  {a0  b2s+l 


(A. 18 .b) 


k=i [2  (s_k+1)a2 (s-k+1) b2k“a2kb2 (s-k) +11 }/2 (s  +  1) 


b2 (s+1)  '  fk=lbk  b2s-k+2  +  2  bs  +  l]  (A . 19 .b) 

Of  course,  when  actually  calculating  the  inner  potential  on 
the  computer  the  a.  and  b.  coefficients  are  not  found  expli¬ 
citly,  but  approximated  numerically. 


III.  Outer  Field  Determination 


Now  we  must  find  an  appropriate  expansion  for  the 
Laplacian  magnetic  potential  outside  the  cylinder.  Laplace' 
equation  in  spherical  coordinates  (o,9,0)  is 


a 2 ^  +  _i  dz'i>  +  1  32^  2  at  +  cote  at 

3p 2  p2  3 9  2  p2  sin  2  6  3 g2  p  dp  p2  90 


0 


(A. 20) 


which  has  as  a  solution 


t  =  (C1  pn  +  C2  p"(n  +  1))  (cos  m<*> )  P™(cos  9)  (A. 21) 


If  we  assume  no  t  variation  (i.e.  m  =  0)  and  transform 
this  solution  to  cylindrical  coordinates  (where  of  course 
it  is  still  a  valid  solution)  then  we  have  that 


t  =  ICX  (r2  +  z2)n/2  +  C2  (r2  +  z2)“(n+1)/2]P°(cos  9) 


Ignoring  the  second  solution,  we  finally  arrive  at 

t  =  C1  (r2  +  z2)n/2  P°(cos  9)  (A. 22) 

where  tan  6  =  ^  and  P°(cos  9)  is  the  associated  Legendre 
function  of  degree  n  and  order  zero.  If  we  then  define 

P*  =  (r2  +  z2)n/2  P° (cos  0)  (A. 23) 

n  n 

A 

we  can  expand  the  outer  potential  t  in  terms  of  these  poly¬ 
nomials  as 


-  *  =  C"  V"*’ 


(A. 24) 


For  reference,  the  first  five 


* 

P 

n 


functions  are 


listed 
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★ 

below.  Notice  that  since  the  functions  are 
for  all  2  (i.e.  have  their  sources  at  infinity 
be  used  for  regions  outside  the  cylinder  which 
z-axis . 


non-sineular 
,  they  r.ust 
contain  the 


* 


In  addition,  since  in  general  the  desired  on  axis  field 
to  be  synthesized  is  a  polynomial  expression,  and  since 
the  inner  potential  was  expanded  in  a  polynomial  in  r  and 
z,  it  makes  sense  to  expand  the  outer  potential  in  the  above 
polynomial  expressions  rather  than  the  usual  Bessel  functions. 

In  regions  outside  the  cylinder  which  do  not  include 
the  z-axis  (or  where  the  z-axis  is  singular) ,  solutions  with 
a  singularity  at  r  =  0  may  (or  must)  be  added  to  equation 
(A. 24).  Instead  of  using  the  second  solution  to  Laplace's 
equation  in  spherical  coordinates  which  would  yield  a  func¬ 
tion  singular  only  at  a  point,  equation  (A. 22)  is  used  along 
with  Laplace's  equation  (ignoring  0  variation)  in  cylindri¬ 
cal  coordinates, 

ill  +  k  +  ill  =  o  (a. 25) 

3r2  r  3r  9z2 

to  obtain  a  second  set  of  solutions  singular  on  a  line  of 
the  form 

Q*(r,z)  =  P*(r,z)  •  In  r  +  N*(r,z)  (A. 26) 

n  n 


-44- 


where 


N*  =  zn  +  f  i  (2  (2n-P  * )  -  n(n-l)  •  \  r  N*  dr] 
n  J  r  n  J  n-2 


Q  (r,z)  is  more  usefully  given  by  the  recurrence  formula 


*  *  n-1  2  2  *  2  *  * 

0  =  =£-=■  2  Q  .  -  —  (z  +  r  )Q  0  -  -  P  -zP  .) 

vn  n  n-1  n  n-2  n  n  n-1 


(A. 27) 


* 

with  Qq  =  (In  r  +  1) .  Therefore,  in  general  the  outer 
potential  y  may  be  expanded  as 

N  *  N-1 

-  i{*  =  Cn  +  Z  C  P(r,z)  +  Z  C '  Q  (r ,  z )  (A. 28) 

0  n  n  n  n 


where  represents  a  reference  potential  and  nay  be  set 
to  zero. 

★ 

For  reference,  the  first  five  Q  functions  are  listed 

*  n  * 

below.  Notice  that  the  Q  functions,  unlike  the  P^  func- 

n  n 

tions ,  have  their  sources  at  both  zero  and  infinity. 


Q*  =  In  r  +  1 

* 

Q2  =  z  In  r  +  z 

*  2  12  2 
Q3=z  In  r  -  j  r  lnr+z 

*  3  3  2  3 

Q^  =  z  In  r  '  j  z  r  lnr+z 

*4  22  34  434 

Q^  =  z  lnr  -  3  z  r  lnr+g-r  lnr+z  -  jg-  r 

*  * 

Although  the  functions  can  be  approximated  by  the  P^ 

functions  in  many  cases,  this  would  require  many  terms  in 

* 

the  expansion  and  P^  functions  with  large  powers  of  z. 

Since  in  general  it  is  better  to  expand  in  lower  order 

polynomials,  if  possible  0n  functions  should  be  included 

in  the  expansion  of  the  outer  potential.  Expansions  using 
*  * 
thf  first  five  P^  functions  and  the  first  five  functions 

yielded  results  consistent  to  better  than  an  Oersted  for  the 

cases  examined  in  this  thesis. 


IV.  Boundary  Conditions 


This  discussion  is  restricted  to  the  case  of  a  cylinder 
or  disk,  so  it  is  beneficial  to  subdivide  the  outer  region 
at  the  corners  and  carry  out  three  separate  expansions  for 

A 

<i/.  Those  outer  regions  with  z  <  0  and  z  >  L  cannot  con- 

*  am  3  •') 

tain  Q  functions.  The  matching  of  p  =  b  and  (1+=— )  ^  = 
n  J  3n 

is  carried  out  in  a  least  squares  sense.  That  is,  if  we 
assume  the  outer  potential  in  a  given  region  has  been  spec¬ 
ified  as 


A 


N 

£  C.  f  .  (r ,  z } 
j=l  3  3 


(A. 28) 


then  we  are  trying  to  minimize  the  least  square  error  de¬ 
fined  as 


E  = 


X 


Boundary 

Surface 


i*  - 


N 

Z 

j  =  l 


C..  f  (r , z ) ]2  dS 


(A. 29. a) 


An  equivalent  expression  is 


E  =  <Af  (ip  - 


N 

Z  C  .  f  .  ) 

j-i  3  3 


>  + 


<i2"1+W 


3f  . 


(A. 29 . b) 


where  <  >  denotes  averages  over  the  boundary  radius,  and 


A^  and  weight,  respectively,  the  relative  importance  of 


the  tangential  and  normal  components  of  the  field  subject 
2  2  SE 

to  A.  +  A_  =  1.  Requiring  =  0  yields  the  set  of  equations 

L  Z  ~  dC^ 

N 
Z 


3-1 


2  2  9fi  9fi 

[<a:  f.  f.>  +  <\t  ^-i>)  c. 

1  1  j  2  3n  3n  3 


=  <X1  fi 


ip  >  +  <  A 


2  9fi 


2  3n 


(1  + 


JL) 

V'p  3  n 


(A. 30) 


for  i  =  1,  2,  3,  ...,  N.  The  solution  of  equation  (A. 30) 


for  the  C.  coefficients  vields  an  outer  field  that  satisfies 
J 

the  boundary  conditions. 

The  right  circular  cylindrical  shape  of  the  samples 
considered  can  be  exploited  to  simplify  the  solution  of 
equation  (A. 30).  The  radial  match  at  r  =  R,  and  the  end 
match  at  z  =  L  are  considered  separately. 

IV. A.  Radial  Match 


In  order  to  obtain  a  good  match  with  the  outer  field 
at  the  radius  of  the  cylinder,  the  inner  field  is  first  ap¬ 
proximated  as 


hz(2,r  ) 


M-l 

I 

n=0 


a 


n 


(£)n 


(A. 31) 


and 


[Hr (z,R  )  +  Mr (Z,R  ) ] 


N-l 

Z 

n=0 


n 


(A. 32) 


where  N  is  an  integer  yielding  tolerable  error.  This 

approximation  is  necessary  since  the  inner  field  at  the 

radius  no  longer  has  a  polynomial  form,  while  the  outer 

field  at  the  radius  does  have  a  polynomial  form. 

The  a  and  6  coefficients  are  found  from  a  least 
n  n 

squares  fit  using  an  inverse  Hilbert  matrix  of  order  N. 

This  arises  from  solving  for  the  least  squares  coefficients 
such  that  the  inner  potential  expansion, 

N 

f  (x)  =  Z  A.  g.  (x)  (A.  33) 

k=l  K  * 

is  approximated  as  a  polynomial 

N 

f°(x)  =  E  aP  xk_1  (A. 34) 

k=l  K 


The  matrix  least  squares  equation  is  of  the  form 
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<g  l> 


'1  y2' 


<g2  V  <g2> 


<gl  93" 


<g2  g3> 


<A°> 

<*2> 


glf> 


<92f> 


(A. 35) 


L<gN 

gl> 

<9N93>... 

A 

> 

z  c 

V 

A 

V 

mJ  \ _  — J 

L_  — 2 

where 

<gi 

V  -  sf  0  gi(x)  gj 

(x)  dx 

and 

<gi 

f>  =  tSo  9i(x)  f<x> 

dx 

If  the  length  of  the  cylinder  is  normalized  to  unity  then 


<g .  g . > 
1  3 


-li*1 


i+ j  -2 


dx  = 


ITFT 


Substituting  into  equation  (A. 35)  yields 


!  1  i  I 

2  3  4 

VI 

I  f  { :•; )  dx 

1111 

2  3  4  5 

*2 

I  x  f (x)  dx 

1111 

3  4  5  6  •  •  ' 

•  •  •  • 

*3 

• 

= 

I x2  f (x)  dx 

•  •  •  • 

i 

A;i 

Jx'1  f.  (x)  dx 

_  — w 

—  — 

(A. 36) 


where  the  matrix  multiplying  the  A°  coefficients  is  the 
Hilbert  matrix  of  order  N. 

The  solution  to  equation  (A. 36)  is  then 


I  f (x)  dx 

=  h”1 

i  x  f (x)  dx 

N 

* 

• 

S  x^  f (x)  dx 

_  _ 

_ 

(A. 37) 
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where  is  the  inverse  Hilbert  rutri.x  of  order  >; .  However, 

the  accuracy  obtained  with  a  comouter  solution  to  enuation 
(A.  37)  can  be  deceptive.  The  Hilbert  matrix  is  an  ill-cor.Ji- 
tioned  matrix,  whose  inverse  rapidly  overflows  the  capacity 
of  a  digital  computer.  Fortunately,  from  an  engineering  stand¬ 
point  this  inaccuracy  is  largely  unimportant,  changing  the 
pole  piece  design  only  at  distances  far  from  the  crystal  sam¬ 
ple.  Since  solutions  to  Maxwell's  equations  typically  decay 
exponentially,  varying  the  pole  piece  design  far  from  the 

sample  has  little  effect  on  the  axial  magnetic  field.  There- 

1 

fore,  the  ill-conditioned  property  of  the  Hilbert  matrix  ac¬ 
tually  allows  greater  flexibility  in  the  pole  piece  design 
and  implementation. 

In  terms  of  «  and  g  ,  the  C  and  C'  coefficients  of 
n  ~n  n  n 

equation  (A. 28)  are  found  to  satisfy 


C  = 


an— 1 

-grr-  -  (InR  +  1)C' 
P  D  P 


*  (P  ih  D 
N~p  2  k 

-  R  Vu  •  Vx  +  (lnR+*k)Cp^] 

c  -  Lje  -  Nep  R2k  a 

P  LP  k=l  p+k'k 


[A. 37. a) 


{2klcp+k  +  (In  R  +4>k)  Cp+k]  +  C’+k>  (A.  37  ,b) 


where 


An,k  =  ^ 


n: 


(n-k) ! [ (|) ! ]2  2k 


and 


k=0 


1  - 


-i+L 

2  3 


■s>  k>0 


These  equations  are  not  general,  but  hold  only  for  the  special 
case  of  a  right  circular  cylinder  under  discussion  Since 
by  design  =  0,  equation  (A.  37. a)  immediately  gives  CfJ. 
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Alternating  between  equations  (A.37.b)  and  (A. 37. a)  produces 
an  unraveling  that  calculates  in  order  C^;  ->•  ; 
. ..;  Cq.  The  value  of  is  immaterial  and  may  be  set  to 

zero. 


IV. B.  End  Match 


The  outer  potentials  for  z  <  0  and  z  >  L  can  be  found 
in  a  similar  fashion.  For  the  z  =  0  end  face,  matching  y 
and  Hz  +  Mz  yields 


C 


n 


r 

< 

v. 


an<°> 

A 

n ,  n 


a  .  (0)  +  M  b  .  (0) 
n-l  n-i 


n,n-l 


n  even 


(A. 38) 


n  odd 


where  a  and  b  satisfy  equations  (A. 17),  (A. 18),  and  (A. 19). 

n  n 

Similar  expressions  result  for  the  z  =  L  end  face. 


V.  Design  Example 

To  illustrate  the  use  of  the  synthesis  technique,  it  is 
helpful  to  work  through  a  low  order  example  by  hand.  We  shall 
attempt  the  synthesis  of  a  linear  axial  magnetic  field  of  the 
form 


H  (0,Z)  =  A  +  B*  z  (A. 39) 

z 

in  a  YIG  cylinder  of  length  L  and  radius  R.  For  simplicity 
anisotropy  will  be  ignored.  Choosing 

a q  =  A»z  +  j  B.z^  (A. 40) 

satisfies  the  requirement  that 


-so- 


7  <  -  Jo  a2"U>  r2n  Mr-0  *  A  +  B'z 


(A. 41) 


Using  equations  (A. 17),  (A.18) ,  and  (A. 19)  we  have  that 

1  A 


bl  " 


a2  = 


b2  “ 


2  B  •  2  +  A  +  M 

1  A  (A  +  B  •  2  ) 

4  B  ■  Z  +  A  +  M 


1 

8 


(B  •  Z  +  A  +  M) 


(A. 42 .a) 
(A. 42  .b) 
(A.42.C) 


This  allows  the  inner  magnetic  potential  to  be  approximated 
as 


\l>  -  -  (A  2  +  z^)  + 


1  (  A  (A  +  B  z)  \ 
4  \  B  z  +  A  +  M  J 

The  inner  magnetic  field  is  then  given  by 

A  B 


(A. 43) 


H  (r ,  z)  =  [A  +  Bz  -  4  M  z  +  A  +  M) 


r2  ]  t. 


+  [  - 


1  A  B  z  +  A' 


M(r,z)  =  M  [  1  - 


:  B  z  + 

A  +  M  r 

by 

1 

A2 

8  (B  z 

+  A  +M) 2 

1  A  (A 

+  B  z) 

: — i - : — —  r 

(A. 44) 


2  •.  A 

r  ]  l. 


(A. 45) 

i-j  i*  i  r\  ~r  r*  X. 

Using  equations  (A. 31)  and  (A. 32)  we  then  have  that 

(A .  4  6  .  a ) 

[  Hr(R,z)  +  Mr(R,z)  J  *  60  +  (£)  +  B2  (£)2 


Hz(R,z)  *  aQ  +  (£)  +  a2  (£) 2 


(A. 46 .b) 
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For  simplicity  of  calculation,  we  now  assume  that  we  are 
only  concerned  with  matching  tangential  H,  and  we  ignore  the 
matching  of  normal  B.  This  corresponds  to  setting  =  1  and 
*2  =  0  in  equation  (A. 30).  Expanding  the  outer  potential  as 

2  *  *  * 

~  -  E  C  P  (r ,  z)  +  C'  Qn(r,z)  +  Cl  Q,  (r,z)  (A. 47) 

Ann  uu  i.1 

n=0 

allows  the  least  squares  coefficients  to  be  found.  In  terms 

we  find  that: 

{if1  ?  +  lnR> 


°i) 


More  complete  solutions  can  be  obtained  using  the  com¬ 
puter  methods  described  in  section  IV.  For  computer  solutions 
an  expansion  of  the  form  of  equation  (A. 20)  is  assumed.  For 
the  case  of  A  =  B  =  300  in  equation  (A. 39)  we  have  as  the 
coefficients : 


of  the  a  coefficients 
n 


C„  = 


°0  "  a2 


C,  =  ~ 


1  al 


2  L 


C„  «  * 


1  a2 


C ' 

‘“O 


C' 

^1 


l(°_ 1 
2  V  L 


■  ■>  (f/ 


End  Match: 


Coefficient 

Value 

C1 

-  2160.0000 

C2 

-  21.6346 

C3 

-  9.2571 

C4 

-  0.2525 

0.1111 

5 
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Radial  Match: 


A  plot  of  the  equipotential  contours  implementing  both  the 
radial  and  end  magnetic  potential  matches  is  shown  in  figure 
(A . 1 ) .  Here,  the  magnetic  potential  end  match  at  z  =  L  lies 
above  that  of  the  radial  match  at  z  "  L,  so  the  area  enclosed 
by  the  two  equipotentials  would  have  to  be  negative.  This  is 
not  physically  realizable,  and  the  end  match  at  z  =  L  is  not 
implemented.  Since  the  magnetoelastic  interactions  occur  near 
z  =  0,  this  does  not  normally  present  any  real  problems. 


Synthesized  Magnetostatic  Modes  and  Waves  Using  Nonuniform  Magnetic  Bias  Fields 
Our  progress  on  employing  nonuniform  dc  field  synthesis  to  magnetostatic 
modes  and/or  waves  is  described  in  the  attached  set  of  preprints  and  reprints 
of  papers  that  were  presented  at  the  following  conferences: 

-1980  Intermag  Conference,  Boston,  Massachusetts,  April  21  -  24,  1980 

-1980  International  Conference  on  Ferrites,  Kyoto,  Japan,  September  29  - 
October  2,  1980. 

-1980  Ultrasonics  Symposium,  Boston,  Massachusetts,  November  5-7,  1980 

-1980  Conference  on  Magnetism  and  Magnetic  Materials,  Dallas,  Texas, 

November  11  -  14,  1980 

For  completeness,  reprints  are  also  included  of  our  earlier  work  described  at  the 
1977  Conference  on  Magnetism  and  Magnetic  Materials  and  the  1978  Intermag, 
Florence,  Italy. 
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C.  R.  Morgentiialcr 

M.I.T. 

Cambridge,  Massachusetts 


Abstract 


V!e  discuss  the  manner  in  which  gradients  in  either  the  bias  field  mauni t..;,' . , 
direction  or  both  can  be  employed  to  synthesize  magnetostatic  wave  dispersion 
characteristics  or  mode  spectra.  This  is  done  to  control  pre- specif ied  character¬ 
istics  such  as  frequency,  rf  energy  distribution,  impedance,  velocity  of  energy 
circulation  and  the  threshold  governing  the  onset  of  nonlinear  effects  due  to 
parametrical ly-unstable  spin  waves. 

Previous  theoretical  models  are  reviewed,  extended,  and  utilized  to  predict 
the  behavior  or  a  new  class  of  microwave  resonators.  In  addition,  we  predict 
that  frequency  selective  filter-iimiters  can  be  constructed  with  gradient- 
controlled  limiting  levels. 

We  also  review  experimental  evidence  of  high-Q  gradient-localized  resonances 
in  platelets  and  thin  films  of  single  crystal  YIG.  Included  are  recent  observa¬ 
tions  of  magnetostatic  surface  waves  propagating  in  a  thin  film  bint  ed  with  a 


nonunifonn  in-plane  field  caused  by  movable  permalloy  strips. 
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and  the  Research  Laboratory  of  Electronics 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts ,  USA 


Introduction 

Our  interest  in  controlling  magnetostatic  waves  and  inodes  by  means  of  dc 
field  gradients  dates  from  the  experimental  observation  of  localized  high-Q 
resonance  in  single  crystal  yttrium  iron  garnet  (V1G)  reported  by  Zeskind  and 
Morgenthaler J 

Our  interpretation  of  these  resonances  is  that  localized  magnetic  mode 
•patterns  are  formed  in  which  the  resonant  energies  are  guided  or  confined  by 
regions  of  high  dc  field  gradient  within  the  crystal.  By  appropriately 
designing  internal  magnetic  field  profiles,  one  can  create  apparent  "surfaces" 
or  "tracks"  of  magnetic  field  discontinuity.  Magnetic  waves  bound  or  guided 
by  such  gradients  can  be  made  to  follow  appropriate  propagation  paths  with 
controlled  group  velcoity.  If  the  mode  amplitudes  are  very  small  at  the  edges 
and  corners  of  the  sample,  the  surface  scattering  (which  one  would  expect  to 
be  enormous)  is  largely  prevented.  Consequently,  the  Q  of  the  resonance 
governed  primarily  by  the  intrinsic.  1  inewidth  of  the  bulk  crystal  together 
with  normal  circuit  loading  considerations. 

In  the  first  experiments  such  gradients  arose  naturally  from  the  nonuniform 

2  3 

shape  demagnetizing  fields;  in  subsequent  work  by  Cooley  and  Horowitz  ,  we 
have  created  pre-speci fied  gradients  with  shaped  pole  pieces  designed  by  field 
synthesis  techniques.  A  similar  approach  has  recently  been  followed  by  Tsutsumi 

at  al 

Wo  first  review  these  early  experiments  and  the  theoretical  treatment 
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quasi-two-dimensional  magnetostatic  modes  of  single  domain  thin  ferrite  cir¬ 
cular  disks  or  annular  rings,  when  the  dc  magnetic  field  is  normal  to  the 

plane  and  varies  radially. 

We  find  solid  theoretical  justification  for  the  'track"  concept  and 

enumerate  mode  patterns  for  both  gradient-modified  boundary  modes  and  gradient- 
5 

dominated  waves  .  The  important  effects  from  rf  fringing  fields  at  the  edge  of 
the  disk  are  also  considered. 

In  the  absence  of  a  radial  gradient,  all  of  the  modes  of  a  solid  disk 
have  circularly-polarized  rf  h-fields  with  zero  volume-divergence.  These  modes 
are  stongly  influenced  by  the  magnetic  pole  distribution  on  the  edge  of  the 
disk  and  the  rf  energy  becomes  progressively  concentrated  near  the  rim  as  the 
mode  index  increases. 

When  the  field  gradient  is  modest,  the  modes  retain  these  general  char¬ 
acteristics  but  develop  non-zero  values  of  v-h  throughout  the  volume  which 
change  both  the  state  of  polarization  of  the  field  and  the  distribution  of  rf 
energy,  the  latter  can  then  predominate  either  in  the  central  portion  or  near 
the  edge  of  the  disk.  In  addition,  both  the  mode  frequency  and  the  velocity 
with  which  the  mode  energy  circulates  are  found  to  be  altered. 

For  increased  gradient  strength,  the  volume  divergence  of  certain  modes 
can  change  so  dramatically  that  selective  localization  or  expulsion  of  the 
energy  occurs.  The  sense  of  polarization  can  also  actually  reverse.  I:'.  cer¬ 
tain  cases,  the  volume  divergence  of  the  rf  magnetization  can  become  infinite 
(in  the  lossless  exchangeless  approximation)  at  a  certain  interior  radius,  r^  . 
The  magnetic  pole  distribution  at  this  "virtual-surface"  thus  resemble  that  of 
a  true  surface  and  can  serve  to  guide  and  localize  the  mode. 

These  previously  discussed  "virtual-surface"  modes,  are  reanalyzed  in 
terms  of  the  polarization  factor  of  the  rf  magnetization.  Analytic  solutions 
of  that  factor  are  continuous  through  the  "vi rtual -surface"  and  are  given  for 
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the  two-dimensional  modes  associated  with  a  special  class  of  field  profiles. 

We  next  discuss  the  synthesis  of  modes  with  pre-specified  characteristics, 
such  as  velocity  of  energy  circulation  and  rf  energy  distribution.  In  addi¬ 
tion,  boundary  conditions  that  have  been  imposed  to  make  the  mathematical  analy¬ 
sis  more  tractable  (namely  placing  the  thin  film  or  disk  between  perfectly 
conducting  plates)  are  removed;  the  bias  field  is  still  normal  to  the  plane  but 
radial  gradients  although  now  generalized  to  arbitrary  form  must  be  only  weak  to 
moderate. 

For  mathematical  convenience,  we  choose  to  model  the  thin  film  disk  as  a 
very  oblate  spheroid  with  semi-axes  a  and  b.  The  thickness  of  the  "disk"  at 
its  center  is  2b  and  the  radius  is  a  where  b«a. 

The  modes  analyzed  are  quasi-two-dimensional  in  that  the  rf  magnetization 
is  assumed  to  be  without  appreciable  thickness  variation;  the  fringing  magnetic 
fields  are,  of  course,  three-dimensional  because  the  spheriod  is  surrounded  by 
free  space. 

The  results  can  be  utilized  to  predict  the  behavior  of  a  new  class  of  micro- 
wave  resonator.  In  addition,  we  predict  that  frequency  selective  filter-limiters 
can  tie  constructed  with  gradient-controlled  limiting  levels. 

Finally,  we  consider  the  guidinq  of  magnetostatic  surface  waves  (MSSW)  on 
films  with  an  in-plane  bias  field  that  may  be  a  function  of  those  coordinates 
transverse  tc  the  propagation  direction.  In  particular,  the  recent  experimental 
observation^  of  MSSWs  in  a  rectangular  YIC  film  place  between  strips  of  permalloy 
and  in  the  plane  of  the  straps  is  reviewed. 

Localized  '■  [i  a  h-Q  Res  or  iancc  i  n  Me  nun  i  form  DC  Mel  ds 

While  measuring  the  microwave  coupling  between  two  antennas  closely  spaced 
on  the  surface  of  a  bulk  single  crystal  YIG  slab  (with  dc  magnetic  field  applied 
perpendicular  to  t.he  slab).  Surprisingly,  despite  the  very  nonuniform  demagnetiz¬ 
ing  field  caused  by  the  nonel 1 ipsoidal  sample  geometry,  Zeskind^  observed  spatially 


localized  hinh-0  resonant  nodes.  The  convenient  planar  geometry  alone  with 
the  hieh-0  and  rianuotic  tunability  of  these  nodes  suggested  the  possibility  ot 
a  totally  intearated  monolithic  YIG  filter. 

The  basic  experimental  configuration  annoars  in  Figure  1.  Gold  thin  film 
antennas  were  photographically  etched  on  the  surface  of  the  bottom  alumina 
substrate.  The  30  um  wide  and  4.2  inm  lonq  antennas  are  separated  by  an  80  pm 
wide  around  strip  to  improve  electro-magnetic  isolation.  The  bulk  single 
crystal  YIG  slab  cut  along  the  (110)  plane  (shown  in  place  on  top  of  the 
antennas)  is  glued  into  a  slot  milled  in  the  surface  of  the  top  alumina  substrate. 
The  top  substrate  and  YIG,  toqether,  move  with  respect  to  the  bottom  substrate 
(antennas)  thus  allowinq  plots  of  spatially  dependent  microwave  properties  of 
the  YIG.  Both  sides  of  the  slab  are  polished  to  a  surface  finish  of  about  one 
half  micrometer.  The  dc  magnetic  field  is  applied  oerpendicular  to  the  plane 
of  the  slab. 

Figure  2  is  a  plot  of  microwave  couni  i  ng  between  the  two  antennas  as  a 
function  of  frequency  with  antennas  located  along  the  center  line  of  the  slab 

and  H,  -  2300  Oe.  The  sharn  resonance  at  2.4f  GHz  is  characterized  bv  a  half- 

dc 

power  bandwidth  of  1.8  MHz  corresponding  to  a  loaded-0  of  approximately  1400. 
Insertion  loss  at  resonance  was  initially  on  the  order  of  lOdB  but  was 
subsequently  reduced  to  under  20b.  The  resonance  is  magnetically  tunable  over 
the  range  2.2  to  3.0  GHz.  No  apparent  power  limiting  was  observed  with  input 
signals  of  up  to  +10  dBm. 

In  Figure  3  we  plot  microwave  coupling  as  a  function  of  antennas  displace¬ 
ment  along  the  surface  of  the  YIG  slab  for  fixed  frequency  (2.4G  GHz)  and  fixed 
field  (2300  Oe).  A  displacement  of  0.0  mm  corresponds  to  antennas  located  along 
the  center  line  of  the  slab.  It  is  evidently  that  spatially  localized  points 
of  resonance  exhibits  symmetry  about  the  middle  of  the  crystal  with  maximum 


-60- 


response  occurring  at  the  midplane.  Half-power  spatial  linewidths  are 
approximately  60  mu. 

Formulation  of  the  Basic  Two-Dimensional  Equations 

We  consider  a  ferrimagnetic  film  disk  magnetized  to  saturation  in  the  z- 
direction  by  a  dc  field  H;  radial  components  are  assumed  negligible  through¬ 
out  the  disk.  Since  we  are  interested  in  magnetic  modes  of  frequency  u>  that 
have  negligible  electric  field  energy,  the  magnetostatic  approximation  is 
suitable  and  the  rf  magnetic  field  may  be  expressed  as  h  =  -Vit.  Because  the 
ferrite  is  assumed  to  be  saturated  by  a  cyl indrical ly,  symmetric  field  H  (r) 
and  modes  without  z-variation  are  sought,  we  take  the  complex  magnetostatic 
potential  (with  suppressed  exp  [jut]  variation)  to  be  of  the  form 


i p  =  exp  ( 1 m | 


' y(r)dr 
r 


3N>) 


(1) 


where  ;r.  is  a  positive  or  negative  integer  and  y(r)  determines  the  radial  varia¬ 
tion  of  the  mode.  The  latter  also  describes  the  polarization  of  the  rf  F-field 

Vhg  =  j  ^  y(r)  {?) 

The  other  field  quantities  follow  from  the  Polder  susceptibility  tensor 


and  the  Maxwell  Equation  governing  the  curl  of  the  electric  e-field. 
They  are 


"V  *  -  (xy  -  * 

(3a) 

■V  3?(x-^y)^ 

(3b) 

and 

ez  =  -wp0|ni|y°(r)«r 

(4) 

where 

y°  =  (l+x)  y  -k  ]S|- 
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Here  x  =  Z/(Z^-sT^)»  k  =  with  Z(r)  =  (r)/M  and  f1  =  The  latter 

are,  respectively,  field  and  frequency  normal  i zed  to  the  saturation  rMqo'-H/itinn 
M;  =  -Tii  M  where  7  (negative)  is  the  gyromagnetic  ratio.  It  is  also  convenient 
to  take  S>0  and  define  a  =  Sim/ 1 m | . 

The  polarization  factors  of  the  rf  magnetization  and  rf  flux  were  ignored 
in  our  previous  work  but  have  been  found  to  be  very  useful  quantities.  Consequently 
we  define 

Vm<|>  =  j  fif  p(r)  and  br/by  =  j  ]¥T  P°(r)  (5a)’  (5b) 

where  from  Eqs.  (3a,  3b)  and  (4) 


Zy-t-o  Zy°+g 
p  Z+oy  Z+l+oy°" 


(Z-t-1  )p-o 
Z+l -op 


(6a),  (6b) 


The  wave  inpendances  are  defined  by 


wuory°(r) 


(7a) 


(3b) 
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and 


do  m  / ,  2,  ''(n  -1  )+p_  dZ 

"  I-'  ' 1  “+  >  2  V  rlv 


dr 


,2.7  2  dr 

L  +7.-0 


(Be) 


In  regions  where  dZ/dr-O,  y  and  p  are  each  of  the  form  (r^ ^ -C)/r^ 'm ! +C) 
where  C  is  a  constant;  in  air  regions  y°  =  y. 

In  regions  where  dZ/dr  is  very  large,  the  first  term  on  the  RHS  may  be 

neglected  and  Eg.  (8c)  integrated  directly. 

Bounda ry  Cond it ions  at  the  Edge  of  the  Parallel  Plate  Region 

As  a  concrete  example,  consider  a  thin  ferrite  disk  of  thickness  d  and 

radius  R  placed  between  perfectly  conducting  circular  plates  that  extend  to 

the  radius  R  .  The  cylindrical  surface  r=R  is  assumed  either  conducting  for 
o  J  o 

| z |  <  d/2,  (a)  or  for  \z\  > d/2,  (b). 

For  (a),  nr=0  at  the  radius  r=RQ  and  from  Eq.  (7a),  y°(R  )=0. 

For  (b)  a  reasonable  approximation  valid  when  d/R  <<1  is 

0 


[y°(",Rjr 


1  m  |  d  .  , 


2v  R  T 

0.,.,/  o)  o.„2 

r4ltZ/T  +~d~Ti" 

ZJmJ+l 


(9) 


Solutions  When  Z(r)  =  A- (A+Br2n)/(A^+Br2n) 

In  order  to  develoo  insight  into  the  character  of  the  modes  that  result 
when  Z(r)  forms  a  ootcrtial  v/e 11,  we  consider  this  general  class  of  profile 
for  which  Eg.  (8c)  reduces  to 


~  JuL 

dZ  2n 


(A  -A  )  o(p'  ^  <p+  T-) 
00  o  ±  x 

\7-KJ0r-  V )  TZ-Zxl(Ii?x+l) 


(10) 
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Once  a  solution  of  Eg.  12  is  found,  it  can  be  applied  to  profiles  with  any 
desi red  value  of  B. 

0_  r 

When  this  family  reduces  to  Z=A+Br  analyzed  previously. 

Sketches  of  possible  rf  magnetic  field  contours  are  given  in  Fig.  3  of 
Ref.  5. 

For  weak  gradients,  |BR^nj«l  and  non-"vi rtual-surface"  modes  exist  for 
m>0  when 


A  + 


m 

n+m 


BR' 


2n  ,  1 


,  l+y°(R0)  R  2m 

r  [1  +  T-THf^T  ] 


OD 


Solutions  for  p(Z) 

A  series  solution  for  p(Z)  valid  in  the  vicinity  of  Z=ZQ  can  be  obtained  by 
assuming 


p(Z)  =  5:  C.(Z-Z  )K 
k-0  K  0 


M2) 


and  using  Eq.  (13)  to  find  the  recurrence  relationship  among  the  various 
coefficients,  CK,  given  that  Co  -  p(ZQ). 

The  special  cases  Z  =AQ  or  A^  must  be  handled  separately  and  the  series 
representation  generally  fails  when  7 ( r  )  -  Z  =  (/H4  M-l  )/2  because  at  a 

X  X 

"virtual -surface"  y(r  )  is  singular  and  the  radius  of  convergence  for  the  power 
series  goes  to  zero. 

On  the  other  hand,  from  Eq.  (10)  it  follows  that  p(Z  )  =  Z Jo  (or  possibly 
-o/Z^)  and  the  finite  value  allows  one  to  integrate  p(Z)  through  the  virtual 
surface. 

The  result  of  a  detailed  analysis  reveals  that  with  u  =  (Z-Zx)/(Zx-A  ) 


p(u)  = 


7 

_x_ 

o 


u  dW  ,  r  a  }mj 
w  du  '  LTi-Zx+l  "  2n 


(A  -A  ) (7-Z  ) 
o  x  -] 

TzmTTCa  -  z)  j 


i  j 


) 


where  for  A  satisfies  the  follcwinq  form  of  linear  equation  with 

nonconstant  corf f i i ients : 
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r,  0  r  [> 

n,  v  o  Ki  d  r  v  „  '  dV,'  ,  ,  n  K.  W 

K- 1  K  di/  K=  1  ^  u  du  K--1  K  (/ 


04) 


using  the  method  of  Frobenius,  the  solution  of  Lo  (14)  is  of  the  form 


W  ~  (C.+tr, ! u  i )  b„uK+  T  d  uK  (15) 

h  K=1  K  k=l  K 

where  and  dK  are  constants  related  to  the  RK,  P,,,  0^  coefficients  and  is 
an  inteqration  constant  used  to  match  between  the  power  series  expansions  of  the 
form  of  Eg.  (Ik!)  used  on  either  side  of  Z^.  The  solution  is  valid  when  juj-l 
but  it  is  also  possible  to  expand  W  iri  powers  of  u-^  and  hence  obtain  the 
asymptotic  behavior  for  larqe  Z. 

9  10 

It  should  be  noted  that  exchange  effects  ’  at  the  "virtual -surface"  are 
ignored  in  this  model. 

The  properties  of  an  m=l  "vi rtual -surface"  mode  are  illustrated  with  the 
profile  Z  ~  .3  +  1.095(r/R)^  for  which,  if  y  (m=l,P.)  =  -1,  such  a  mode  occurs 
when  M  -.655  .  The  associated  in-field,  h-field  and  F-field  loci  are  Dlotted, 
respectively,  in  Figure  (4a,b,c)  with  the  "virtual-surface"  radius  (r  =  .3P5R) 
shown  by  the  dotted  circle.  Notice  the  reversed  direction  of  m  leadinq  to 
"surface"  magnetic  poles  near  r-r  .  Notice  too,  that  p(r/R  ••  .59)  =  0.  Along 
with  changes  in  larger  values  of  BR^  will  cause  r^  to  shrink;  smaller  values 

to  expand. 

Synthesized  Magnetostatic  Resonances  in  a  Nonuniformlv  Biased  Thin  Disk  Without 
Conduc Li nq  Bout. Juries 

If  the  conducting  plates  bounding  the  ferrite  are  either  separated  or  removed 
entirely,  the  rf  field  inside  the  disk  will  fringe.  Then  field  variations  with 
respect  to  z  may  be  important  even  for  the  q.;asi-two-dimensional  modes,  as  the 
outer  and  inner  fields  interact  along  the  entire  disk  surface.  If  the  ferrite 
is  comparatively  thick,  z-variations  in  both  H^c,  and  hence  \  ,  will  occur  and 
mode  localization  effects  due  to  these  vertical  gradients  can  be  expected. 


i 


-65- 

When  the  disk  is  thin,  the  gradients  weak  to  moderate,  so  that  the  modes 
in  question  do  not  have  "virtual-surfaces,  another  form  of  analysis  is 
convenient.  We  here  review  that  alternate  approach.^ 

The  complex  ngnetostatic  potential  within  the  ferrite  is  approximately 
given  hy  t,!j=R(r)e  where  r  and  ■*>  are  cylindrical  coordinates  the  factor 
is  again  suppressed  and  m  is  a  positive  mode  integer. 

We  again  choose  to  work  with  dc  fields  and  frequencies  that  are  normalized 
to  the  saturation  magnetization  M  of  the  material.  However,  we  now  Generalize 
the  dc  bias  to  be 

N  ?n 

Hz/M  =  Z(r)  =  A  +  Z  Bnr  n  (16) 

£  n=l 

In  terms  of  circularly-polarized  small-signal  magnetization  vectors 

m  -  =  m+(r)(T+jT,)  (17a,  b) 

and  their  associated  scalar  susceptibilities  x~  -  1/Z+ft  the  equation  7*b=0 
and  Vxfi=0  become 

V  •  [(l+Z-ft)m  +  +  (l+Z+fl)m"]  =  0  (18a) 

V  x  L(Z-!d)m  +  +  (z+$;)m  ]  =  0  (18b) 

Because  the  field  gradients  are  assumed  to  he  modest,  the  mode  frequencies 

for  low  order  m  will  satisfy  |Z-fi|  <<  1  and  m"  will  be  small.  Therefore  m  +  is 
expected  to  be  approximately  Laplacian  in  character  and,  in  fact,  detailed 
perturbation  analysis  reveals  that  when  jm  /m+j<<!  the  associated  potential 
inside  the  spheroid  is 


i 


mo[ 


71-  A 


N 

-  Z 
n=l 


B_n 

n+in  -'a 


(19) 


whereas  outside  it,  the  potential  satisfies  Laplace's  Equation.  Both  and  the 
the  normal  component,  of  5"  are  continuous  at  the  boundary. 
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Tho  application  of  the  boundary  condition  v1-  -a0  over  the  boundary  surface 
is  straiqht-forward,  but  calculation  of  b^=b°  involves  evaluation  of 
3^/3z*0. 

Fortunately,  there  is  an  alternate  wav  to  Droceed  that  circumvents  this 
difficulty;  we  calculate  the  dominant  contribution  to  the  far  field  potential 
from  our  knowledne  of  the  form  of  m  inside  the  spheroid  and  match  it  to  the 
asymptotic  value  of  y;0  .  The  net  result  is  that  a  mode  exists  whenever 


n  =  A  +  N(m) + 


1+3/2 


(C/ 


1 


1+3/2 

m+1 


<v 


1 


1+3/2 

m+2 


(V 


)))..) 


(20) 


where 


N(m) 


b 


(m-1 ) !  2m 


(21) 


and  r-^ — 

Fm(m)U)  =  (2m-l)!!fm  sin'1(J-)  C(2m-2)!!  +  (2m-l )  (2m-4) ! !  r} 

+  ( 2i'.i- 1 )  ( 2i;i-3 )  ( i  iii~6 ) !  !  +  . ..+  (2m- 1 ) !  !  c,'~m  ]  (‘-1-; 

where  n!!  =  n(n-2) (n-4). . . 
which  both  Vi-0  and  b/a-^0. 

2 

In  such  cases,  retaininq  terms  in  h  that  are  first-order  in  kQ  ,  leads  to 
a  modification  of  the  normalised  frequencies  that  can  be  incorporated  by 


lettinq 

(kQa)2 

N(m)  '  fUm)  - 


As  expected,  this  correction  is  la  roost  for  m-1. 
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Synthosi  zed  _Di  spersion  Rel  a  L i  ons 

If  it  is  desired  to  create  a  spectrum  with  the  resonance  frequencies 

separated  by  pre-specif ied  amounts  or  if  one  wishes  to  control  the  velocity 

of  enerav  circulation  v2  of  individual  modes,  the  independent  constants  B 
"  <p  ■  n 

can  be  adjusted  and  the  required  field  H  (r)  synthesized.  In  the  former  case, 

Eq.  (27)  is  used  alone  whereas  in  the  latter,  the  relation  v^(m)- r (m)/'  m 

if 

(analaaous  to  the  croup  velocity  for  a  plane  wave)  is  also  employed. 

It  is  useful  to  realize  that  v2  can  be  forced  to  be  independent  of  m  over 

some  ranqe  of  m.  This  also  suqqests  that  maqnetostatic  plane  wave  propaqation 

can  be  made  precisely  nondispersive  over  a  predetermined  bandwidth. 

Control  of  v2  not  only  affects  the  qroup  delay  of  signals  propagatinq  through 

the  mode  but  also  the  total  energy,  E,  of  the  mode  in  terns  of  the  siqnal  power 

P..  The  Governing  relationship  is  E=2~rPVv2.  The  important  point  to  be 

realized  is  that  the  normally  slow  energy  circulation  that  occurs  when  B  =0  can 

either  be  speeded  up  or  slowed  down.  In  the  latter  event,  the  direction  of  net 

power  flow  can  even  be  reversed.  Near  the  balance  noint  where  V  •  0,  E  becomes 

very  large,  for  a  fixed  value  of  P^.  Because  nonlinear  behavior,  due  to 

1 3 

parametric  spin  wave  instabilities  ,  occurs  when  the  energy  density  of  the 

mode  reaches  a  critical  value,  it  follows  that  the  threshold  power  of  limitinn 
level  should  be  gradient  controllable. 

Consider  two  examples  in  which  both  kQa  and  N(m)  are  negligible. 

In  the  first,  we  require  that  B-j  and  B2  force  rh?/'»m-0  for  m-I  and  .01  for 
m=2.  The  required  values  are 


=  -.267  B2a4  =  .344 


In  the  second,  we  require  that  B-j,  B^  an.1  B^  force  3f)/"sm~.01  for  r 


i~l  and  0 


for  m=2  and  m=3. 

The  required  values  are 


Bja  -  .477 


$2<V  *  -1.121 


B3a  “  .739 
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In  Figures  (5a, b)  and  6a, b)  the  values  of  P.(m)  and  Sfi/am  are  plotted  vs. 
m  for  the  two  cases;  in  Figures  (5c)  and  (6c)  the  required  field  profiles  of 
Z-A  are  plotted  vs.  r/a. 

Field  Synthesis 

In  order  to  normally  magnetize  a  very  thin  disk  of  magnetization  M  so  as 
to  produce  within  it  the  normalized  field  HZ~M  =  A+Br  ,  we  first  create  the 
free  space  dc  magnetic  potential. 

$  =  CQ  +  (A+l )Mz  +  (-1  )n(§HTTI  BMP*2n+l(r’z)  (23) 

where  P*(r,z)  =  (z2+r3)n//2  P°(cosG),  tanG  =  r/z  and  P°  is  the  associated 
n  n  n 

Legendre  function  of  degree  n  and  order  zero.  These  polynomial  satisfy  the 
recurrence  formula 

p;  *  ^ 2  pV\  -  t1  <z2+r2>pVz  <2<> 

with  P*=l  . 
o 

The  field  associated  with  can  be  generated  utilizing  high  permeability 
pole-pieces  desiqned  to  follow  suitable  equipotential s  for  Eq.  23. 

If  the  thin  disk  is  located  at  the  plane  z=0,  and  centered  at  r=0,  the 
boundary  conditions  will  be  matched  (neglecting  fringing  at  the  rim)  and  the 
required  Z(r)  generated  inside  the  disk.  Notice  that  through  the  use  of  super¬ 
position,  ajrvy  field  of  the  form  H  (r)  =  E  a  r3n  can  be  synthesized.^ 

2  n=0  n 

Magnetostatic  Modes  and  Waves  in  Films  with  Nonuniform  In-Plane  Bias 

The  analysis  of  three-dimensional  modes  in  films  with  nonuniform  in-plane 

bias  by  means  of  a  coupled  integral  equation  approach  will  be  presented  at  the 

1 5 

1980  Conference  on  Maqnetism  and  Magnetic  Materials. 

The  prospect  of  guiding  magnetostatic  waves  is  of  considerable  interest 
because  of  possible  device  applications.^  Such  guided  waves  might  be  used  to 
increase  the  delay  time  realizable  on  a  given  size  sample  by  meadering  the  path. 
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or  to  make  a  resonator  by  guiding  the  waves  alonq  a  closed  loop.  In  addition, 
controlling  the  coupling  between  adjacent  waveguides  could  make  possible  signal 
routing  devices  such  as  directional  couplers. 

Guiding  such  a  wave  is  complicated  by  the  fact  that  MSSW  propagation  is 
only  possible  in  one  direction  on  a  given  surface  when  the  applied  in-nlane 
bias  field  is  uniform;^  turns  of  90°  would  normally  require  conversion  to  back¬ 
ward  volume  waves. 

It  should  be  possible  to  overcome  this  difficulty  by  employing  gradients 
that  arise  from  a  change  in  the  direction  of  the  bias  field.  As  an  example, 
consider  a  YIG  film  that  is  covered  with  permalloy  containing  a  slot  of  controlled 
width.  If  the  permalloy  is  at  a  different  magnetostatic  potential  on  either 
side  of  the  slot,  the  dc  magnetic  field  will  be  parallel  to  the  film  in  the 
region  underneath  the  slot  but  normal  to  the  film  surface  in  those  portions 
directly  underneath  the  permalloy. 

The  in-plane  fields  permit  magnetostatic  surface  wave  propagation;  the 
normal  fields  do  not  .  Therefore  the  surface  wave  energy  should  be  localized 
under  the  air-filled  slot  rather  than  under  the  conducting  permalloy;  oddv 
current  dissipation  is  thereby  minimized. 

Notice  also,  that  if  the  entire  bias  field  is  normally  directed,  the  surface 
wave  disappears  completely.  Taken  together,  these  factors  should  allow  novel 
control  of  the  surface  wave  channel . 

We  here  review  Stancil's  experimental  observations^  of  MSSW  propagation  in 
a  rectangular  YIG  film  centered  in  the  slot  between  strips  of  hiah  permeability 
metal  foil.  The  results  indicate  that  surface  waves  can  propagate  in  such  a 
geometry  but  with  modified  dispersion  characteristics,  as  expected. 

The  aecmetry  of  the  experiment  to  be  described  is  shown  in  Fig.  (7).  The 

sample  is  a  4.5  micron  thick  film  of  YIG  approximately  .28  cm  wide  and  1.09  cm 

long.  The  fine  wire  antennas  are  50  microns  wide  and  separated  by  1  cm. 


i 
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The  results  of  the  experiments  can  be  summarized  as  follows: 

1.  As  the  qradient  is  increased,  by  decreasinq  the  snace  d  between  the 
permalloy  strips,  a  series  of  discrete  modes  appears  on  the  low  frequency  side 
of  the  MSSVI  band.  This  behavior  is  illustrated  in  Figure  8. 

2.  Reversing  the  polarity  of  the  bias  field  results  in  si  ghtly  reducing 
coupling  to  both  the  continuous  and  discrete  modes  implying  that  these  modes 
have  a  nonreciprocal  surface  localization  similar  to  conventional  MSSVI's. 

3.  Although  the  discrete  modes  show  evidence  of  nonreciprocal  surface 
localization,  they  appear  at  frequencies  below  the  bottom  of  the  surface  wave 
band  at  the  center  of  the  film  as  calculated  from  the  corrected  fields  at  y=x=0. 
The  arrows  shown  in  Fiq.  8  indicate  the  bottom  of  the  MSSU  band  for  each  d  as 
calculated  in  this  way. 

4.  The  slope  of  the  nhase  was  used  to  calculate  approximate  group  velocities. 
The  discrete  modes  propagate  several  times  faster  than  the  normal  MSSU  modes. 
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Fiqurc  Captions 

figure  1.  Experimental  configuration  (a)  partially  exploded  view 

(b)  cross-sectional  view 

Figure  2.  Microwave  couplinn  between  two  antennas  as  a  function  of  frequency. 

Antennas  located  along  the  crystal  center  line  with  =  2300  Oe. 

Figure  3.  Coupling  between  antennas  as  a  function  of  crystal  position  with 
frequency  =  2.46  GHz  and  =  2300  Oe. 

Figure  4.  Locus  of  field  lines  for  the  "virtual  surface"  mode  n(m=l )  -  .655 

when  Z  =  .3+1.095  (r/R)^.  The  rf  m-field  is  shown  in  (a)  the 
la-field  in  (b)  and  the  b - field  in  (c).  All  patterns  rotate  at  the 
normalized  frequency  $2.  The  "virtual-surface"  (shown  dotted)  occurs 
at  rx/R  =  .385  . 

Figure  5.  The  normalized  mode  frequency  (a)  and  pre-speci fied  energy  frequency 
(b)  both  plotted  vs.  mode  number  toqether  with  the  required  field 
profile  (c)  plotted  vs.  radius. 

Figure  6.  The  normalized  mode  frequency  (a)  and  pre-specified  energy  frequency 
(b)  both  plotted  vs.  mode  number  toqether  with  the  required  field 
profile  (c)  plotted  vs.  radius. 

Figure  7.  Basic  experimental  configuration. 

Figure  8.  MSSW  transmission  spectra  as  a  function  of  the  spacing  between  the 
permalloy  strips.  The  fine  structure  near  the  high  frequency  end 
is  due  to  interference  with  the  EM  feedthrough. 
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Figure  8 


Effects  of  Exchange  in  Nonuniform  Fields 


I.  Introduction 


In  this  Chapter  we  will 

discuss  how  nonuniform  fields  along  with  the  exchange  interaction  can 
localize  volume  wave  excitations  in  a  normally  magnetized  thin  film  between 
conducting  plates.  Although  the  presence  of  the  conducting  plates  is 
undesirable  from  an  experimental  point  of  view,  a  significant  mathematical 
simplification  results  which  in  some  cases  allows  closed  form  solutions  to 
be  obtained.  These  solutions  give  insight  into  the  role  of  exchange  in  the 
localization  of  magnetostatic  waves  with  nonuniform  fields. 

The  simplest  field  nonuniformity  is,  of  course,  that  of  a 
one-dimensional  linear  profile.  Such  a  profile  has  been  used  to  synthesize 
linear  delay-vs. -frequency  characteristics  in  single  crystal  YIG  rods  (29). 
In  this  Chapter  we  discuss  the  possibility  of  a  similar  application  of 
linear  profiles  in  thin  ferrite  films.  Also,  since  an  arbitrary 
one-dimensional  profile  can  always  be  linearized  in  a  small  region  about  a 
point  of  interest,  the  understanding  gained  from  the  study  of  a  linear 
profile  has  very  wide  applications. 

A  profile  which  could  be  used  in  signal  routing  applications  is  the 
one-dimensional  quadratic  profile.  We  will  show  that  a  magnetostatic  wave 
can  be  confined  by  such  a  profile  and  made  to  propagate  along  the  direction 


transverse  to  the  gradient,  much  as  water  is  guided  by  a  trough. 

Third,  we  will  show  that  a  magnetostatic  mode  can  be  totally  confined 
by  a  two-dimensional  quadratic  field  profile.  Such  a  confined  mode,  or 
resonance,  could  form  the  basis  of  a  tunable  microwave  filter. 

As  a  bonus,  the  calculations  predict  the  existence  of  a  continuum  of 
spin  wave  modes  at  frequencies  above  the  top  of  the  volume  wave  manifold 
confirming  an  earlier  theory  described  by  Eshbach  (30)  based  on  a 
quasi-uniform  plane  wave  approach. 

Although  we  have  motivated  the  study  of  these  particular  profiles  using 
possible  device  applications,  we  should  point  out  that  the  primary  value  in 
the  solutions  to  be  described  is  in  understanding  the  role  of  exchange  in 
nonuniform  fields;  it  is  doubtful  that  modes  dominated  by  exchange  will  be 
of  great  practical  interest  due  to  hign  propagation  losses.  It  seems 
certain,  however,  that  analogous  mode  confinement  will  result  if  the  wave 
dispersion  is  dominated  instead  by  the  dipolar  interaction.  This  would  be 
the  case  in  the  more  practical  geometry  of  a  ferrite  slab  without 
conducting  plates. 

Although  the  present  study  was  begun  because  of  questions  raised  by  the 
virtual  surface  theory  (12,15),  our  geometry  differs  from  thac  of 
Morgenthaler ' s  in  one  important  respect.  The  properties  of  the  virtual 
surface  modes  depend  crucially  on  the  boundary  conditions  at  the  edges  of  a 
finite  width  or  diameter  thin  film.  In  our  geometry  the  film  is  of 
infinite  extent  in  the  plane;  hence  edge  effects  are  completely  neglected. 
It  is  clear  that  the  modes  obtained  by  the  two  methods  will  not  in  general 


be  in  one-to-one  correspondence.  However,  it  will  in  some  cases  be 
possible  to  draw  connections  between  the  two  theories,  and  these  will  be 
stressed  where  appropriate. 

II.  Expansion  in  Normal  Modes  of  the  Magnetization 

It  is  well  known  that  the  linearized  torque  equation  reduces  to  a 
Schrodinger  equation  when  dipolar  interactions  are  neglected  (31).  When 
solutions  to  this  are  known,  it  is  sometimes  possible  to  construct  a  set  of 
basis  vector-functions  from  which  solutions  to  Maxwell’s  equations  in  the 
magnetostatic  limit  can  be  constructed.  We  construct  a  set  of  basis 
vector-functions  in  a  manner  similar  to  that  used  by  Vendik,  et  al. 

(17,18). 

The  geometry  considered  is  that  of  a  thin  sheet  of  ferrite  between 
perfectly  conducting  plates  (Figure  2.1).  The  ferrite  is  assumed  to  be 
unbounded  in  the  x  and  y  directions.  The  material  is  magnetized  by  an  H 
field  whose  only  component  is  in  the  z  direction  but  which  has  some 
arbitrary  variation  along  x  and  y.  This  is  strictly  not  allowed  by  the 
requirement  7  x  H  *  0,  but  can  be  a  reasonable  approximation  if  the 
thickness  is  sufficiently  small.  Although  the  applied  field  is  nonuniform, 
the  material  is  assumed  to  be  saturated  everywhere.  The  presence  of  the 
conducting  plates  allows  us  to  consider  fields  which  depend  only  on  x  and 
y.  All  fields  depend  on  time  through  the  factor  exp(-iwt). 

As  shown  in  Appendix  A,  the  linearized  equation  of  motion  for  the 
magnetization  can  be  written 
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where 


'op 


z  “  *ex7t2 

-m 


Z  -  XexVt2 


v  2  »  a2  +  32 

"  3^2  3^2 


n  =  M  ,  2  -  H0eff/  Mg  , 

— YUoMg 

h  is  the  rf  magnetic  field,  H0e^  is  the  total  effective  dc  magnetic  field, 
Ms  is  the  saturation  magnetization,  Aex  is  a  phenomenological  exchange 
constant,  and  y  is  the  gyromagnetic  ratio  (negative). 

m 

Suppose  that  the  eigenvectors  and-  eigenvalues  of  Aop  exist  and  are 
known;  i.e.,  assume  the  equation 

Aop’S*  -  (2.2) 


has  been  solved.  Since  Aop  is  Hermitian,  the  eigenvalues  (a^)  will  be 
real. 

Let  us  assume  that  the  magnetization  can  be  expanded  in  terms  of  these 
eigenvectors  (the  success — or  failure — of  the  method  will  determine  the 
validity  of  this  assumption).  The  magnetization  can  then  be  written 

m  *  [  .  (2.3) 

l 


The  magnetic  field  can  now  be  obtained  easily  from  m  in  this 
representation: 

h  «  £  cjajSj  .  (2.4) 

l 

We  can  now  substitute  (2.4)  into  Maxwell's  equations  in  the 
magnetostatic  limit  and  look  for  the  coefficients  (c$)  and  eigenvalues  (ajj) 
such  that  the  equations  are  satisfied. 

Having  outlined  our  approach,  let  us  return  to  the  eigenvalue  problem 

(2.2). 

at 

Note  that  Aop  can  be  written  as  the  sum  of  a  scalar  operator  and  a 
matrix  operator 

-  -  f  0  i  ' 

Aoo  *  U  -  Xex7t2)I  +  a  (2.5a) 

L-i  0  . 

Aop  “  ^op  1  "*■  Pop  (2.5b) 

ac 

where  Rop  acts  on  the  space  of  arbitrary  functions  of  x  and  y  and  Pop  acts 
on  the  space  of  two-dimensional  complex  vectors.  As  a  trial  form  for  the 
eigenvectors,  let  us  consider 

S i  -  f(x,y)p  (2.6) 

where  f(x,y)  is  an  eigenfunction  in  the  space  of  Rop  given  by 
R0p  *  «  f(x,y) 

and  p  is  an  eigenvector  in  the  space  of  Pop  given  by 


(2.7) 
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(2.8) 


(The  reason  for  the  minus  sign  will  be  discussed  shortly.)  The  eigenvalue 
equation  (2.2)  becomes 

A0p  *  Sg  *  (Rgp  ^  ^  P0p)f(x,y)  P 

■  P  R0p  f(**y)  +  &  f(x,y)  Pop  ‘  P 

-  (p  -  Q  a)  f (x,y )  p 

-  (p  -  Q  o)  St  (2.9) 

which  implies  ag  *  p  -  fl  a.  The  problem  therefore  reduces  to  solving  the 
two  eigenvalue  equations  (2.7)  and  (2.8). 

Equation  (2.8)  is  readily  solved  and  the  solutions  are 


Pa 


1 

/~ 2 


1 

ai 


;  a  -  ±  1. 


(2.10) 


We  have  been  assuming  chat  Che  subscript  l  on  Che  eigenvector  Sg  is  a 
general  counting  index  which  ranges  over  all  possible  eigenvectors.  At 
this  point  it  is  advantageous  to  modify  the  notation  slightly  and  to 
substitute  the  subscripts  n  and  a.  It  is  understood  that  n  ranges  over  the 
solutions  to  (2.7)  and  a  ranges  over  the  solutions  to  (2.8).  Hence 


1  cl  s£  *  L  cn a  sna 
l  n,  a 


(2.11) 


The  vectors  pCT  are  the  polarization  eigenvectors  of  the  magnetization. 
The  vector  p+i  corresponds  to  right-hand  circular  polarization  (thumb  in 
the  z  direction)  and  p_i  corresponds  to  left-hand  circular  polarization. 
The  choice  of  the  minus  sign  in  equation  (2.8)  establishes  the  convention 
that  a  *  +1  corresponds  to  resonant  precession  (RCP)  of  the  magnetization. 

We  now  turn  our  consideration  to  (2.7).  We  have 

(Z  —  ^t‘)  fn(x,y)  ■  pn  fri(x,y) 

which  can  be  written 

Vc2  fn(x,y)  +  1  [pn  -  Z(x,y ) ]  fn(x,y)  -  0  (2.12) 

*ex 

Solutions  to  this  Schrodinger  equation  will  be  discussed  in  the  following 
sections. 

The  eigenvector-functions  CSna)  discussed  in  this  Section  can  be 
thought  of  as  the  normal  modes  of  the  magnetization  or  spin  wave  modes 
(17). 

Mote  that  if  Z(x,y)  were  a  constant,  (2.12)  would  be  a  simple  wave 
equation.  If  the  dimensions  of  a  finite  sample  were  taken  into  account, 
only  waves  with  certain  discrete  wavevectors  would  be  allowed.  This  was, 
in  fact,  the  case  in  the  problems  discussed  by  Vendik  and  Chartorizhskii , 
and  would  also  be  the  case  if  this  method  were  used  to  describe  the  Walker 
modes  (69)  of  a  uniformly  magnetized  sphere.  If  the  sample  were  infinite 
but  the  field  were  nonuniforra,  the  wavevectors  could  still  be  required  to 
be  discrete  just  as  the  energy  levels  of  a  quantum  mechanical  harmonic 
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oscillator  are  discrete.  In  this  manner  we  shall  see  that  the  gradient 
can  confine  a  mode  even  in  the  absence  of  boundaries. 

III.  Magnetostatic  Waves  Guided  by  One-dimensional  Bias  Field  Profiles 
A.  General  Formalism 

In  this  section  we  will  assume  that  the  bias  field  is  given  by  a 
function  of  x  only  and  that  fn(x,y)  is  of  the  form 

fn(x,y)  -  <frn(x)  e1*^  (2.13) 

Equation  (2.12)  becomes 

l£<t>n(x)  +_i_[Pn  -  lex  S2  -  Z(x)l  ^>n(x)  -  0  (2.14) 

3x2  Aex 

For  the  moment  we  will  assume  this  equaion  has  been  solved  and  consider  the 
constraints  imposed  by  Maxwell’s  equations  when  taken  in  the  magnetostatic 
limit.  These  equations  are  given  by 

7txh-0  (2.15a) 

7t  •  (Z  +  h)  -  0.  (2.15b) 

The  magnetization  is  given  by 

m  *  [  cno  $n(x)  ei6y  Po  (2.16) 

n,  a 

from  which  h  is  easily  obtained: 

h  -  [  ana  cn0  <frn(x)  eiSy  p0 
n,  a 


(2.17) 
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The  x  and  y  components  of  m  and  h  can  be  obtained  with  the  aid  of 
equation  (2.10): 


1 


®x  *  ~  l  cno  4>n<*)  eiSy 
v2  n,o 


(2, 


“y  “  ~  I  a  cnc  M*)  «iBy 
/  2  n,  a 


(2, 


hx  ”  _  [  ano  cno  'f'n C x ^  e*By 

Jl  n,  a 


(2. 


hy  *  ~  l  °  ana  cna  <)>n(x)  e^y 
fl  n,o 


(2. 


The  curl  equation  (2.15a)  requires 

illy  -  ill*  *  o 

3x  3y 


(2. 


which  becomes 


I  lano  cno  (o  iln  “  S<pn) ]  -  0 
n,  a  3x 


(2. 


In  a  similar  way  the  divergence  equation  (2.15b)  requires 


i_(a  +  h  )  +  i_(m  +  h  )  *  0 

3x  *  *  3y  y  y 


which  becomes 


(2. 


I  (d  +  ano>  cno  [ii"  "  8  a  >n  ] !  -  0 
n,  a  3x 


.18a) 


18b) 


19a) 


19b) 


20) 


21) 


22) 


(2.23) 
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B.  Solutions  for  8  *  0 

A  significant  simplification  occurs  for  the  limiting  case  8*0  (this 
is  the  cutoff  condition  for  propagation  in  the  y  direction).  Since  the 
derivatives  in  (2.21)  and  (2.23)  cannot  vanish  for  all  x  (except  in  the 
trivial  case),  we  have  the  conditions 

l  cno  a  -  0  (2.24a) 

n,  a 

l  0  +  *na)  cno  *  0  (2.24b) 

n,  a 

for  each  n.  Expanding  the  summation  in  (2.24a)  gives  the  relation 

cn+  „  an-  (2.25) 

cn-  an+ 

where  the  subscripts  t  are  shorthand  for  o  *  ±1.  Substituting  this  into 
(2.24b)  and  making  use  of  the  expression  for  ana  in  (2.9)  gives 

On2  "  ®n  (On  +  O  (2.26) 

where  pn  is  the  nch  eigenvalue  of  the  equation  (cf.  equation  (2.14)) 

+  _^_(pn  -  Z(x) )  -  0  (2.27) 

3x2  Xex 

Note  that  if  we  can  find  a  solution  to  (2.27)  for  a  given  Z(x)  we  have 
found  an  exact  solution  to  the  magnetostatic  wave  problem  for  the  case 
3  =  0.  In  such  a  case  the  frequency  is  given  by  (2.26),  the  polarization 
of  the  magnetization  by  (2.25)  and  h  and  m  are  determined  to  within  an 
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arbitrary  constant  by  (2.16)  and  (2.17)  evaluated  for  6  *  0. 

C.  Linear  Bias  Field  Profile 

Consider  Che  following  normalized  bias  field  profile: 

Z(x)  -  Z0  +  B  x.  (2.28) 

In  the  calculations  which  follow,  we  consider  only  the  case  8*0  for 
simplicity.  Equation  (2.27)  becomes 

32*n  +  1  (pn  -  Z0  -  B  x)  -  0  (2.29) 

3x2  *ex 

Let  us  make  the  following  change  of  variables 

x  -  a  £  -  (ZQ  +  pn)/  B  (2.30) 

where  a  is  a  characteristic  length  to  be  determined  and  5  is  a 
dimensionless  variable.  Upon  this  substitution  (2.29)  becomes 

324>n(0  -  B  a3  £  <t>n(£)  -  0  (2.31) 

3C2  Xex 

We  can  now  choose  che  characteristic  length  (a)  so  as  to  simplify  the 
equation.  We  therefore  choose 

a  -  Uex  /  B)l/3  (2.32) 

which  gives 

324>n(£)  .  C  *n<5)  -  0.  (2.33) 

The  solutions  to  this  equation  are  the  Airy  functions  A^(£)  and  B^(£). 

If  we  were  considering  a  general  boundary  value  problem,  a  linear 
combination  of  both  functions  would  be  necessary.  However,  for  the  present 
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infinite  medium  problem  only  Ai(C)  is  permitted1.  The  solution  in  terms  of 
x  is  given  by 


<j>(x)  *  A^x/a  +  (Z0  -  p)/Ba).  (2.34) 

The  subscript  n  is  no  longer  necessary  since  only  one  solution  is  being 
considered. 

According  to  tlorgenthaler  (15),  the  "virtual  surface"  in  the  absence  of 
exchange  would  be  located  at  the  point  of  inflection  (  *  0  of  the  Airy 
function  A^(5).  The  location  of  the  virtual  surface,  or  turning  point,  in 
the  x  coordinate  is  then  given  by  equation  (2.30): 

*tp  *  (P  -  Z0)/B  (2.35) 

where  the  subscript  "tp"  refers  to  the  turning  point.  The  frequency  is 
given  by 

SI2  -  (B  xcp  +  Z0)(B  xtp  +  Z0  +  1);  B  xtp  >  -Zo  (2.36) 

which,  when  xtp  ■  0,  reduces  to 

SI2  -  Z0  (Z0  +  l)  (2.37) 

which  is  the  same  frequency  obtained  by  tlorgenthaler  (15).  Although  the 

present  mode  does  not  have  precisely  the  same  functional  form  as  a  virtual 

surface  mode  corrected  to  Include  exchange,  we  see  that  both  theories 

predict  that  the  turning  point  occurs  at  the  position  where  the  field 

places  the  frequency  at  the  top  of  the  volume  wave  manifold,  which  is  where 

the  virtual  surface  occurs  in  the  absence  of  exchange.  (These  results  are 

lStrictly  speaking,  our  assumptions  of  a  linear  bias  field  profile,  an 
infinite  medium  and  saturation  of  the  magnetization  throughout  the  material 
are  not  self consistent ;  there  will  be  some  negative  x  for  which  the  field 
is  not  strong  enough  to  saturate  the  magnetization.  We  can  imagine, 
hovever,  that  the  material  is  terminated  with  a  perfectly  absorbing 
boundary  at  some  negative  x  at  which  the  magnetization  is  still  saturated. 


also  in  exact  agreement  with  those  of  Eshbach  (30)  even  though  he  only 
considered  dipolar  effects  in  an  approximate  way.) 

Note  that  in  the  present  case  p  can  assume  a  continuum  of  values  each 
corresponding  to  a  different  location  of  the  turning  point. 

This  solution  can  be  thought  of  as  a  standing  wave  resulting  from  a 
spin  wave  incident  from  the  left  and  being  totally  reflected  at  the  turning 
point.  Equation  (2.36)  indicates  that  the  higher  the  frequency  of  the 
wave,  the  further  it  can  penetrate  into  the  material.  This  can  be 
understood  qualitatively  through  a  consideration  of  the  strong  field  limit 
(B  xtp  +  Z0  »  1)  from  the  quasi-particle  viewpoint.  In  this  case, 
equation  (2.36)  states  chat  che  magnon  is  reflected  from  the  point  at  which 
the  strength  of  the  potential  energy  barrier  («  B  x  +  ZQ)  is  equal  to  the 
energy  of  the  magnon  (<*  ft). 

The  possibility  that  different  frequency  components  may  be  reflected  at 
different  depths  into  the  crystal  suggests  chat  a  linearly  dispersive  delay 
line  could  be  made  usign  this  principle.  A  major  problem  with  building 
such  a  device  would  be  coupling  to  the  very  high  wavenumber  spin  waves 
involved.  Based  upon  earlier  work  by  Strauss  (32)  and  others,  Platzker  and 
Morgenchaler  (29)  have  described  a  linearly  dispersive  delay  line  using  a 
linear  field  profile  in  a  YIG  rod.  Coupling  to  the  high  wavenumber  spin 
waves  in  their  case  was  accomplished  with  a  clever  scheme  making  use  of  a 
backward-wave  region  of  the  dispersion  diagram  (present  as  a  result  of  che 
boundary  conditions  involved)  and  phonon-magnon  coupling.  Perhaps  a 
similar  scheme  could  be  employed  to 
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fabrlcaCe  such  a  delay  line  using  chin  films. 

If  Che  problem  of  coupling  to  these  modes  can  be  solved,  their  presence 
might  be  directly  confirmed  using  optical  probing.  The  one-dimensional 
standing  wave  pattern  described  above  would  appear  as  alternating  stripes 
of  large  and  small  cone  angles.  This  could  act  as  a  diffraction  pattern 
for  properly  polarized  and  analyzed  light  passing  perpendicularly  through 
the  film  (33,34). 

We  have  mentioned  that  the  standing  wave  pattern  resulting  from  a 
linear  field  profile  would  involve  very  short  wavelength  (high  wavenumber) 
spin  waves.  It  is  instructive  to  consider  a  numerical  example. 

For  large  negative  arguments,  Che  function  Ai(-|?|)  is  approximately 

A1('l5l)  ’  /  2~  J  g  |  ~4  Sln  [  \  l5|3/2  +  11/4  1  (2'38) 

Because  of  the  Id322  in  the  argument,  the  period  of  the  sinusoid  will  vary 
with  posicion.  We  can  estimate  the  sinusoidal  period  near  a  position  in 
the  following  way: 

-(U0  +  ACl3/2  -  U013/2)  -  2tt 
3 

which  gives  AC  3  2ir//?^  .  But  C  *  x/a,  so 

Ax  =>  2n  j  q3  *  2*  /  ^ex  .  (2.39) 

/  Xq  /  Bx0 

where  x0  is  Che  distance  from  the  "virtual  surface."  For  a  gradient  such 
chat  B  »  1  cm-1  and  assuming  »  3  x  10“i2  cm2,  the  wavelength  1  cm  from 
the  virtual  surface  is  approximately  O.II  um. 


\ 

k 


D.  Quadratic  Bias  Field  Profile 


1. Basis  Eigenvector-Functions 

We  next  consider  a  normalized  field  profile  of  the  form* 

Z(x)  -  Z0  +  Bx2  (2. AO) 

We  begin  by  considering  the  problem  of  an  arbitrary  S.  In  this  case 
equation  (2.1A)  becomes 

a2«n(x)  +  J_(pn  -  Xex82  -  Z0  -  Bx2)  (j>n(x)  -  0  .  (2.A1) 

3^2  Xex 

As  in  the  linear  profile  case,  it  is  desirable  to  express  this  equation  in 
natural  units.  We  do  this  by  making  the  change  of  variable 

*  *  a0  5  (2.42) 

where  aQ  is  a  characteristic  length  to  be  determined  and  (  is  a 
dimensionless  parameter.  Equation  (2.41)  becomes 

lfV5)  +  r  (pn  -  Z0  -  Xexs2)  -  c2  ]  *„«•)  -  0.  (2.43) 

3C2  Xex  Xex 

It  is  now  clear  that  we  should  choose 

«o4  *  Xex  /  B  (2.44) 

which  gives 

IS  +  r  (pn  ~  zo  -  Xex82)  _  c2  1  -  o.  (2.45) 

m — * 

ex 

The  solutions  to  this  equation  can  be  immediately  obtained  by  noting 

chat  it  is  of  the  same  form  as  the  one-dimensional  harmonic  oscillator 

‘The  units  of  the  generic  gradient  parameter,  B,  vary  but  should  be 
apparent  from  the  context  (cf.  equations  (2.28)  and  (2.40)). 
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equation 


Z2*  +  (2c  -  52)  <l>  -  0 

w 


(2.46) 


Acceptable  solutions  to  this  equation  exist  only  for  discrete  eigenvalues 
which  are 


e  ■  n  +  1/2, 


n  ■  0,  1  .  .  . 


Identifying  terms  between  equations  (2.45)  and  (2.46)  gives 
pn  -  (2n  +  D/BX^  +  Z0  +  \exB2 


(2.47) 


(2.48) 


an0  -  (2n  +  l)/BXex  +  ZQ  +  Xex82  -  a,  o  ■  ±1. 


(2.49) 


The  solutions  to  (2.45)  are  the  Hermite-Gaussian  functions 


♦„(«  .  '*'2/2 

ftff'ni  /T)  1/2 


(2.50a) 


„  ,  ,  v  -x2/2a»z 

<0n(x)  .  Hn(x/a0)  e 

(a0  2n  n!  /“7)1/2 

where  Hn(£)  is  the  nc^  order  Herraite  polynomial  (see  Appendix  B). 
Combining  (2.6),  (2. 10) , (2. 13)  and  (2.50b)  gives 


(2.50b) 


sno  ■  “  eUy 

/  2  l  io 


(2.51) 


Using  the  orthogonality  properties  of  $n(x)  (see  Appendix  B)  it  is 
straightforward  to  show  that  these  vector-functions  are  orthonormal: 


(^n^Oi "  ‘  ^nl°l  *  ®O202^x  ”  ^Oia2  (2.52) 
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2.  General  Solutions 

Using  equation  (B.13)  for  the'  derivative  of  $n(x)  allows  us  to  write 
Maxwell's  equations  (2.21)  and  (2.23)  as 

CM 

l  l  {  g  an°  cn°  [$n-i  j  n  -  f  n+l]  -  ana  cno  S  d»n  }  .  o  (2.53) 
a  n«o  ®o  /  2  /  2 

l  l  ((I^nq)  cna  [  <l>n-i  /~tT  -  $n+1  j  tvK  ] 
a  n«o  ao  /  2  /  2 


8  a  (1  +  ana)  cna  d>n  }  m 


(2.54) 


Finding  the  coefficients  (cna)  and  the  frequencies  (fln)  which  satisfy  these 
two  equations  will  give  the  general  solutions  for  the  problem.  It  is 
useful,  however,  to  consider  approximate  solutions  constructed  with  only 
the  first  in  +  1  eigenfunctions.  The  summations  in  (2.53)  and  (2.54)  can 
then  be  reordered  to  give 


ra-1 

£  I  [  q  an-H,q  cn+l,q  j  n+1 
a  n-0  an  /  2 


-  a  a-- 


n-l,q  cn~l , o  fT 
aQ  /  2 


*  ®  ano  cno  ]  $n 


l  (  1 

o 


anr»l,j  cm-l»o  /n~  +  8  a, 
<>o  ✓  2 


na  '-m a 


]  *m 


\  CT  amq  cmo 
a 


m+J.  =  0 

2 


(2.55) 
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l  l  cn+l , a  /  n+1  ~  (1  1  an-l , cn-l,rr  /n~ 

a  n*0  a0  /  2  a0  /  2 

-  S  a  (1  +  ana)  cn0  ]  <frn 

-  £  [(1  +  ap-^j)  Cm.i>a  /m  +  6  a  (1  +  ama)  cma  ] 

a  a0  /  2 

-  I  ^  +  ama)  cma  /m+1  =0  (2.56) 

a  a0  /  2 

In  order  Co  satisfy  these  truncated  equations  exactly,  the  coefficient  of 
each  <(>n  must  vanish  because  of  orthogonality  giving  a  total  of  2m  +  4 
equations.  However,  there  are  only  2m  +  2  unknowns  so  that  in  general  the 
equadons  cannot  be  satisfied  exactly.  To  pursue  this  approach  further 
either  of  Che  following  procedures  can  be  attempted. 

a.  Minimize  the  space  average  of  the  equations. 

Instead  of  requiring  (2.55)  and  (2.56)  to  vanish  everywhere,  we  can  try 
to  minimize  the  average  of  their  absolute  squares  or  absolute  values  over 
all  space.  The  equations  so  obtained  would,  of  course,  be  nonlinear.  In 
addition,  although  we  constrain  S2n  to  be  real,  the  cna  coefficients  may  be 
complex  resulting  in  a  search  for  minima  in  a  4m  +  3  dimensional  space! 

b.  Satisfy  the  equations  exactly  at  a  finite  number  of  points. 

Alternatively,  we  can  require  the  equations  to  be  satisfied  at  2m  +  2 


points.  Although  this  method  results  in  a  set  of  simultaneous  linear 
equations,  it  is  not  clear  which  points  we  should  choose  in  order  to 
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optimize  the  solution. 

Because  of  these  difficulties,  we  will  not  consider  the  general 
solution  further.  Instead,  we  will  examine  these  equations  for  several 
limiting  cases  in  which  significant  simplifications  occur.  The  case 
for  B  ■  0  has  been  mentioned  already.  In  addition,  the  weak  gradient 
( ^BXex  <<  1)  and  strong  gradient  C*^BXex  >>  1  )  limits  will  be  examined. 

3.  Solutions  for  6*0 

Using  the  results  of  Section  III.B  we  can  immediately  write  down 
solutions  for  6*0.  There  is  a  solution  for  each  n  given  by 


mn  "  cn+  ^n+  I  +  cn-  Sn-  I  (2.57a) 

6*0  8-0 

hn  “  an+  cn+  sn+  I  _  +  an-  cn-  sn-  I  _  (2.57b) 

8*0  8*0 

where  the  polarization  of  the  magnetization  is 

cn+  -  .  (2.25) 

‘-n-  an+ 

an±  -  (2n  +  1)/  1\ex  +  ZQ  *  ^  ,  (2.58) 

and  the  frequency  of  the  nc^  mode  is  given  by 

:ln2  -  [(2n+l)/BAgX  +  Z0 ] [ ( 2n+l )/BXex  +  Z0  +  1 ]  (2.59) 


Note  in  this  case  that  although  the  magnetization  is  elliptically 
polarized,  the  small  signal  h  field  is  linearly  polarized  in  the  x 
direction  (this  can  be  seen  by  separating  (2.57)  into  x  and  y  components 
and  substituting  (2.25)  in  for  the  coefficients  cn0). 
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It  Is  worth  noting  that  localized  non-propagating  resonances  have  been 
observed  in  tangentially  magnetized  thin  films  (35,36).  In  such  cases  it 
has  been  proposed  (36)  that  the  inhomogeneity  in  the  saturation 
magnetization  near  the  interface  between  the  gadolinium  gallium  garnet  and 
yttrium  iron  garnet  could  give  rise  to  an  effective  "potential  well”  having 
one  bound  state.  Although  the  present  analysis  deals  with  an  inhomogeneous 
magnetic  field  rather  than  magnetization,  the  two  problems  appear  to  be 
very  closely  related. 

There  are  an  unlimited  number  of  modes  in  the  present  case,  however, 
since  we  have  assumed  an  infinite  paraboloidal  bias  field  profile.  It  is 
of  interest  to  examine  (2.59)  in  order  to  get  an  idea  of  the  frequency 
separation  of  these  modes.  Expanding  (2.59)  to  first  order  in  /  BXex  gives 

»  /Z0(Z0+1)  +  (2n+I )(2Z0+1 )^BXex  (2.60) 

2/  Z  (Z  +1)1 
o  o 

Using  the  parameter  values  (10)  of  Z0  ■  .78,  B  *  13  cm-2,  Xex  *  3xl0-12  cm2 
gives 

All  ■  Hn+1  -  fln  =  l.AxlO”5. 

This  is  at  least  an  order  of  magnitude  smaller  than  the  line-width  of  the 
best  materials  so  that  it  should  not  be  possible  to  experimentally  resolve 
individual  resonances.  Instead,  they  should  manifest  themselves  as  a 
quasi-continuous  absorption  or  transmission  band  immediately  above  Che 
volume  wave  manifold. 

It  is  instructive  to  find  the  locations  of  virtual  surfaces  for  these 


modes  as  we  did  in  the  case  of  the  linear  profile.  Using  Che  definition  of 


(2.61) 


a0  given  by  equation  (2.44),  the  mode  frequencies  can  be  written 

8n  -  [(2n+l)Bo02  +  Z0][(2n+l)Ba02  +  Z0  +  l] 

The  frequency  which  would  locate  a  virtual  surface  at  a  position  x  is  given 
by 

«2(x)  -  Z(x)[Z(x)+l] 

or  n2(x)  -  (Z0  +  Bx2]{Z0  +  Bx2  +  i]  (2.62) 

Given  the  frequency  of  the  nc^  mode  from  (2.61),  the  location  of  the 
virtual  surface  can  be  obtained  from  (2.62).  Clearly  the  expressions  are 
equal  when  (2n+l)a02  «  x2,  or 

£  -  x/a0  *  /2n+l  (2.63) 

But  this  is  simply  the  location  of  the  turning  point  of  the  nch 
Hermite-Gaussian  function  (see  Appendix  B).  Thus  we  find,  as  before  with 
the  linear  profile,  that  the  location  of  the  virtual  surface  is  coincident 
wich  the  turning  point  when  exchange  is  included. 

The  previous  analysis  of  propagation  in  a  linear  gradient  also  showed 
that  the  location  of  the  turning  point  could  be  varied  continuously  by 
changing  the  frequency.  The  higher  the  frequency,  the  further  into  the 
gradient  the  wave  could  penetrate.  Said  in  another  way,  Che  spectrum  for 
these  waves  is  not  a  single  resonance  but  a  continuum.  The  situation  in 
the  case  of  a  parabolic  profile  is  very  similar.  In  this  case,  however, 
there  are  two  turning  points  symmetrically  located  about  Che  center  of  the 
well.  The  additional  requirement  made  by  this  geometry  is  that  the 
osscillations  must  be  continuous  at  the  center.  This  requirement  splits 
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Che  continuum  into  a  series  of  resonances  which,  as  we  have  seen,  are 
very  closely  spaced  for  commonly  realized  gradients.  The  result  is  again  a 
quasi-continuous  mode  spectum.  We  shall  see  that  this  is  a  characteristic 
of  essentially  all  of  the  modes  accessible  with  the  present  formalism. 

Such  a  spin  wave  continuum  was  first  discussed  and  observed  by  Eshbach 
(30).  Eshbach  (30)  also  was  first  to  discuss  the  behavior  of  spin  waves  in 
linear  and  quadratic  field  profiles,  although  he  did  not  solve  for  analytic 
solutions  which  also  satisfy  Maxwell's  equations  in  the  magnetostatic  limit 
as  presented  here. 

Finally,  let  us  consider  the  spatial  localization  for  the  present 
example.  The  characteristic  length  is  given  by 
ao  “  (*ex/B)1/4  3  7 

Thus  we  find  that  the  modes  near  the  top  of  the  manifold  are  very  highly 
localized. 

4.  Strong  Gradient  Limit,  ^8Xex  »  1 

We  have  already  seen  that  a  relatively  small  gradient  can  give  rise  to 
a  highly  localized  resonance.  For  very  strong  gradients  we  can  Imagine 
that  the  localization  would  be  such  that  the  exchange  interaction  would 
completely  dominate  the  resonance  and  the  dipolar  interactions  could  be 
neglected.  In  this  case  equation  (2.4)  gives 

I  ana  cna  Bn a  3  (2.64) 

n,  o 

If  all  of  the  cna's  vanished,  then  so  would  the  magnetization  and  we  would 


have  the  null  solution.  Since  we  have  assumed  /BAex  >>  1,  all  of  the 
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an0's  cannot  simultaneously  vanish— even  in  an  approximate  sense— but  they 
can  vanish  individually.  From  (2.49)  we  see  that  for  0n  >  0  only  the 
eigenvalues  ana  for  a  •  +1  (resonant  precession)  can  vanish.  We  conclude 
chat  each  of  the  vectors  Sn+  represents  a  solution  with  the  frequency 

On  -  (2n+l)/BXex  +  Z0  +  Aex02.  (2.65) 

The  magnetization  for  each  mode  is  right  hand  circularly  polarized. 

5.  Weak  Gradient  Limit,  /BAex  <<  i 

Because  of  the  small  size  of  Aex,  the  assumption  that  /BAex  <<1  is 
well  satisfied  in  most  experimental  situations.  In  this  approximation  the 
fact  that  adjacent  an0's  are  degenerate  introduces  a  redundancy  which 
reduces  che  2m  +  4  equations  of  Section  11I.D.2  to  an  effective  set  of 
2m  +  2  which  can  be  solved  exactly. 

Expanding  equations  (2.55)  and  (2.56)  for  an  arbitrary  m  quickly 
becomes  cumbersome,  so  let  us  first  consider  a  graphical  representation  of 
the  equations  which  can  quickly  give  us  some  insight  into  which  and  how 
many  eigenfunctions  are  needed  in  order  to  construct  a  solution.  A  more 
rigorous  treatment  of  this  method  is  given  in  Appendix  C. 

Consider  a  representation  of  the  curl  equation 

_  _^x  »  o  (2.20) 

ax  ay 

illustrated  in  Figure  2.2(a).  The  indices  of  the  Hermite-Caussian 
functions  we  would  like  to  include  in  a  trial  solution  are  listed  in  che 
far  Left  and  right  columns  as  terras  in  hx  and  hy,  respectively.  The 
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indices  of  terms  resulting  after  taking  the  appropriate  derivatives  are 
listed  in  the  inner  two  columns.  By  assumption,  taking  a  derivative  with 
respect  to  y  simply  gives  the  same  function  times  a  constant  ( i B ) ,  whereas 
taking  an  x  derivative  gives  two  terms  with  raised  and  lowered  indices 
according  to  equation  (B.13). 

Equation  (2.20)  says  that  the  algebraic  sum  of  all  of  the  terms 
represented  in  the  inner  two  columns  must  vanish.  Orthogonality  of  the 
functions  ♦n(x)  further  requires  that  the  sum  of  all  terms  with  a  given 
index  must  vanish.  Thus  the  oval  enclosing  the  terms  representing  n  ■  1 
indicates  that  the  difference  between  these  two  terms  must  vanish.  Note 
that  there  is  only  one  term  with  n  *  2  so  that  its  coefficient  must  vanish. 
This  requires  all  of  the  elements  enclosed  by  the  rectangular  box  to  vanish 
since  they  are  all  related  through  non-zero  constants.  This  leaves  only 
one  n  »  0  term  so  its  coefficient  must  also  identically  vanish.  A  similar 
diagram  can  be  constructed  for  the  V  •  b  •  0  equation  as  shown  in  Figure 
2.2(b). 

These  diagrams  indicate  that  a  solution  may  exist  when  hx  «  6  and 
hy  *  <p0 .  If,  on  the  other  hand,  all  of  the  terms  in  both  hx  and  hy  were 
required  to  vanish  (as  symbolized  by  enclosing  all  terms  with  rectangular 
boxes),  there  would  be  no  nontrivial  solution  for  the  particular  set  of 
<t>n's  Included. 

We  have  said  that  the  diagrams  of  Figure  2.2  may  represent  a  solution. 
Experience  indicates  that  solutions  will  exist  only  under  the  following 


conditions: 
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1.  The  eigenvalues  (an0)  of  all  included  eigenvector-functions  (Sn0) 
must  be  degenerate  in  n. 

2.  The  fields  in  the  curl  (or  divergence)  diagram  oust  not  identically 
vanish. 

In  the  present  case,  both  conditions  are  satisfied  since  /BXex  «1.  We 
conclude  chat  the  diagrams  of  Figure  2.2  do  Indeed  represent  a  solution. 

Before  actually  calculating  the  fields  for  this  solution,  it  is  useful 
to  consider  a  slightly  more  general  case.  Let  us  begin  by  considering  the 
nc^  eigenfunction  and  attempt  to  construct  a  solution  by  adding  adjacent 
eigenfunctions  one  by  one.  The  diagrams  of  Figure  2.3  result.  Clearly,  a 
solution  exists  if  three  consecutive  eigenfunctions  are  included. 

We  will  now  calculate  the  field  quantities  for  this  mode.  From  (2.18) 
and  (2.19)  the  coefficients  which  are  required  to  vanish  by  the  rectangular 
boxes  in  parts  (c)  and  (d)  of  Figure  2.3  give 


cn-l ,+ 

.  an- 1 , - 

(2.66a) 

cn-l ,- 

an-l ,+ 

^n+  m 
cn- 

_  an- 

(2.66b) 

cn+l,+ 
cn+l ,- 

m  an+l,- 
an+l ,+ 

(2.66c) 

cn-l ,+ 

_  _  (1  +  an_it_) 

(2.67a) 

cn~l ,- 

+  an-l,+> 

Cn+  m 

(1  +  an_) 

(2.67b) 

cn- 

(1  +  an+) 

cn+l ,+  ,  _  ( 1  +  an+l 

cn+l,-  (1  +  an+l>+) 


(2.67c) 
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®  cn+l,+  8(an+»an-)  +  1  /n"**^  cn+  g(an^.,an.)  * 

“o/  2 


(2.72b) 


where  g(an+,an-)  is  a  common  function  of  each  term  and  can  be  divided  out. 
We  obtain,  finally: 

cn-l ,+  m  1  /  n 

^n+  0ao  /  2  (2.73a) 

cn+l,+  -  -  1  /n+T  (2.73b) 

cn+  6a0  /  2 

It  is  easily  shown  that  under  the  condition  /  BAex  <<  1  equations  (2.73) 
also  satisfy  the  remaining  conditions  of  Figure  2.3(d). 

Having  obtained  the  frequency  and  the  coefficients  to  within  an 
arbitrary  constant  we  can  now  write  down  the  fields.  We  have 


mx(n)  -  Cn+  {/7[l+l2±U  .  +  3a  flfl-HstU 

28n0  1  1  an-  n-1  o  1  an.J  n 

°y(n)  -  {/7[l-J“£]*  +  B«0f2[i-p±]*a 

20ao  an-  n-1  0  an-  n 

(2.74b) 

hx(n)  .  fn+fn+  [/7  ^  -  /Th  *n+l  ]  (2.74c) 

(2. 74d) 

(2. 74e) 


hy(°^  «  icn4an+/2  $n 
e,(°)  m  cn+UJu0  f  |_an+l 

s/T  1  ^rl 


n 
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where  we  have  obtained  the  above  approximation  to  the  first  order  electric 
field  from  the  Maxwell  equation 


7  x  e  -  -  Ik  *  (2.75) 

at 

Substituting  n  *  0  in  the  above  equations  allows  us  to  obtain  the 
fields  for  the  previously  discussed  example: 

mx(o)  -  c°-+  {  8a  /7[l-io±4  "  ]  (2.76a) 

26a0  0  ao-‘  0  "  ao- 

■»y(o)  -lf2±  {  8a  *  h-ao+U  1  (2.76b) 

28a0  1  0  1  0  1  ao-  1  J 

hx^°^  -  -  c0'^30'f  <t> ,  (2.76c) 

Ba0 

hy<«)  -  ic0+a0+/2  *  (2. 76d) 

1  o 


(2. 76e) 


Note  that  the  net  electromagnetic  power  flow  in  the  y  direction  is  zero 
for  all  of  the  waves  described  by  (2.74)  because  of  the  orthogonality  of 
the  Hermite-Gaussian  functions: 


p  (em)  >  _  pe|  /<exiT*)*iy  dx  f 


1 

2 


CO 

/  e2hx  dx 
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«  /  ij»n(x)i(>n_i(x)  dx  +  /  $n(x)$n+i (x)  dx 

»■  0.  (2.77) 

Thus  the  power  flow  in  these  modes  is  carried  entirely  by  the  exchange 
channel.  The  resulting  group  velocity  can  be  estimated  by  expanding 
(2.70a)  to  lowest  order  in  AexS2  and  taking  the  derivative  with  respect  to 
B.  The  result  for  n  small  is 

22  o  +  l 

Vg  -  ~ (2.78) 

'20<zo  +  i> 

It  should  also  be  pointed  out  that  the  waves  described  by  equations 
(2.74)  are  not  all  mutually  orthogonal  with  respect  to  the  inner  product 

CO 

(Fi.Fj)  -  /  Ti*-Fj  dx 

where  is  any  field  quantity  of  the  ich  solution.  The  solution  for  n  is 
orthogonal  to  those  for  n-1  and  n+1 ,  but  it  is  not  orthogonal  to  the 
solutions  for  n+2  and  n-2.  Hence  the  modes  for  n«0  and  n»l  can  be  taken  to 
be  the  first  two  solutions  of  an  orthogonal  set,  but  solutions  for  n>3  must 
be  carefully  constructed  so  as  to  be  orthogonal  to  all  lower  solutions. 
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IV.  Magnetostatic  Resonances  Bound  to  Two-dimensional  Quadratic  Bias  Field 
Profiles 


A.  Basis  Eigenvector-functions 

We  would  now  like  to  consider  normalized  bias  field  profiles  of  the 
form 

Z(x,y)  *  Z0  +  Bxx2  +  Byy2.  (2.79) 

The  eigenvalue  equation  (2.12)  becomes 

+  [a  +  Wo  -  Z  -  B  x2  -  B  y2]  *  -  0.  (2.80) 

W  W  0  x  * 


which  is  of  the  same  form  as  the  two-dimensional  harmonic  oscillator 
equation.  Following  the  usual  procedure  of  separation  of  variables,  the 
solutions  can  be  easily  shown  to  be  products  of  one-dimensional 
Hermite-Caussian  functions  (Appendix  B).  The  eigenvalues  are  calculated  to 
be 


an  n„o  "  ,/BxAex  <2nx+l)  +  /ByXex  (2ny+l)  +  Z 0  -  flo  . 


x  y 


(2.81) 


and  the  eigenvector-functions  are  given  by 


sn  n  a  * 
x  y 


$n  *n  <y> 
x  y 


(2.82) 


where  bn  (x)  and  bn  (y)  are  properly  normalized  one-dimensional  harmonic 
x  y 

oscillator  wavef unctions  as  defined  by  equation  (2.50b).  These 
eigenvector-functions  are  orthonormal  in  the  sense  that 


-114- 


/  /  dxdy  Sn  „  „*•  S„ 


n  a 
x  y 


_  m  a,m 
x  y 


,  m 

x  x 


Jn  ,m 

y  y 


o' 


(2.83) 


B.  General  Solutions 

We  assume  that  the  eigenvector-functions  (2.82)  form  a  complete  set  so 
that  m  and  h  can  be  expanded  in  terms  of  them.  We  can  then  write 


I  cn  n  o  ^n  n  o  (2.84a) 

nx,ny  *  y  X  y 

a 


^  "  I  annacnno^nna  (2.84b) 

nx,ny  x  y  *  y  x  y 

a 


The  individual  components  of  these  fields  are  found  to  be 


mx  “  -n  l  c n  n  o  *n  *n  (2.85a) 

/  2  nx,ny  ^  y  x  y 

a 


i  r 

®y  ■  1  0  cn  n  a  «n  <t>n  (2.85b) 

/  2  nx,ny  x  y  x  y 

o 

1  _ 

hx  "  1  an  n  a  cn  n  a  $n  ♦n  (2.86a) 

/  2  nx,ny  *  y  *  y  x  y 

a 


i 

7l  n 


r 

t 

x>n\ 


a  a 


n  n  a 
x  y 


-n  n  a 
x  y 


?n  vn 
x  y 


(2.86b) 


Applying  the  curl  equation  (2.20)  to  (2.86)  and  using  equation  (B.13) 


for  the  derivatives  gives 


I  ' y  {if.  a  c  [$  I  /**  ”  4  ,  /nx+i ] 4 

n  n  o  o-  nxny°  nxaya  nx-1/—  nx+1  / -  ny 

x’  y  /  2  /  2 

-  _L-a  c  *  [4  .  /v  -  4  /ny+1  ] }  -  0  (2.87) 

^  nxnya  nxnya  nxL  ny-l  J  _  ny+l  J  ... _ J 1 


Similarly,  the  divergence  equation  (2.22)  becomes 


I  I  { 


(  *+anynva) 


n  ,n  a 

*  y  a. 


Cnxnya[4>nx-1y^i  ~*nx+1  ^/nx+n 


+  ia(^+anvnva^ 


nynyg^nyf^ny-l ^*y  ~  ^ny+l  ^ny+l ] }-  0  (2.88) 


If  we  wish  to  keep  only  a  finite  number  of  terms  in  these  equations  we 
must  specify  precisely  how  we  are  going  to  count  terms.  In  order  to  do 


this,  let  us  define  a  principle  eigenfunction  index  n  such  that 
n  *  nx  ay  * 


(2.89) 


We  will  truncate  the  series  by  keeping  all  of  the  terms  (nx,ny)  such  that 
nx  +  ny  n.  Equations  (2.87)  and  (2.88)  can  then  be  written 


n  n-nx 

!  I  I  {  if  /f*  an  n  acn  n  a  V n  -i,n 

x-o  ny*o  o  1  ®x  /  2  x  y  x  y  x  > 


if  /f*il  an  n  ocn  n  a  *^n  +l,n 

0^/2  *  y  ^  y  x  y 


~  _i  /_Z  an  n  ocn  n  a  ’I'n  ,n  -I 

av  /  2  xy  xy  xy 


+  _i  /J-}.  an  n  ocn  n  a  4n  ,  n  +1 
av  /  2  x  y  x  x  y 


(2.90) 


-116- 


C.  SoluClons  for  Bx  ”  By 

In  order  to  avoid  the  complications  of  expanding  equations  (2.90)  and 
(2.91)  without  a  good  idea  of  how  many  and  whi.'h  terms  to  include,  we  will 
first  explore  different  combinations  graphically  using  a  straightforward 
extension  of  the  method  discussed  in  Section  III. 0.5.  The  rules  indicating 
the  existence  of  solutions  generalized  to  the  present  case  of  a  finite 
number  of  two-dimensional  Herraite-Gaussian  functions  are: 


I 
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1.  The  eigenvalues  (an  n  a)  of  all  included  eigenvector-functions 

x  y 

(Sn  n  a)  must  be  degenerate  in  (nx,ny). 
x  y  7 

2.  The  fields  in  the  curl  (or  divergence)  diagram  must  not  identically 
vanish. 

The  first  condition  is  met  for  all  neighboring  eigenvector-functions  in  the 
weak  gradient  limit  (/  Bi*ex  «1,  i“x,y).  It  is  also  met  for  all 
eigenvector-functions  belonging  to  the  same  principle  eigenfunction  index 
when  Bx  »  By.  This  latter  case  is  the  one  we  will  examine  in  this  section. 

The  curl  equation  diagrams  for  the  first  four  principle  eigenfunction 
indices  are  shown  in  Figure  2.4.  It  appears  that  no  solutions  are  possible 
for  n  even,  but  solutions  do  exist  for  n  odd. 

The  fields  can  be  obtained  from  the  diagrams  in  a  manner  very  similar 
to  the  one-dimensional  case  discussed  earlier.  We  first  make  the  formal 
identifications : 

if  index  is  decreased  from  nu, 

(2.92) 

if  index  is  increased  from  nv. 

We  must  also  remember  that  the  y  component  of  each  field  quantity  has  a 
factor  of  i  (this  factor  was  cancelled  by  the  i  in  3/3B  *  iS  in  the 
previous  case). 

Let  us  first  consider  the  solution  for  n-1.  The  components  of  hx  and 
hy  which  are  required  to  vanish  by  Figure  2.4(b)  imply 
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hx 

3/3y 

3/3x 

hy 

00 

♦ 

01 

10 

00 

(a) 

Fields 

vanish  for  n“0. 

01 

00,02 

01 

10 

* 

00,20 

10 

(b)  A 

solution  exists  for  n“l. 

02 

01,03 

1  1 

12 

02 

1! 

10,12 

01,21 

li 

20 

-► 

21 

10,30 

20 

(c)  Fields  vanish  for  n«2. 


(d)  A  solution  exists  for  n«3. 


Figure  2.4.  Curl  Diagrams. 
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col+  m 

_  *1- 

(2.93a) 

1 

1 

Col“ 

al+ 

i 

\ 

clo+  . 

a  1- 

(2.93b) 

clo- 

al+ 

Referring  to  Figure 

2.5,  the  components  of 

bx  and  by  which  are  required  to 

vanish  give 

col+  m 

( 1+a  i_) 

(2.94a) 

col- 

(l+ai+) 

clo+  . 

-  (1+ai-) 

(2.94b) 

clo- 

(l+ai+) 

Equations  (2.93)  and  (2.94)  can  be  satisfied  only  if 

nl2  “  [4/  B\ex  +  Z0][4/  BXex  +  Z0  +  1 ]  (2.95) 

where  we  have  used 

ano  -  2/  BAex  (n+1)  +  ZQ  -  tto  (2.96) 

and  B  ■  Bx  »  By. 

Requiring  the  difference  between  the  terras  enclosed  by  the  oval  in 
either  the  curl  or  divergence  diagrams  to  vanish  gives  the  relationship 
between  the  (01)  and  (10)  eigenfunctions.  From  the  curl  equation  we  have 

-  _ ]_  ci0+  g(ai+,a1_)  +  _ i_  col+  g(a1+,a!_)  -  0 

a /  2  a/  2 

where  g(ai+,ai-)  is  a  function  common  to  both  terras.  This  equation 
simplifies  to 


clo+  ,  j. 
co  1+ 


(2.97) 
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(a)  n-1 


(b)  n»3 


Figure  2.5.  Divergence  Diagrams 
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We  now  have  all  of  Che  information  needed  eo  write  down  the  fields  for  the 


We  have 

mxO>  -  CoH  + 

rT  3i-  01 

i[ni!+]*lo| 

aj-  10 

(2.98a) 

(2.98b) 

hX(1^  ■  icol+  ^  al+  *10 

(2.98c) 

hy(1)  -  icol+  /2  a1+  ii>0i 

(2.98d) 

ez<l)  -  -ia)ua2c0i+(l+ai+)a 

’('oo 

(2.98e) 

Spatial  maps  of  the  field  strengths  and  polarizations  of  m  and  h  are 
presented  in  Figure  2.6.  Note  Chat  h  is  linearly  polarized. 

Repeating  this  procedure  for  n-3  (see  Figures  2.4(d)  and  2.5(b))  we 
obtain  the  following  fields: 


(3) 


Dy 


{[i-12t]*03  +4_[  »£2±]*12  +  __[i -12±]*21  +  i[i+±2± 


n 


a3_ 


rr  a3. 


12 


n  a3- 


a3- 


30  1 


(2.99a) 


(3) 


1£2+|[m«+]*03  ♦  ‘  [i-12+]*12  * _[i+12+]*21  *  i( i-i2+l*30 1 
rr  a 3-  n>  33-  n,  a3-  a3- 


(2.99b) 

ico3+/2  a3+[  *1_2  +  *  ] 

n  30 

(2.99c) 

ico3+/2  a 3+  [  i|>  *21  ] 

03  n, 

(2. 99d) 
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(c) 


Figure  2.6.  Fields  of  Che  n“l  solution  in  the  presence  of  a  z  directed 
bias  field  which  is  a  symmetri£  paraboloid  of  revolution.  (a)  Relative 
magnitude  and  polarization  of  n>.  The  ellipses  indicate  the  locus  of  points 
traced  by  the  tip  of  the  small  signal  magnetization  vector  as  ic_  rotates 
counterclockwise.  (b)  Relative  magnitude  and  polarization  of  h.  __ 

(c)  Relative  magnitude  of  the  z  directed  electric  field  (ez).  Unlike  ra, 
both  h  and  e  are  linearly  polarized. 
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e  (3)  .  -i(Du0c03+  (l+a3+)a  2  f  /  2  +  \J>  +  ]  (2.99e) 

oo  02  20 

and 

«32  -[  8/  BAex  +  Z0][  8/  8Aex  +  Z0  +  1  ]  .  (2.100) 

SpaClal  maps  of  Che  field  strengths  and  polarizations  for  this  solution  are 
shown  in  Figure  2.7. 

In  general,  a  solution  can  be  constructed  from  the  n+l 
eigenvector-functions  belonging  to  the  principle  eigenfunction  index  n  when 
n  is  odd.  The  nc^  solution  has  rotational  symmetry  about  the  z  axis  and 
the  magnetic  field  has  (n+l)/2  radial  nodes.  The  frequency  of  the  nch 
solution  is  given  by 

-  [  2/  BXex  (n+l)  +  Z0][  y  BAex  (n+l)  +  Z0  +  1 ]  (2.101) 


D.  Strong  Gradient  Limit,  /  BXex  >>l 

If  the  magnetic  field  gradient  is  so  large  that  the  exchange 
interaction  completely  dominates  the  behavior  of  the  magnetization,  the 
magnetization  can  be  approximated  by  any  of  the  eigenvector-functions  which 
can  have  a  zero  eigenvalue.  The  solutions  are 

®n  n  +  *  x  y  4>n  (xHn  (y)  (2.102) 

x  y  rr  L  i  x  y 
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Figure  2.7.  Fields  for  n»3  solution.  (a)  Magnitude  and  polarization  of  m, 
(b)  Magnitude  and  polarization  of  h,  (c)  Magnitude  of  ez. 
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Here,  as  in  the  one-dimensional  field  profile  of  Section  I1I.D.3,  the 
magnetization  of  each  mode  is  right  hand  circularly  polarized. 

It  is  interesting  to  note  that  when  Bx  and  By  differ  slightly,  the 
spectrum  of  resonances  given  by  equation  (2.103)  consists  of  a  series  of 
multiplets  qualitatively  very  similar  to  those  observed  by  Cooley  in  a  disk. 
(37).  However,  in  order  for  the  resonances  to  differ  by  only  a  few 
megahertz  for  YIG  (f,  «  5  GHz)  we  must  have  /BAex  *  10“3,  or  B  3  3(105) 
cm-2.  We  conclude  that  although  the  qualitative  behavior  of  equation 
(2.103)  is  very  suggestive,  it  is  extremely  unlikely  that  field  gradients 
of  this  magnitude  existed  in  the  experiment  reported  by  Cooley!  This 
disagreement  results  in  part  from  our  neglect  of  surface  magnetic  poles 
which  in  the  magnetostatic  modes  of  a  disk,  as  in  the  Walker  modes  of  a 
sphere  (69),  play  a  dominant  role  in  determining  the  mode  characteristics. 
It  is  therefore  possible  that  Cooley's  results  could  be  explained  by  a 
similar  analysis  if  the  effects  of  finite  sample  edges  were  included. 

E.  Weak  Gradient  Limit  /BAex  <<  1 

When  terms  of  order  /BAex  can  be  neglected,  combinations  of 
eigenvector-functions  from  adjacent  principle  eigenfunction  indices  may  be 
considered  as  a  trial  solution.  Clearly  the  "exact”  solutions  of  Section 
IV. C  are  also  valid  in  this  limit  when  Bx  ■  By.  In  addition,  similar 
solutions  can  be  found  when  Bx  *  By  if  che  fields  are  rederived  taking 
into  consideration  that  ax  *  oiy  (this  will  destroy  the  symmetry  and 
elongate  the  mode  patterns  in  che  direction  of  the  weakest  gradient). 

Since  we  are  no  longer  restricted  to  eigenvector-functions  belonging  to 
che  same  principle  eigenfunction  index,  the  diagrams  shown  in  Figure  2.8 
may  now  be  considered.  Notice  that  the  diagrams  shown  in  Figure  2.8(a)  and 
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2.8(b)  transform  into  one  another  under  the  interchange  of  x  and  y.  Thus 
the  two  solutions  are  equivalent  under  a  rotation  of  it/2  about  the  z 
axis1.  It  is  sufficient,  then,  to  consider  the  solution  represented  by 
Figure  2.8(a). 

Following  a  now  familiar  procedure,  the  fields  corresponding  to  Figure 


2.8(a)  are  found  to  be 

■*  ■  y?  U‘-^l*oo  '  ^  n  I'-flV 

(2.104a) 

my  -  iCoo+{[l+a+]  -  iaj,  [!-•+]♦  -  n  [1+*+]*  } 

/2  a-  Ox  a_  1 1  a_  u *■ 

(2.104b) 

hx  *  -iCoo+/I  a+ 

(2.104c) 

«x 

hy  -  ic00+  /2  a+(<i»oo  -  /I  *02> 

(2. 104d) 

ez  -  io>u0c00+2(  l+a+)ay  * 

(2.104e) 

and 

U2  -  Z0(Z0  +  1), 

(2.105) 

where  a0  3  ZQ  -  fla.  Spatial  maps  of  the  field  strengths  and  polarizations 
for  this  solution  are  shown  in  Figure  2.9. 

In  the  weak  gradient  limit,  all  solutions  near  a  given  n  have 
approximately  the  same  frequency  which  for  large  n  is  approximately 

On2  -  [  2n/  BAex  +  Z0][2n/  8Aex  +  Z0  +  1 ]  (2.106) 

Again,  we  find  that  the  frequencies  of  che  modes  are  so  closely  spaced  thac 
they  form  a  quasi-comtinuum. 

‘The  diagrams  do  not  explicitly  indicate  the  signs  of  the  quantities 
involved.  The  actual  transformation  is  x  *  y,  y  *  -x. 
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Flgure  2.9.  Fields  of  the  lowest  order  solution  which  does  not  possess 
rotational  symmetry  about  the  z  axis  even  when  the  bias  field  is  symmetric. 
The  format  is  similar  to  that  of  Figure  2.6.  (a)  Magnitude  and 

polarization  of  ra.  (b)  Magnitude  and  polarization  of  h.  (c)  Contour  map 
of  ez. 
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It  should  be  possible  to  construct  additional  solutions  in  such  a  way 
that  the  solution  of  this  Section  and  the  solutions  of  Section  IV. C  in  the 
Unit  /BXex  <<  l  are  three  of  the  lowest  order  solutions  of  a  complete 
orthogonal  set. 

V.  Discussion 

If  exchange  were  neglected,  the  dispersion  relation  for  k  J_  H  in  an 
infinite  medium  (or  in  a  normally  magnetized  slab  between  conducting 
plates)  is  simply  fl2  »  Z0(Z0+1),  independent  of  S.  Hence  in  the  virtual 
surface  theory  (1,2),  at  the  point  in  the  material  where  the  field  and 
frequency  have  this  relationship  all  wavelengths  are  possible. 
Qualitatively,  this  gives  an  idea  of  how  a  singularity  in  the  field  might 
be  possible  at  this  point. 

The  addition  of  exchange,  on  the  other  hand,  lifts  this  degeneracy  and 
results  in  a  one-to-one  correspondence  between  11  and  8,  thus  eliminating 
Che  possibility  of  a  singularity.  In  the  case  of  a  nonuniform  field,  the 
resulting  mode  is  oscillatory  in  regions  of  the  sample  where  propagation  is 
permitted,  and  decays  rapidly  outside  of  these  regions. 

This  suggests  a  mechanism  whereby  a  gradient  could  localize  a  mode  in  a 
Chin  slab  with  Che  conducting  plates  removed.  As  shown  in  Appendix  A,  che 
effect  of  Che  new  boundary  conditions  is  to  introduce  a  component  of  e  x  h 
power  flow  which,  in  most  cases,  dominates  the  power  flow  in  the  exchange 
channel.  Hence  one  can  imagine  a  mode  which  again  can  propagate  in  some 
regions  and  not  in  others,  but  whose  propagation  characteristics  are 
determined  completely  by  the  magnetostatic  boundary  conditions  rather  than 


exchange.  This  general  concept  of  mode  localization  should  not  be  limited 
to  volume  waves,  but  should  apply  to  other  bias  field  orientations  as  well. 
It  is,  in  fact,  the  conceptual  basis  for  the  surface  wave  theory  described 
in  Chapter  Three. 

Finally,  two  limitations  of  the  present  theory  deserve  mention: 

1.  In  many  cases  (including  that  of  Eshbach  (30))  the  magnetoelastic 
crossover  point  occurs  within  frequencies  of  interest  and  cannot  be 
neglected.  The  present  theory  is  strictly  valid  only  where  this  is  not  the 
case. 

2.  As  pointed  out  in  the  introduction,  these  solutions  lose  their 
validity  if  the  turning  point  occurs  near  the  edge  of  a  finite 
width/diameter  film.  In  such  a  case,  the  second  solution  to  the  harmonic 
oscillator  equation  (which  does  not  vanish  at  infinity)  must  be  included  in 
the  eigenfunction  expansion. 

VI.  Summary  of  Chapter 

We  have  extended  the  method  of  expansion  in  normal  spin  wave  modes 
(17,18)  to  the  case  of  a  two-dimensional  non-unif ormly  magnetized  ferrite 
with  the  effects  of  exchange  included.  Using  this  method,  a  number  of 
solutions  (both  exact  and  approximate)  to  the  linearized  magnetostatic 
equations  have  been  obtained. 

A  consideration  of  the  special  case  of  a  one-dimensional  linear  bias 
field  profile  and  propagation  parallel  to  the  gradient  showed  that  the 
profile  of  the  magnetization  is  given  by  an  Airy  function.  The  turning 
point  of  the  Airy  function  was  shown  to  correspond  to  the  location  of  a 


virtual  surface  in  qualitative  agreement  with  previous  work  by  tlorgenthaler 
(15).  In  addition,  it  was  shown  that  the  distance  penetrated  into  the 
gradient  is  roughly  proportional  to  the  frequency  of  the  spin  wave. 

Next,  spin  wave  propagation  down  a  "trough”  formed  by  a  one-dimensional 
parabolic  bias  field  was  considered.  It  was  shown  that  solutions  to  the 
magnetostatic  problem  can  be  constructed  from  basis  vector-functions 
consisting  of  the  product  of  a  polarization  vector  and  a  Hermite-Gaussian 
function.  In  the  weak  gradient  limit,  /BXex  «  1,  it  was  shown  that  the 
wave  can  be  locallized  to  a  very  narrow  “trough”  (as  narrow  as  15  um  for 
B  *  12  cm“2)  when  the  frequency  is  very  near  the  top  of  the  volume  wave 
manifold.  The  net  Poynting  flux  was  found  to  be  zero,  however,  indicating 
chat  power  can  flow  only  through  the  exchange  channel. 

Finally,  localized  resonances  confined  by  two-dimensional  parabolic 
bias  field  profiles  were  considered.  In  this  case,  basis  vector-functions 
consisting  of  Che  product  of  a  polarization  vector  and  a  two-dimensional 
Hermite-Gaussian  function  can  be  used  to  construct  solutions  to  the 
magnetostatic  problem.  A  class  of  exact  solutions  (good  for  any  strength 
quadratic  gradient)  was  found  for  the  symmetric  case  Bx  »  By. 

Boch  the  frequencies  and  field  distributions  of  the  solutions  are  in 
agreement  with  the  spin  wave  continuum  discussed  and  observed  by  Eshbach 
(30). 

These  solutions,  free  of  the  singularities  encountered  in  the  virtual 
surface  theory,  demonstrate  that  exchange  can  play  an  important  role  in  the 
localization  of  magnetostatic  waves  in  nonuniform  bias  fields. 
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Dispersion  Relations  for  tlagnetostacic  Waves  in  Thin  Ferrite  Films 

I .  Background 

Magnetostatic  waves  in  a  ferrite  slab  with  an  in-plane  magnetic  field 
were  first  discussed  by  Danon  and  Eshbach  (54)  in  1961,  and  waves  in  a 
normally  magnetized  slab  were  first  discussed  by  Damon  and  Van  de  Vaart 
(55)  in  1965.  since  these  initial  studies,  various  modifications  have  been 
made  to  the  basic  theories.  Some  of  the  most  notable  of  these  are: 

1.  Effects  of  a  nearby  parallel  ground  plane  (56-59) 

2.  Effects  of  crystalline  anisotropy  (51,60-64) 

3.  Effects  of  finite  slab  width  (38,44,65,66) 

In  the  following  sections  we  derive  the  general  susceptibility  tensor 
including  anisotropy  and  exchange  and  present  the  basic  ecuations  of 
magnetostatics.  We  then  derive  the  basic  dispersion  relations  for  the 
three  principle  normal  modes  of  an  isolated,  isotropic,  infinite  width 
slab.  Dispersion  relations  for  ferrite-dielectric-metal  and  finite  width 
geometries  are  also  given  but  without  derivation. 

1 1.  The  Gene ral  Susceptibil ity  Tensor 

The  equation  of  motion  for  the  magnetization  is  given  by 

3M  _  _  _  _ 

—  =  y  p0  M  x  (H  -t-  Ha  +  H„x) 

3 1 


(A.  1) 
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where  y  m  -  g|el/2m  Is  Che  gyronagnetic  ratio,  M  is  the  magnetization,  H  is 
the  Maxwellian  field,  and  Ha  and  Hex  are  effective  fields  due  to  anisotropy 
and  exchange,  respectively.  The  effective  fields  are  defined  by 


_  1  3W3 

Ha  -  -  -  —  (A.  2a) 

Uo  3M 


H 


ex 


'ex 


V2*? 


(A. 2b) 


where  Wa  is  the  anisotropy  energy  density  (Section  III)  and  Xex  is  a 
phenomenological  constant  equal  to  approximately  3(10-12)  cm2  for  YIG. 

Following  the  usual  linearization  procedure,  these  fields  are  divided 
into  static  and  time  varying  components  as  follows: 


M  -  fi0  +  m 


H  -  H0  +  h 


*  ^oa  +  ^a 


^ex  ”  ^oex  +  *^ex 


(A. 3a) 

(A. 3b) 


(A. 3c) 


(A. 3d) 


The  field  H0  should  understood  to  be  the  vector  sum  of  the  actual 
externally  applied  magnetic  field  and  the  internal  demagnetizing  field,  but 
we  will  not  explicitly  indicate  this  in  order  to  avoid  further  notational 
complexity,  flaking  these  substitutions  into  (A.l)  and  assuming  the  time 
variation  exp(-iujt)  gives 
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— iwn  _  _  _  _  _  _ 

“  x  [So  +  ^oa  +  ^oex  J  +  ^o  x  H'  +  +  ^ex  1 

™0 

+  m  x  [Ho  +  Ha  +  Hgx  ]  +  m  x  [h  +  ha  +  hex  ]  (A.  4) 

The  equilibrium  direction  of  the  magnetization  is  determined  by  the 
e  qua  tion 

M0  x  [H0^0>  +  Ha(M0)  +  Hoex(MQ)  ]  -  0  (A.5) 

where  we  have  explicitly  indicated  the  dependence  of  H0,  HQa,  and  HQex  on 
M0.  The  dependence  of  H0  on  M0  results  from  the  demagnetizing  field 
mentioned  previously.  Equation  (A.5)  represents  a  set  of  three  coupled, 
nonlinear,  simultanious  equations  which  are,  in  general,  quite  difficult  to 
solve.  They  have  a  rich  family  of  solutions,  however,  including  Bloch 
walls  and  bubble  domains. 

Ac  present  we  are  primarily  interested  in  single  domain  materials  in 
which  case  Hoex  vanishes  by  (A. 2b).  If  in  addition,  |H0|  >>  I Ha I  ,  we  can 
assume  M0  II  HQ.  This  approximation  is  valid  for  many  cases  of  interest, 
but  should  be  examined  closely  for  lower  microwave  frequencies  and  strong 
anisotropies. 

Making  the  strong  field  approxina Cion  and  neglecting  terms  second  order 
in  small  quantities  enables  us  to  write  (A. 4)  in  Che  form 
_  ___  * 

ifai  »  z  x  +  Na ■ n  +  h  -  (Z0  +  Za)m]  (A.b) 

where  we  have  introduced  several  new  quantities: 
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(ii 

(A. 7a) 

^0*0 

ha  -  Na  •  m 

(A. 7b) 

Zo  “  Ho  /  «o 

(A. 7c) 

Za  ”  Hoa  *  2  /  ^o 

(A. 7d) 

.  32  32 

7.2  >  —  +  — 

C  3x2  3y2 

(A. 7e) 

and  we  have  chosen  the  coordinate  system  so  that  the  equilibrium 
magnetization  points  along  the  z  direction.  The  form  of  the  small  signal 
anisotropy  field  is  justified  by  noting  that  this  field  will  only  be 
present  if  there  is  some  deviation  from  the  equilibrium  direction  of  Tt. 

A 

Keeping  only  the  z  component  of  Hoa  is  justified  since  the  other  components 
give  rise  to  an  mz  which  we  neglect  as  a  seconi  order  quantity. 

Solving  for  h  in  (A. 6)  gives 

h  -  Aop  •  m  (A. 3a) 


where 


'op 


7  +7  -  M  -l  7  2 

4o  'xx  Aex  t 


■ii)  -  N, 


in  -  n 


‘’xv 


a  - 

7+7-M  -lqZ 

uo  6a  v  Aex  t 


(A. 8b) 


This  operacor  (with  anisotropy  set  to  zero)  is  the  starting  point  for  the 
exchange  calculations  of  Chapter  Two. 


1 

ll 
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If  exchange  is  neglected,  equation  (A. 8)  can  be  easily  inverted  to 
obtain  the  Polder  susceptibility  tensor 


m“X‘h  (A. 9a) 

where 

Xxx  *  <zo  +  za  ~  Nyy ) /D  (A. 9b) 

Xxy  *  (~iQ  +  NXy)/D  (A. 9c) 

Xyx  -  (ifl  +  N^)/D  (A.9d) 

Xyy  -  (z0  +  za  -  Nxx)/D  (A.9e) 

D  -  (Z0  +  Za  -  Syy)(Z0  +  Za  -  Mxx) 

-  (iQ  +  +  4y)  (A.  9f ) 


Bajpai,  et  al.  (63)  have  derived  an  expression  analogous  to  (A. 9)  for 

a  a 

the  permeability  censor  (I  +  x)  for  an  arbitrary  orientation  of  the  dc  bias 
field.  Their  result  is  incorrect,  however,  since  their  tensor  is  not 
Hermitian  (this  is  a  requirement  since  their  calculation  does  not  include 
1 oss) . 


If  anisotropy  is  neglected,  the  susceptibily  reduces  to  the  familiar 

tom 

-  X  "it 

X  ■  (A. 9g) 

-  i<  X  ■ 


whe  re 


X 


Z 

Z2  -  * 


<  * 


■2 
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III.  The  Effective  Anisotropy  Field  (67) 

We  consider  here  the  anisotropy  In  a  cubic  ferromagnet.  Such  a  model 
is  appropriate  for  single  crystal  bulk  YIG,  but  inevitably  YIG  thin  films 
also  exhibit  a  stress-induced  uniaxial  anisotropy  due  to  the  slight  lattice 
mismatch  between  the  YIG  and  GGG  substrate.  This  mismatch  can  be 
minimized,  however,  by  doping  with  suitable  nonmagnetic  impurities  (68). 

For  a  cubic  ferromagnet,  the  anisotropy  energy  density  to  lowest  order 
is  given  by 

K 

Wa  -  -  [Mi2M22  +  M22M32]  (A. 10) 

where  the  subscripts  (1,2,3)  refer  to  the  principle  crystal  axes.  The 
effective  anisotropy  field  is  defined  by 

1  3W® 

«a  - = 

uo  an 

Ml  (M2  2  +  M32) 

2K 

- - -  M2(Mj  2  +  M32)  (A. 11) 

U0N0 

M3(Mj 2  +  M22) 

separating  the  static  and  time  varying  components  and  keeping  terns  to 
first  order  in  m  gives 

Mqi(Mq22  +  M032) 

2K  =  _ 

Ha  - - -  M0t(M01  2  +  M() 3 2 )  +  Na  •  m  (A. 12a) 

Uo-Mo 

M03(M0i 2  +  m022) 
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where 


a  2K  ,  , 

Hi - <Mo2  -  Hoi  > 

UoM0 


4K 

^o 


’ij  “  ~  M°i  M°j  ;  L  *  j 


(A. 12b) 


(A. 12c) 


Again,  these  quantities  are  expressed  in  terms  of  the  major  cubic  axes  of 
the  crystal.  In  order  to  obtain  the  corresponding  quantities  in  the  xyz 
system  defined  in  terms  of  the  direction  of  the  applied  magnetic  field,  we 
use  the  transformation  matrix  T: 


Ha '  *  T  •  Ha  (A. 13a) 

Na.  ,  T  .  Na  .  -r1  (A.  13b) 

where 

Tn  »  cos0  cos<$  cos?  -  sin$  sin? 

^21  *  ~sin<>  cos?  ~  cosQ  cos$  sin? 

T3 ^  *  sinQ  cosi> 

T^2  ■  cos6  sin$  cos?  +  cos$  sin? 

T22  *  cosifi  cos?  -  cos6  sin$  sin?  (A.  13c) 

T32  *  sinQ  sin$ 

T 1 3  *  -  sin6  cos? 

T23  *  sine  sin? 

T33  »  cos9 


and  8,b,?  are  the  Euler  angles  relating  the  two  coordinate  systems. 
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IV.  Wal ke r 1 s  Equation  and  Magnetostatics 

Maxwell's  equations  In  the  magnetostatic  limit  are 

V  x  h  »  0  (A. 14a) 

V  •  (m  +  h)  *  0  (A. 14b) 

Using  (A. 9)  and  introducing  a  magnetostatic  scalar  potential  allows  us  to 
combine  (A. 14)  into  a  single  equation: 

3i[^ij  +  Xij]3^  "  0  (A. 15) 

where  h  ■  -  Vty,  =  3/3xi,  and  sums  over  repeated  indices  are  implied. 
Expanding  this  equation  gives 

{3x[(l+Xxx)3x  +  Xxy3y]  +  3y[(l+Xyy)3y  +  Xyx3x]  +  3Z2}*  “  0  (A.16) 

m 

If  x  is  not  a  function  of  position  this  can  be  written 

a 

2N 

[<l  +  Xxx)3x2  +  (l  +  Xyy)3y2  +  3Z2  +  ~~  3x3y  I*  *  0  (A.17) 

a  a 

where  we  have  used  the  fact  that  NXy  ■  N^.  Note  that  the  presence  of  the 

cross  term  in  (A.17)  significantly  complicates  the  solution  of  this 

equation.  All  of  the  correct  existing  calculations  involving  anisotropy  of 

which  the  author  is  aware  have  been  done  in  selected  coordinate  systems 
a 

in  which  NXy  vanished.  Unfortunately,  there  are  many  coordinate  systems  of 
practical  interest  in  which  this  is  not  the  case. 

Because  of  the  computational  difficulties  involved  with  anisotropy  we 
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will  not  consider  it  further.  In  the  absence  of  anisotropy,  (A. 17)  reduces 
to  the  Walker  equation  (69) 

[(l+x)Ox2  +  3y2)  +  3z2]*  ”  0  (A. 18) 

This  result,  combined  with  the  requirement  that  tangential  h  and  normal  b 
be  continuous  across  all  boundaries,  completes  the  formulation  of  the 
nagne cos ta tic  boundary  value  problem. 

V.  Dispersion  Relations  for  Magnetostatic  Waves  in  Ferrite  Slabs 

The  following  calculations  neglect  both  exchange  and  anisotropy. 

A.  Surface  Waves 

Consider  Che  geometry  of  Figure  A.l.  In  air,  Walker's  equation  reduces 
to  Laplace's  equation  so  the  potential  in  the  three  regions  can  be 
expressed 

-kx  +■  i  vky 

i|r  *  A  e  (A.  1 9a ) 


TT  ,  -kx  kx  ,  ivky 

'“-[Be  +  C  e  ]  e 

(A. 19b) 

kx  +  ivky 

,H1  •  d  e 

(A. 19c) 

where  v  *  ±1,  and  k  is  taken  to  be  positive  definite. 

Continuity  of  tangential  h  requires  the  potential  to  be  continuous 
across  the  surfaces  of  the  ferrite  giving  rise  to  the  equations 


-kd/2 

-kd/2 

kd/2 

A  e 

*  B  e 

+  C  e 

(A. 20) 

-kd/2 

kd/2 

-kd/2 

D  e 

=  8  e 

+■  C  e 

( A .  2  1 ) 
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Figure  A.i.  Surface  and  backward  volume  wave  geometry. 
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Similarly,  continuity  of  normal  b  gives 

-kd/2  kd/2  -kd/2  -kd/2  kd/2 

-  Ae  ”  (l+x)(Ce  -  Be  )  +  vtc(Be  +  Ce  )  (A. 22) 

-kd/2  -kd/2  kd/2  kd/2  -kd/2 

De  **  (l+x)(Ce  -  Be  )  +  v<(Be  +  Ce  )  (A. 23) 

The  dispersion  relation  is  obtained  by  setting  the  determinant  of  the 
coefficient  matrix  of  the  simultaneous  homogeneous  equations  (A.20)-(A.23) 
equal  to  zero.  The  result  is 

,  -2  I k I d 

n2  -  Z(Z+1)  +  ( l-e  )/4  (A. 24a) 

or  Ikl  -  -  ~  tn[4Z(Z+l)  -  4a2  +  l]  (A. 24b) 

The  absolute  values  emphasize  the  fact  that  although  surface  waves  exhibit 
field  displacement  nonreciprocity,  the  dispersion  relation  is  reciprocal 
for  propagation  in  Che  ±  v  directions.  Equation  (A. 24)  is  plotted 
qualitatively  in  Figure  A. 2. 

B.  Backward  Volume  Waves 

The  geometry  for  backward  volume  waves  is  also  given  by  Figure  A. 1 , 
except  chat  propagation  is  assumed  along  the  i  z  directions  rather  than 
along  i  v.  Let  us  first  consider  even  modes  (ie.,  modes  in  which  ij/(x)  is 
even).  The  potential  in  the  three  regions  can  be  written 


I  -kx  +•  i  vkz 

=  A  e 


(A. 25a) 
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II  ivkz 

•  B  cosk^x  e 

(A. 25b) 

III  kx  +  ivkz 

i(i  •  C  e 

(A. 25c) 

Substituting 

(A. 25b)  into  Walker's  equation  (A. 18)  gives  the 

following 

relationship 

between  k^  and  k: 

k 

* 

/  -( i+x) 

(A. 26) 

The  radical 

is  real  since  1+x  <  0  in  the  volume  wave  manifold 

• 

Requiring  \j/  to  be  continuous  at  x  *  ±  d/2  gives 

-kd/2 

A  e  -  B  cos(kxd/2) 

(A. 27) 

-kd/2 

C  e  *  B  cos(kxd/2) 

(A. 28) 

while  Che  boundary  condition  on  normal  b  gives 

-kd/2 

-k  A  e  •  -B  kx  (1+x)  sin(kxd/2) 

(A. 29) 

-kd/2 

k  C  e  -  B  (1+x)  sin(kxd/2) 

(A. 30) 

The  dispersion  relation  can  be  obtained  by  combining  either  (A. 27)  and 
(A. 29)  or  (A. 28)  and  (A. 30).  The  result  for  even  modes 

cot(kxd/2)  *  -  /  -(1+x)  (A. 31) 

Uote  Chat  although  there  is  only  one  surface  wave  mode,  there  are  a 
multiplicity  of  backward  volume  wave  modes  given  by  the  roots  of  (A. 31). 
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1  _ _ 

2cot(kxd)  »  •  -  /-( 1+x)  (A. 35a) 

Ml+x) 

This  is  the  form  first  given  by  Damon  and  Eshbach.  Equations  (A. 31)  and 
(A. 33)  can  be  combined  in  an  even  simpler  way,  however,  by  using  the 
identity  tan(9+it/2)  »  -cot9.  The  result  is 

t  .  1 

can  [  (kxd  +  nir ) / 2  ]  »  -  -  ■  ■  -  ;  n  »  0,1,2,  .  .  (A.  35b) 

'“<1+X> 

Clearly  for  all  even  integers  this  is  equivalent  to  (A. 31)  while  for  odd 
integers  it  gives  (A. 33).  The  backward  volume  wave  spectrum  as  given  bv 
*  \*35)  is  illustrated  qualitatively  in  Figure  A. 3. 
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C.  Forward  Volume  Waves 

Forward  volume  waves  propagate  in  a  normally  magnetized  slab  as 
illustrated  in  Figure  A. 4.  In  the  absence  of  anisotropy,  the  propagation 
of  forward  volume  waves  is  isotropic  in  the  plane  of  the  slab.  We 
therefore  arbitrarily  choose  propagation  in  the  t  y  directions. 

For  even  modes,  the  potential  in  the  three  regions  is  of  the  form 


I  -kz  +  i vky 

<(i  •  A  e 

(A. 36a) 

11  i vky 

<(i  *  B  cos(kzz )  e 

(A. 36b) 

111  kz  +  i vky 

iji  =  C  e 

(A. 36c) 

Requiring  to  be  continuous  at  z  -  t  d/2  gives 

-kd/2 

A  e  »  B  cos(kzd/2) 

(A. 37 ) 

CJ 

C  e  *  B  cos(kzd/2) 

(A. 38) 

and  matching  normal  b  at  these  boundaries  gives 


-kd/2 

-kAe  =»  -B  kz  sin(kzd/2)  (A. 39) 

-kd/2 

kCe  =  B  kz  sin(kzd/2)  (A. 40) 

Combining  either  (A. 37)  and  (A. 39)  or  (A. 38)  and  (A. 40)  gives  Che 
dispersion  relation  for  even  modes 

1 

tan(kzd/2)  * 


/-O+x) 


( A. 4  l ) 


z 


Figure  A. 4.  Forward  volume  wave  geometry. 


where  from  Walker's  equation 

k2  3  k  /-( 1+x)  •  (A. 42) 

As  before,  the  odd  modes  can  be  obtained  by  redefining 

4.11  -  B  sin(k2z)  eivky  (A. 43) 

The  resul t  is 

1 

-cot(kzd/2)  3  —  -  (A. 44) 

/  -(l+x) 

As  in  the  backward  volume  wave  case,  the  even  and  odd  modes  can  be  combined 
in  a  single  equation  using  tan((H-ir/2)  3  -cot9  to  obtain 

tanfkj.d/2  +  nn/2  ]  3  — (A. 45) 

/  -d+x) 

The  magnetostatic  forward  volume  wave  spectrum  is  illustrated  in  Figure 

A  •  5. 

D.  Effects  of  an  Adjacent  Parallel  Ground  Plane 

Some  physical  insight  into  the  effects  of  placing  metal  in  contact  with 
one  face  of  a  ferrite  slab  can  be  obtained  through  the  use  of  the  image 
theorem.  As  illustrated  in  Figure  A. 6,  magnetic  dipoles  parallel  to  a 
perfect  electric  conductor  are  imaged  without  inversion,  whereas  the  images 
of  dipoles  normal  to  the  surface  are  inverted  (70). 

In  the  case  of  forward  volume  waves,  the  rf  magnetic  dipoles  are  always 
nirailei  to  the  ground  plane  and  so  they  are  imaged  without  inversion, 
hence  me  nodes  are  identical  to  the  even  nodes  of  1  slab  twice  as  thick 
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without  a  ground  plane. 

Reference  to  (A. 24),  (A. 35)  and  (A. 45)  reveals  that  for  a  fixed 
frequency,  k  ■  1/d  for  all  three  types  of  modes.  Since  the  band  edge  at 
k**0  does  not  depend  on  d,  this  implies  that  thicker  slabs  sustain  faster 
modes.  Hence  the  presence  of  a  conducting  plate  increases  the  velocities 
of  forward  volume  waves. 

In  the  case  of  surface  waves  and  backward  volume  waves,  the  rf  magnetic 
dipoles  have  components  both  perpendicular  and  parallel  to  the  metal 
surface.  As  a  result,  the  effect  on  the  wave  propagation  is  more  complex 
than  simply  doubling  the  thickness.  However,  the  qualitative  result  that 
the  presence  of  the  conducting  plate  increases  mode  velocities  is  still 
correct. 

The  case  of  a  metal  surface  placed  a  finite  distance  away  from  the  slab 
can  be  understood  in  terms  of  a  transition  between  the  limiting  cases  of  a 
conducting  surface  in  contact  with  the  slab,  and  the  conducting  surface 
completely  removed. 

Due  to  the  consequences  of  Laplace's  equation,  the  magnitude  of  the 
potential  outside  of  the  slab  always  decays  as  exp(-kL),  where  k  is  the 
wavenumber  and  L  is  the  distance  along  the  slab  normal1.  As  a  consequence, 
for  small  k  the  metal  will  appear  electrically  as  if  it  were  in  contact 
with  the  slab,  while  for  large  k  the  mode  will  not  be  affected  by  the 
metal.  Intermediate  values  of  k  will  be  characterized  by  a  smooth 
transition  between  these  two  limits. 

‘The  form  of  the  decay  away  from  the  surface  becomes  more  complicated  if 
the  bias  field  is  not  uniform  (see  Chapter  Three). 


If  d/L  <<  1  where  L  is  che  is  the  separation  between  the  slab  and  the 
metal,  only  waves  with  small  k  will  be  affected.  The  results  are 
illustrated  qualitatively  in  Figure  A.  7.  The  actual  dispersion  relations 
for  the  three  modes  are  given  below,  with  the  corresponding  geometries 
illustrated  in  Figure  A. 8. 


HSSW  (57): 


,2|k|d 


1  rl+(flv+-Z)[l+tanh(-lklL): 
2(flvfZ)+l  Ll-(nv-Z)  (l-tanh(-lklL) : 


(A. 46) 


where  <*  e-^-  ^ i  1  y  ;  v  •  ±1. 


MSBVW  (59): 


lk|d  j  2+y(  1+e-^  1  ^  ^) 

M'l+x)  2  A(  1+x) 


(A. 47) 


where  <(/  <*  eiv|klz  ;  v  «  tl. 


MSFVW  (58): 


-cot[-/-0+x)lkld 


I+X+tanh( I k I L) 

/-( 1+x) [ tanh( I  k  10+1 ' 


(A. 48) 


i  v  Ik  lx  ^ 

where  «  e  ;  v  ■  ±1 ,  x ^  ■  (x , /) . 


I  —  L  finite 


(b) 


60 


Figure  A. 7.  Qualitative  propagation  characteristics  for  ferrite- 
d ieiectric-metal  layered  structure.  (a)  Dispersion  relation  for  forward 
waves.  (b)  Delay  characteristic  for  forward  waves.  (c)  Delay 
characteristics  for  backward  waves. 
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Figure  A. 8.  Geometry  and  coordinate  systems  for  ferrite-dieLectric-metal 
layered  structure. 


E.  Effects  of  Finite  Slab  Width 


Generally  speaking,  introducing  a  finite  width  to  the  ferrite  slab  has 
the  effect  of  adding  a  transverse  component  to  the  total  k  vector.  In 
ocher  words,  we  view  the  mode  as  consisting  of  a  plane  wave  propagating  at 
an  oblique  angle  and  bouncing  back  and  forth  down  the  slab. 

As  in  Section  V.D,  the  forward  volume  waves  are  again  easier  to 
understand.  In  the  present  case  this  is  because  the  propagation  is 
isotropic  in  the  plane  of  the  slab.  Hence  the  frequency  depends  on  the 
transverse  component  of  k  (kc)  in  exactly  the  same  way  as  the  longitudinal 
component  (kj).  The  transverse  component  is  quantized  and  held  relatively 
constant  by  the  boundary  conditions  while  the  longitudinal  component 
depends  strongly  on  frequency.  Hence  kt  is  negligible  when  kj  is  large, 
but  kt  dominates  for  kj  +  0. 

The  various  quantized  values  of  kc  give  rise  to  an  infinite  number  of 
nodes  which,  when  combined  with  the  various  modes  corresponding  to 
thickness  variations,  comprise  a  mode  family  entirely  analogous  to  the 
modes  of  a  rectangular  metallic  electromagnetic  waveguide. 

The  effect  on  surface  waves  is  more  complicated,  however,  since 
propagation  in  the  plane  is  not  isotropic.  As  a  result,  when  k^  is  small 
and  kt  dominates,  the  wave  has  a  volume  wave  character,  whereas  when  kg 
dominates  the  surface  wave  character  is  rescored.  Hence  the  finite  width 
induces  a  volume  wave  band  near  the  band  edge  where  k^  0. 

To  the  author's  knowledge,  no  calculation  of  the  width  effects  on 
hackward  volume  waves  has  been  done  to  date.  However,  it  seems  very  1 ikel 
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that  a  similar  argument  to  Che  above  could  be  made  in  favor  of  the 
existence  of  width-induced  surface  waves  near  the  low-k  band  edge.  It 
would  be  expected  that  these  surface  waves  should  show  backward  wave 
character. 

For  d/w  <<  1,  where  w  is  the  slab  width,  only  the  spectrum  near  kj-*0 
will  be  affected.  The  qualitative  propagation  behavior  for  surface  and 
forward  volume  waves  for  this  case  is  illustrated  in  Figure  A. 9.  The 
presence  of  a  ground  plane  such  that  d/L  <<  I  is  also  assumed.  The 
dispersion  relations  for  these  cases  are  given  below.  It  should  be  noted 
that  both  of  these  results  are  based  on  the  somewhat  artificial  assumption 
of  spin-pinning  at  the  edges  of  the  sample  (44,65);  the  electromagnetic 
boundary  conditions  at  the  edges  have  not  been  rigorously  satisfied.  An 
integral  equation  formulation  which  rigorously  satisfies  these  edge 
boundary  conditions  has  been  described  by  Morgenthaler  (38,66). 
Unfortunately,  the  solution  to  the  integral  equation  is  not  available  in 
closed  form. 

HSSW  (44)  (width-induced  volume  waves  excluded): 

2Md  ZfH-nv|k|+(Z2-Q2)(M-N) 

ZM-(5vlkl+(Z2-a2)(M+N) 

ZM-nv|kl+(Z2-fl2)[M-Ntanh(NU  1 

x  - - - r -  (A. 49) 

Zrt+flv|kl+(Z2-!J2)[M-Wtanh(NL)  J 


where  M2  *  [  nir/w  ]2/(l+x)  +  k2  ,  N2  *  [nu/wj2  +  k2 


lowest  width 
mode 


infinite  width 
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and  \J/  *  ;  v  ■  il* 

MSFVW  (65): 

af 2  —  a^2tanh(a^L) 

COt<“fd)  ’ 

where  af2  “  -(1+x)  [k^nir/w)  2]  ,  2  *  k2  +  (nw/w)2  , 


(A. 50) 


and  i|i 


eivlkly  .  v  ,  si. 
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APPENDIX  B 

Harmonic  Oscillator  Mathematics 

A.  One-dimensional  Harmonic  Oscillator 

1.  The  equation  for  a  one-dimensional  quantum  mechanical  harmonic 
oscillator  can  be  written 

_3£4>(x)  +  r2e  -  xi]*(x)  -  0  (B.l) 

3x2  a2  a1* J 

where  a  is  a  characteristic  length  defined  in  the  context  of  a  particular 
problem.  In  natural  or  normalized  units  the  equation  assumes  the  form 

32$(0  +  (2e  -  S2H(0  »  0  (B.2) 

where  £  =  x/a. 

2.  The  eigenfunctions  of  (B.l)  which  vanish  at  infinity  are  given 
by  the  He rrai te-Gaussiar.  functions 

-xz/2o2 

Hn(x/oO  e 

4>n(x)  -  -  (B.3) 

(a  2n  n  !  /if)  *  /2 

where  Hn(5)  is  the  nc^  order  Hermite  polynomial  (see  (8.9) ) .  The 
corresponding  solutions  to  (8.2)  are  given  by 

-Sz/2 

«n<Oe 

4>nU)  -  - — -  (B.4) 

( 2n  n!  /ir  )  1/2 
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The  eigenvalues  belonging  Co  either  eigenfunction  are 

£  *  n  +  1/2  ,  n  »  0,1,2,  .  .  .  (B.5) 

3.  Hermite  polynomials  are  solutions  to  the  equation 

d^  "  2^Hn  +  2nHn  -  0  (B.6) 

d£2  d£ 


4.  The  Hermite  polynomials  can  be  expressed 


where 


Hn(0  -  fd^  e“s2+2s?l 

dsn  S“0 


F(s.O 


-s  2+2s  5 
e 


OB-7) 


(B.8) 


is  che  genera  ting  function  for  Hermite  polynomials.  The  first  few  Hermite 
polynomials  are 


H0U)-  1 

H3(0-  -I2C+8C3 

H i ( 5)"  2$ 

H^U)-  12-48^24.1 6  C 

(B.9) 

H2(C)-  -2+4^2 

H5(0-  1205-160c2f32^5 

5.  Useful  recurrence  relations  for  Hermite  polynomials  are 


dHn  =  2nHn_1 
d£ 

Unfi¬ 


ts.  10a) 


2^»n  -  2"Hn-l 


(B. 10b) 
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6.  The  orthogonality  integral  for  Hermite  polynomials  is 


/  -  5m  2«  n!/n 


CB.ll) 


7.  The  orthonormal ity  of  the  He rmi te-Gaussian  functions  (B.3)  and 
(B.4)  can  be  easily  demonstrated  using  (B.ll): 


/  dC  ^(5)^(0  =*  |  dx  4>nCx) qm(x)  -  6, 


mn 


(B. 1 2) 


3.  The  derivative  of  the  He rmi te-Gaussian  $(x)  can  be  obtained 
using  (B.3)  and  (B.10): 


d<t>nU)  =■  if'On-l  B  ~  ^rrH  BEL 

dx  /  2  /  2 


(8-13) 


9.  A  measure  of  the  spatial  extent  of  the  functions  4>n(x)  is  given 
by  the  location  of  the  turning  point  of  equation  (B.l)  defined  as 


2e  -  xtp2/a2  =  0 
Combining  this  with  (B.5)  gives 


(B. 14a) 


x  tp 


a/  2n+ 1 


(B. 14b) 


3.  Two-dimension  Ha  rmonic  Oscillator 

1.  The  equation  of  a  two-dimensional  quantum  mechanical  harmonic 
oscillator  is 


a  2  3  2 +•  ol,  2  3  2  +  (2t  -  x 2  -  v2  )  =  0 

Tx2  '172  — - 


a*2  «v2 


(B.15) 
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Assuming  a  product  function  for  this  separates  into  two  one-dimensional 
harmonic  oscillator  equations 

iiix  +  [£f*  -  Sl  ]*  -  0  ( B .  1 6a ) 

3x2  <*x  «« 

-  0  ( B.  1 6b) 

3y^  Oy^  dy4  " 

where  e  »  +  £y 

In  the  case  of  eouation  (2. 80) (Chapter  Two),  we  have 

2ex  +  ^  ey  *  (a  +  Go  -  Z0)/Aex  (B.16c) 

^  °y2 

where  “  ( Igx/^x)  1  /l4  and  Oy  *  ( Aex/By)  1 /’1+. 

2.  The  eigenfunctions  of  C B. 15)  which  vanish  at  infinity  are 
simply  products  of  one-dimensional  Herat  te-Gaussians  as  indicated  by 
(8.16).  We  have 


'I'n  n  <x«y>  *  *n  <x>$n  <y> 
x  y  x  y 


( B. 1 7 ) 


where  $n  (x^)  is  given  by  (B.3). 


3.  The  orthonoraal  i  ty  of  the  functions  ij/n  n  is  easily 

x  y 

demonstrated  using  (B.17)  and  (B.L1).  The  result  is 


/dx/dy4inni^nm  ■  6n  m  ®n  ,ra 
■  <»  -m  xyxy  x  v  x  y 


(B.1S) 
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C.  Properties  of  the  General ized  Functions  tfrm(z) 


l.  For  the  analysis  of  Chapter  Three,  it  is  convenient  to  introduce 
the  generalized  one-dimensional  He rmi te-Gaussian  function: 

~z  2/2°tn2 

—  (B. 1 9) 


n  Vz/an>  e 

4»ra(z) 


(a„  2®  «!  /it  )  *  22 


These  functions  represent  the  set  of  m  solutions  to  n  different  harmonic 

oscillator  problems. 

n 

The  functions  $n(z)  are  shown  in  Chapter  Three  to  be  solutions  of  the 
e  qua  t ion 


32<t>n(z)  +  (A-Bz2)A2^n(z)  =  0 

3z2 


(B. 20) 


where  A2  =  B(2n+1)2/A2  ,  and 


°n 


/ 


B(2n+1) 


2.  The  functions  <^(2)  are  not,  in  general,  orthogonal: 


/  dz  (^(z)  $q(z)  *  0 


(B.21) 


This  is  simply  because  for  n  *  r  the  functions  are  solutions  to  different 
harmonic  oscillator  problems.  It  follows  that  functions  having  the  same 
characteristic  lengths  (ie.,  n“r)  are  orthonormal: 


/  *  )  $q( 2  ) 


Jmq 


(B.22) 
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3.  The  derivative  of  ^(z)  is  given  by 


n  .  n  n  _ 

3tm(z)  =  i  fV-l  /®  ~  *nH-I  /Eti] 
3z  ®n  /  2  /  2  J 


(B. 23) 


4.  The  quantity  zzqn(z)  is  needed  for  the  integral  equation  theory 

of  Chapter  Three.  This  quantity  can  be  expressed  in  terms  of  the  second 
n 

derivative  of  ^  by  use  of  the  differential  equation  (B.20a).  The  result 


_  n  n  n 

z  <t>n  *  °n  ii*n  + 

3z2  B 


(B.24) 


The  second  derivative  can  be  obtained  by  differentiating  (B.23).  This 


gives 


3 2 <|>n  “  /n(n-l')  4^-2  "  (2n+i )  $n  +  /(n+l)(n+2)  $^2 
3z2  2on2  2a,,2 


(B. 25) 


Substituting  (B.25)  and  (B.20c)  into  (B.24)  gives  the  final  result 


z2$n  ”  A/n(n-l)  1^-2  +  _A  $n  +  A/(n-H)  (n+2) 
2B(2n+l)  2B  2B(2n+l) 


( B. 26) 


5.  The  Fourier  Transform  of  bm(z)  is  defined  as 


n  -iCz  n 

VO  •  /  dz  e  ^(z) 


(B.27) 


It  is  computationally  convenient  to  define  the  normalized  variables  u’tou, 
and  ^“z/oin.  Making  these  substitutions  and  using  (B.I9)  gives 
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Control  of  the  Frequency  of  Energy  Circulation  of 
Magnetostatic  Modes  in  a  Sphere 

The  SM  thesis  of  Mr.  Daniel  A.  Fishman,  now  nearing  completion,  is  concerned 
with  the  control  of  frequency  at  which  magnetostatic  mode  energy  circulates  in 
a  YIG  sphere;  the  uniform  precession  is  of  particular  interest. 

A  portion  of  the  thesis  proposal  is  given  below: 

The  use  of  ferrimagnetic  materials,  particularly  YIG  spheres,  in  microwave 
magnetically  tunable  filters  is  widespread.  In  general  these  devices  have  good 
power  handling  character! sti cs .  It  is  proposed  that  this  is  due,  in  part,  to 
the  fast  velocity  of  energy  circulation  of  the  modes,  (analogous  to  the  group 
velocity  for  planes  waves),  and  that  this  works  to  prevent  low  power  limiting. 
However,  it  has  been  shown  by  Morgenthaler  that  this  velocity  can  be  altered  by 
changing  the  specific  electromagnetic  boundary  conditions.  The  situation  where 
this  boundary  condition  is  of  the  form  of  a  concentric  conducting  spherical  shell, 
is  studied  theoretically  and  approximated  experimentally.  The  theoretical 
analysis  indicates  that  there  is  a  critical  ratio,  between  the  radius  of  the  YIG 
sample  and  that  of  the  conducting  spherical  cavity,  where  the  velocity  will  become 
zero.  These  calculations  did  not  include  the  electromagnetic  losses  that  arise 
in  the  experimental  situation,  thus  the  velocity  approaches  zero  at  this  ratio. 

In  addition  to  the  effect  of  the  energy  velocity  upon  signal  delay  times,  there 
is  the  effect  of  decreasing  the  critical  power  threshold  for  nonlinear  effects. 

The  onset  of  instabilities  in  the  uniform  precession  of  the  magnetization  is  due 
to  a  level  of  rf  energy  density  inside  the  sample  which  exceeds  a  certain  critical 
value.  A  decrease  in  energy  velocity  could  be  expected  to  decrease  the  amount  of 
incident  power  required  to  obtain  this  critical  energy  density. 

The  object  of  this  work  is  to  observe  the  predicted  effects  the  velocity  of 
energy  circulation  may  have  on  the  onset  of  spin  wave  instabilities  and  other 
nonlinear  effects. 
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